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ABSTRACT 


The  contents  of  this  report  describe  accomplishments  during  the  first  half  of  the  Low-Level 
Critical  Air  Turbulence  (LO-LOCAT)  Phase  III  program.  Data  that  were  obtained  during  flight 
over  routes  near  McConnell  Air  Force  Base,  Kansas,  and  Edwards  Air  Force  Base,  California,  are 
presented.  Also  included  are  Lime  histories  of  some  of  the  larger  magnitude  gust  velocities 
recorded  over  the  Peterson  Field,  Colorado  route.  Approximately  72  hours  of  low-level  (0-1000 
feet)  turbulence  and  associated  meteorological  data  recorded  from  August  16,  1968,  through 
January  8,  1969,  are  presented.  The  turbulence  environment  is  analyzed  in  terms  of  gust  velocity 
primary  peaks,  amplitude  samples,  rms  values,  maximum  values  and  derived  equivalent  gusts. 
Mathematical  expressions  for  turbulence  spectra  and  scale  length  statistics,  and  correlations 
between  atmospheric  gust  velocities  and  meteorological  and  geophysical  phenomena  are  shown. 
The  most  predominant  characteristics  of  these  data  are  the  increases  in  rms  values,  scale  lengths, 
Taylor  micro  scales  and  maximum  gust  values  as  compared  to  the  LO-LOCAT  Phases  1  and  II 
data.  These  differences  are  attributed  to  the  increased  speed  of  the  Phase  III  airplane  which  is 
providing  a  better  definition  of  the  longer  wavelengths  of  atmospheric  turbulence.  This  report 
consists  of  two  volumes. Volume  I  contains  a  discussion  of  data  acquisition  and  quality,  along 
with  a  preliminary  analysis  of  turbulence  and  meteorological  data.  Instrumentation  details, 
calibrations  and  checks,  data  processing,  gust  velocity  data  (power  spectra,  peak  and  amplitude 
count,  etc.)  and  test  logs  are  presented  in  Volume  II. 

(Distribution  of  this  abstract  is  uniimited) 
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Symbols  and  Abbreviations : 

A  Regression  coefficient . 

A  thru  H.  J,  K  Coefficients  in  the  equation  for  angle  of  attack. 

A'  thru  H  ,  J ,  K'  Coefficients  in  the  equation  for  angle  of  sideslip. 

A(n)  Complex  finite  transform. 

A(n)  Complex  conjugate  of  A(n). 

a  “Universal  Constant"  in  the  longitudinal  gust  velocity  component  spectrum 

expression,  constant  in  a  equation;  and  shape  parameter  used  in  spectra 
mathematical  expressions. 

B  Air  stability  ratio. 

b  Constant  in  equation,  also  a  constant  in  gust  velocity  rms  distribution  equation. 

br.  Real  part  of  Fourier  transform. 

C  Constant  in  the  Lumley-Panofsky  equation 

c  Mean  aerodynamic  chord, 

c.g.  Airplane  center  of  gravity. 

Cr  Airplane  lift  curve  slope  in  I  radians, 

cpf  Cycles  per  foot, 

cps  Cycles  per  second. 

D  Distance  from  airplane  c.g  to  measurement  location  in  feet, 

d  Degrees  of  freedom,  distance  traveled  in  statute  miles,  differential, 

db  Decibel 

Exp  Actual  energy  per  tv.it  area  radiated  from  a  real  surface . 

Ebb  Energy  per  unit  area  radiated  from  a  black  body. 
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NOTATIONS 


Eb 

Ea 

c 

ran 

f 

fb 

f. 

fs 

fps 

f. 


Gf 

G- 

g 

H 

H; 

H(0 

Hg 

H  i  .H  j 
HN(f) 
Hp 
h 


Incoming  radiant  energy. 

Total  crossings  of  the  level  a 

F-ratio  or  F-test  of  significance  in  regression  analysis. 

Fourier  transform  modulus. 

Frequency  in  cps;  function. 

Frequency  of  occurrence  in  band  b. 

The  highest  frequency  of  a  low-pass  filter  whose  amplitude  is  passed  with  unity 
gain  (cutoff  frequency)  in  cps. 

Nyquist  or  folding  frequency  in  cps. 

Feet  per  second 

Sampling  frequency  in  samples  per  second. 

A  significant  frequency  present  in  the  spectrum  in  cps. 

Gust  factor. 

Ground  speed  in  feet  per  second 

Earth's  gravitational  constant  at  sea  level  =  32. 1  74  ft. /’sec.2 
Altitude  above  the  earth's  surface  in  feet,  true  altitude  in  feet 
Calibrated  pressure  altitude  in  feet 
Transfer  function  of  the  filter. 

Symbol  for  the  element  Mercury 
Inc*  .ated  pressure  altitude  in  feet. 

Constrained  transfer  function. 

Pressure  altitude  in  feet . 

True  heading  in  degrees. 
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h. 

hn 

K 

h(t,r) 

ho 

i 

ip* 

IAS 

IRIG 

in. 

J 

K*»*> 

KIAS 

K« 

K* 

k 

L 

Lf 

Lh 

1 

M 


Magnetic  heading  in  degrees. 

Filter  weights  at  time  t  *  n  At . 

Constrained  filter  weights. 

Time  domain  weighting  function. 

Constant  in  empirical  scale  length  equation. 

Counter,  time  series  simple.  _ 

Inches  per  second. 

Indicated  airspeed. 

Inter-range  Instrumentation  Group 
Inches . 

Square  root  of  minus  one. 

Constants. 

Knots  indicated  airspeed. 

Gust  alleviation  factor. 

Ram  recovery  factor  for  OAT  probe 

Spatial  frequency  in  cpf;  also  defines  order  of  the  derivative  of  the  transfer 
function  H(0. 

Truncated  standard  deviation  frequency  limits  corresponding  to  0.04,  0.667  and 
10  cps,  respectively,  In  cycles  per  foot. 

Turbulence  scale  length  in  feet. 

Minimum  focal  range  of  the  camera  lens  in  feet. 

Constant  In  empirical  scale  length  equation. 

Distance  from  probe  tip  to  accelerometers  in  feet. 

True  Mach  Number. 


NOTATIONS 


M. 

m 


mb 

mm 

N 

NBFM 

No 

Ncl 
N  p 


N 


p* 


N 


a 


n 


OAT 

Or, 

O', 

PSD 

P. 

P. 

Po 

Pa 

0 

q= 


Magnitude  of  the  center  of  the  step  preiaure. 

Maximum  lag  -  f„  /  A  f,  slope  of  a  line,  number  of  function*  of  the  input 
variables,  and  meters, 

Millibar 

Millimeter. 

Total  number,  condition*,  data  point*,  samples,  peaks,  bands. 

Narrow  band  frequency  modulation, 

Characteriitic  frequency  from  power  »pectral  density  ■  cpf. 

Number  of  crossings  of  the  zero  level. 

Total  number  of  peaks  per  mile. 

(characteristic  frequency). 

Total  number  of  peaks  obtained  using  extrapolation  technique, 

Crossings  per  miie  of  the  level  a. 

Acceleration  in  g  units,  a!:o  a  counter,  number  of  data  pointa,  shape  parameter 
tn  spectra  mathematical  expressions, 

Outside  air  temperature  in  eF. 

Run  test  observation  of  fhe  mean. 

Run  test  observation  of  the  mean  square 

Power  spectral  density. 

Static  pressure  in  inches  of  mercury. 

Indicated  static  pressure  in  inches  of  mercury. 

Static  pressure  at  sea  level  under  standard  conditions  -  29.921  Inches  of  mercury. 
Levtl  crossing  probability 

Number  of  independent  variables  in  the  regression  equation. 

Calibrated  impact  pressure  in  inches  of  mercur> . 
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i 


<1 

Ra 

R 

C*] 

Ri 

rm» 

R. 

R(  r  ) 

R(  r  )/R(o) 
S 

SCO 

St 

S  31 
T 

T. 

Tab 

Tbb 

Temp. 

To 

Tj 

Tl 

To 


Indicated  impact  preuure  in  inchei  of  mercury. 

Radar  altitude,  feet 
Richardson  number. 

Rotation  matrix  uaed  to  tranaform  measurements  from  the  airplane  reference 
axia  to  the  earth  reference  axei. 

Value  of  leaat  square  line  at  a  given  wind  speed. 

Root  mean  square. 

Distance  from  c.g.  to  camera  lens  in  feet. 

Gust  velocity  covariance  function. 

Gust  velocity  autocorrelation  function  where  r  a  o  •  •  •  rL . 

Airplane  planform  wing  area  in  ft.2,  horizontal  distance  in  miles. 

Subcarrier  oscillator 
True  solar  time  of  day  ■  hours, 

Sine  of  solar  elevation. 

Period  of  time  function  in  seconds. 

Ambient  air  temperature  in  degrees  Rankine 
Actual  temperature  of  real  surface  in  degrees  Rankine. 

Temperature  of  black  body  in  degrees  Rankine. 

Temperature. 

Ground  surface  temperature  in  degrees  Fahrenheit. 

Calibrated  outside  air  temperature  in  degrees  Rankine. 

Lower  time  limit  of  weighting  function  in  seconds. 

Temperature  at  sea  level  under  standard  conditions  (518,69  degrees  Rankine). 


xxiv 


NOTATIONS 


Surface  temperature  as  recorded  by  radiometer. 

Upper  time  limit  of  weighting  function  in  seconds. 

Time  in  sec.,  alto  standardized  variable. 

Derived  equivalent  gust  velocity  in  fpa. 

Probe  motion  terms 

Longitudinal  gust  velocity  in  fps,  positive  aft. 

True  airspeed  or  ground  speed  In  fpt. 

Equivalent  airspeed  in  fps 
Confetti  fall  rate  in  fps. 

Ground  speed  in  fps 
Indicated  airspeed, 

Confetti  fall  rate  at  sea  level  in  fps. 

True  airspeed  in  fps. 

Represents  the  matrix  of  true  airspeed  components  and  corrections  for  pitch 
and  yiw  rate 

Lateral  gust  velocity  in  fpa,  positive  to  the  right. 

Airplane  weight  in  pounds,  wind  speed  In  fps. 

Wind  angle,  angle  In  degrees  between  airplane  ground 
track  and  wind  vector. 

Wind  direction  in  degrees  azimuth. 

Average  easterly  winds  in  fps. 


Average  northerly  winds  in  fps 

Vertical  gust  velocity  in  fps,  positive  upward. 


NOTATIONS 


Amplitude. 

Mid-band  value. 

i  th  value  of  x,  predictor*,  coefficient  of  limple  Linear  correlation,  and 
value  of  gust  velocity  rraa  at  a  given  wind  speed. 

Amplitude  values  of  the  time  series. 

k,h  record  included  in  the  ensemble  averaging  scheme. 

Smoothed  values. 

Sampled  value  of  time  series. 

Filter  output  data  at  equi-epaced  intervals. 

Sampled  value  of  low-paw  filtered  time  series. 

Sampled  value  of  high-pats  filtered  time  series. 

Level  of  gust  velocity  •  fps, 

Dependent  variable,  regression  function 

Intercept  of  y  axis,  random  variable. 

Input  variable. 

Percent. 

Angle  of  attack  in  degrees,  levels  of  gust  velocity  used  in  level 
crossing  procedure;  confidence  limit. 

Angle  of  sideslip  in  degrees. 

Lapse  rite  in  degreei  Fahrenheit/ft.  (+  indicates  tempersture  decrease 
with  increasing  altitude) 

Dry  adiabatic  lapse  rste. 

Ratio  of  specific  heats  for  air  ■  1 .40. 

Coherency  function. 


Incremental  or  difference. 


NOTATIONS 


AT 

Ah 

AH 


Aps 

APa/q. 

APOI1 

SP 

AP  a 

*PS 

At 

A  W/A  H 
AWdh 

A  WH 

A  « 

A°rR 

A«PL 

A  a  PR 

A  p 
Aflrp 


A  change  in  frequency  in  cps,  also  the  width  between  where  the  filter  tniufer 
function  equals  one  and  zero. 

Static  pressure  error  in  feet. 

Distance  above  or  below  reference  point  HIT  ,  also  difference  in  terrain 
elevation-in  feet. 

Incremental  load  factor  in  acceleration  units. 

Static  pressure  ereor  in  inches  of  mercury. 

Pressure  coefficient  •  dimensionless 

Static  pressure  position  error  correction  inches  of  Mercury. 

Angle  of  attack  differential  pressure  in  inches  of  mercury. 

Angle  of  sideslip  differential  pressure  in  inches  of  mercury. 

Time  interval  between  data  samples  In  seconds. 

Vertical  wind  gradient  in  fps/ft. 

Number  of  degrees  of  wind  direction  change  between 
start  and  end  of  data  sample. 

Horizontal  wind  gradient  in  ft. /sec. /ml. 

Total  angle  of  attack  correction  in  degrees. 

Angle  of  attack  correction  due  to  confetti  fall  rate  In  degrees. 

Angle  of  attack  correction  due  to  measurement  location  in  degrees. 

Angle  of  attack  correction  due  to  aircraft  rate  of  pitch  in  degrees. 

Total  angle  of  sideslip  correction  in  degrees. 

Angle  of  sideslip  correction  due  to  confetti  fall  rate  In  degrees. 
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NOTATIONS 


A#  PL 

A* 

i 


n 

4 

t 

X 

M 

*1  <*a 

*,(t) 


*o 


9 


D 


W 


9 


L 

a 


Angle  of  sideslip  correction  due  to  measurement  location  in  degrees. 

Angle  of  sideslip  correction  due  to  aircraft  ratv  oi  .  m  degrees. 

Error  in  angle  of  yaw, 

Drift  angle  in  degrees, 

Vlacous  dissipation  rate  in  ft2/sec3  filter  constant,  and  emiaaivity  factor, 
Emisaivlty  of  surface  used  for  calibration. 

Emiaaivity  of  a  real  surface 
Kolmogorov  microscale. 

Pitch  angle  in  degrees. 

Rate  of  pitch  in  degrees/sec. 

Taylor  turbulence  microscale  length. 

Air  viscosity  in  lb.  sec. /ft. also  mean  value,  airplane  mass  ratio. 

Constant  used  in  gust  velocity  rms distribution  equation. 

Ensemble  average  time  function 
Kinematic  viscosity,  degrees  of  freedom. 

3  1416 

Air  density  in  ilugs/ft? 

Air  density  at  sea  level  under  standard  conditions  •  ,002378  slugs/  feet3. 

Standard  deviation  of  a  statistical  sample;  standard  deviation  in  fps  from  gust 
velocity  spectra  between  0,0416  and  10  cps;  air  density  ratio;  Stefan-Boltzman 
constant. 

Standard  deviation  about  a  least  square  line. 

Standard  deviation  of  level  crossing  distribution. 

Standard  deviation  of  noise. 
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9 


t 


a 


U 

O' 


9 


1 


9i.0i 


T 


Tl 

<*> 

X 

i 


w 


w 


0 


Standard  deviation  of  vertical  acceleration. 

Standard  deviation  of  primary  peaks  obtained  using  extrapolated  value  of  Hpt  . 
Standard  deviation  of  recorded  data. 

Standard  dev.ition  in  fps  obtained  from  the  gust  velocity  spectra  between  0.667 
and  10  cps 

Standard  deviation  of  the  gust  velocity  time  aeries  In  fps. 

Minimum  valid  « t 
Standard  deviation  of  Y 
Standard  deviation  of  y. 

Standard  deviation  of  angle  of  attack  differential  pressure  in  inches  of  Hg. 
Standard  deviation  of  sideslip  differential  pressure  in  inches  of  Hg. 

Standard  deviation  (dispersion  of  the  distribution)  of  gust  velocity  rms. 

Corrected  standard  deviation  of  gust  velocity  rms. 

Standard  deviation  of  the  time  derivative  of  the  associated  function. 

Constant  used  in  rms  gust  velocity  distribution  equation. 

Lag  time  for  the  weighting  operation  convolution  in  seconds. 

Maximum  correlation  distance  where  rL  ■  rr.V, 

One-dimensional  gust  velocity  power  spectral  density. 

Roll  angle  ii,  degioes. 

Chi  statistic. 

Yaw  angle  in  degrees. 

Rate  of  yaw  in  degrees  per  second. 

Frequency  in  radians  per  second. 

Break  frequency  radians  per  second. 
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NOTATIONS 


SUBSCRIPTS  AND  SUPERSCRIPTS: 


b 

Corr. 

D 

e 

fUt. 

H 

i 

N 

n 

K 

max 

min 

P 

R 

T 

t 

u 

v 

w 

X 

y 

z 


Band  number. 

Corrected  value. 

Baaed  on  the  Dryden  equation. 

Extrapolated  value. 

Filtered  value. 

Horizontal. 

Sample  number. 

Nolle. 

Counter. 

Baaed  on  the  von  Karman  equation. 

Indicates  maximum  value. 

Indicatei  minimum  value. 

Peak  count  or,  based  on  the  Lumley-  Panofsky  equation. 

Recorded  data. 

Truncated,  obtained  from  wind  tunnel  instrumentation,  obtained  from  tower. 
From  the  guat  velocity  time  series. 

From  the  longitudinal  gust  velocity  comoonent. 

From  the  lateral  gust  velocity  component. 

From  the  vertical  gust  velocity  component. 

Longitudinal, 

Lateral. 

Vertical. 
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NOTATIONS 


SUBSCRIPTS  AND  SUPERSCRIPTS  (CONT.): 

o  Initial  value. 

\  Pertaining  to  quadrant  1 

II  Pertaining  to  quadrant  2 

III  Pertaining  to  quadrant  3 

IV  Pertaining  to  quad  P  it  4 

a  Related  to  level  crossing  count;  also  refers  to  significant  frequency  content 

present  in  the  data. 

Primes  indicate  Hanned  estimates  unless  otherwise  noted. 

Overbars  depict  time  means. 

2  Refers  to  starting  point  of  leg 

3  Refers  to  ending  point  of  leg 

32  Refers  to  difference  between  start  and  end. 


SECTION  I 


INTRODUCTION 


A  mode!  for  low  altitude  atmospheric  turbulence  is  being  developed  by  The  Boemg  Company  in 
cooperation  with  the  United  States  Air  Force.  This  is  being  accomplished  through  a  research 
program  designated  LO-LOCAT  (Low  Altitude  Critical  Air  Turbulence).  LO-LCCAT  is  part  cf  an 
extensive  program  known  as  ALLCAT  (Reference  1.1). 

The  purpose  of  the  LO-LOCAT  Program  is  to  determine  the  turbulence  environment  below  1,000 
feet  above  the  ground  utilizing  statistically  representative  samples  of  turbulence  data  obtained 
over  a  wide  range  of  meteorological,  topographical,  seasonal,  and  time-of-day  conditions.  The 
data  will  be  used  to  improve  turbulence  design  criteria  for  aircraft  that  are  required  to  operate 
below  1000  feet  above  the  terrain  for  extended  periods  of  time. 

The  LO-LOCAT  Program  consists  of  three  phases  of  testing  Phases  I  and  II,  which  have  been 
completed  and  are  discussed  in  Reference  1.2,  used  four  C-I31B  aircraft  as  instrumentation 
platforms  tu  obtain  the  turbulence  data.  Turbulence  wavelengths  up  to  7,000  feet  per  cycle  were 
measured  during  that  program.  The  purpose  of  Phase  III  is  to  extend  the  statistical  definition  of 
the  turbulence  environment  and  define  wavelengths  up  to  14,000  feet  per  cycle  through  the 
utilization  of  a  higher-speed  (T-33A)  airplane  as  the  instrumentation  platform. 

The  T-33A  aircraft  was  instrumented  to  measure  both  meteorological  and  gust  velocity  data.  A 
radar  altimeter,  a  Doppler  radar  system,  a  radiometer,  and  an  outside  air  temperature  (OAT) 
probe  comprise  the  basic  instrumentation  used  to  obtain  the  meteorological  data.  Atmospheric 
turbulence  in  the  0.04  to  10  cps  frequency  range  is  sensed  by  instrumentation  contained  in  a 
Boeing-designed  gust  probe  mounted  on  the  airplane.  The  gust  probe  head  being  used  is  the  one 
that  was  installed  on  a  B-5  2H  probe  during  the  testing  discussed  in  References  1.3  through  1.5. 
This  probe  head  has  been  both  statically  and  dynamically  calibrated  in  a  wind  tunnel  at  Mach 
numbers  including  the  range  of  the  Phase  III  recording  requirements.  The  head  was  mated  with 
the  body  of  a  gust  probe  used  during  Phases  1  and  II.  The  boom,  on  which  the  gust  probe  is 
mounted,  was  designed  fn  optimization  between  air  flow  and  vibration  characteristics. 

Airplane  attitude  and  roll  rate  instrumentation  are  installed  at  the  base  of  the  boorn. 
Instrumentation  contained  inside  the  gust  probe  senses  the  following: 

•  Differential  pressures  due  to 

sideslip  and  attack  angles 

•  Vertical,  lateral,  and  longitudinal  accelerations 

•  Static  and  impact  pressures 

•  Pitch  and  yaw  rates 

A  Narrow  Band  Frequency  Modulated  (NBFM)  system  is  used  to  record  the  data.  Duung  a 
one-year  ’’ight  period,  the  T-33A  is  being  flown  over  test  routes  located  in  New  York,  Kansas, 
Col-  md  California.  These  routes  were  established  to  give  a  wide  range  of  possible 
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topographical  and  climatological  conditions,  and  were  selected  such  that  populated  areas  would 
be  avoided.  The  routes  are  situated  in  the  same  vicinity  of  those  routes  used  during  PhasesI  and 
II.  Each  route  consists  of  eight  straight  legs.  The  legs  are  approximately  30  nautical  miles  in 
length.  They  are  traversed  in  the  same  direction  on  each  flight.  Prior  to  flying  each  leg, 
meteorological  survey  conditions  are  flown  at  100  and  1,000  feet  above  the  terrain.  These 
conditions  are  performed  over  the  start  point  of  each  test  leg.  Normally,  three  missions  are 
scheduled  every  other  day;  one  at  dawn,  one  at  mid-morning,  and  one  at  mid-afternoon.  This  is 
varied  as  necessary  when  weather  conditions,  aircraft  maintenance  problems,  or  other  factors 
interfere.  Two  specific  altitudes,  250  feet  and  750  feet  above  the  terrain,  are  used  for  the  gust 
data  gathering  position  of  the  flights.  The  pilot  follows  the  terrain  contour,  as  closely  as  safety 
will  allow,  using  a  radar  altimeter  to  maintain  a  constant  absolute  altitude. 

A  turbulence  sample  consists  of  four  and  one-half  minutes  of  flight  over  a  leg.  In  the  gust  velocity 
calculations,  airplane  motion  effects  are  removed  giving  three  orthogonal,  space  oriented  gust 
velocity  components.  High-pass  numerical  filtering  of  the  gust  velocity  time  function  and  dynamic 
calibration  of  all  measurements  is  accomplished  to  ensure  valid  gust  data,  especially  at  the  lower 
frequencies. 

Frequency  of  exceedance  and  tire  probability  of  exceeding  given  levels  of  gust  velocity  are 
determined  from  amplitude  count,  peak  count,  and  level  crossing  count  data.  Gust  velocity 
standard  deviations,  determined  from  peak  count,  amplitude  samples,  and  level  crossing  data  are 
compared.  Cumulative  probability  and  probability  density  distributions  are  calculated  for  these 
data.  Their  distributions  are  being  analyzed  and  compared  to  a  normal  distribution.  Magnitudes 
of  the  gust  velocity  rms  values  are  investigated  as  a  function  of  numerous  variables. 

The  properties  of  turbulence  are  evaluated  by  power  spectral  density  (PSD)  techniques.  The  PSD’s 
are  calculated  and  plotted  for  selected  conditions.  The  spectra  are  normalized  by  the  gust  velocity 
variance  a\  and  by  L  .  The  shapes  of  the  spectra  are  then  analyzed.  Mathematical 
representation  of  the  spectra  is  being  investigated  by  comparing  the  experimental  spectra  with 
several  mathematical  expressions  as  suggested  by  von  Karman,  Dryden,  and  others.  The  integral 
scale  lengths  of  the  turbulence  are  also  being  analyzed. 

Certain  characteristics  of  selected  turbulence  samples  such  as  homogeneity,  isotropy,  and 
coherency  are  calculated  and  plotted  versus  spatial  frequency.  Correlations  between  turbulence 
and  meteorological  phenomena  are  being  studied,  and  turbulence  forecasting  techniques  are  being 
investigated. 

This  report  summarizes  the  program  accomplishments  during  the  first  half  of  the  LO-LOCAT 
Phase  III  Program.  The  data  being  analyzed  are  those  obtained  during  flight  over  the  routes  near 
McConnell  Air  Force  Base,  Kansas,  and  Edwards  Air  Force  Base,  California.  This  includes  testing 
through  8  January  1969.  Time  histories  of  some  of  the  larger  magnitude  gust  velocities  being 
recorded  over  the  Peterson  Field,  Colorado,  route  are  also  presented.  Flying  was  initiated  on  16 
August  1968  and  will  be  completed  on  1 5  August  1969. 
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SECTION  II 


DATA  ACQUISITION 


INSTRUMENTATION 

Contractual  authorization  for  the  LO-LOCAT  Phase  Ill  Program  was  received  on  15  April  1968. 
The  test  airplane.  T-33A,  AF  57-G550,  arrived  at  Boeing-Wichit?  from  Holloman  Air  Force  Base, 
New  Mexico,  on  17  April  1968. 

Airplane  modification  and  instrumentation  installation  began  on  23  April  1968.  The  modification 
consisted  of  removal  of  airplane  equipment  not  required  during  this  testing  and  relocation  of 
other  equipment  to  permit  instrumentation  installation.  Among  those  items  removed  were  the 
instructor  pilot’s  seat,  control  stick,  instrument  panel,  and  other  miscellaneous  equipment  in  the 
same  general  area.  This  provided  room  for  the  installation  of  a  magnetic  tape  recording  system, 
signal  conditioning  equipment,  Doppler  receiver  transmitter  unit,  and  other  smaller 
instrumentation  items  The  APX-25  control  was  moved  to  the  pilot’s  left  side  panel  ‘o  provide 
space  for  the  Doppier  radar  control  on  his  center  panel.  The  aircraft  batteries  were  removed  from 
the  gun  bay  area  and  relocated  to  the  engine  compartment  in  partial  compliance  with  TO 
1T-33A-663,  “Conversion  of  T-33A-5LO  Navigational  Trainer  to  AT-33A-20LO  Attack  Trainer.” 
Relocation  of  the  batteries  provided  spare  for  the  gust  boom  and  sideslip  camera  installations. 
Two  additional  inverters  were  installed  in  the  engine  compartment  to  provide  60  and  400  cycle 
instrumentation  electrical  power.  All  modification  and  installation  work  was  completed  by  8 
August  1968. 

In  order  to  increase  “visibility"  of  the  airplane  while  Hying  at  low  altitude,  conspicuity  paint  was 
applied  per  T.O  1-1-4,  "Exterior  Finishes,  Insignia,  and  Markings  Applicable  to  Aircraft  and 
Missiles,” 

The  design  of  the  instrumentation  system  was  simplified  by  using,  primarily,  equipment  which 
had  previously  been  used  on  the  C-131B  airplanes  during  LO-LOCAT  Phases  1  and  II  The  major 
area  of  design  involved  modification  and  installation  of  this  equipment.  However,  a  special  radome 
was  designed  and  fabricated  and  installed  on  the  underside  of  the  airplane  to  house  the  Doppler 
radar  antenna 

Instrumentation  equipment  locations  for  the  LO-LOCAT  Phase  111  Program  are  shown  in  Figure 
2.1.  The  radar  altimeter,  Doppler,  radiometer,  and  outside  air  temperature  (OAT)  probe  are  used 
basically  to  obtain  meteorological  data.  The  sideslip  camera  shown  in  Figure  2.1  was  used  only 
during  the  inflight  calibrations  of  the  probe.  Atmospheric  turbulence  in  the  frequency  range  from 
0.04  to  10  cycles  per  second  is  sensed  by  instrumentation  installed  in  the  probe.  Instrumentation 
contained  inside  the  gust  probe  is  shown  in  Figure  2.2.  Note  that  fine  and  coarse  measurements 
are  provided  to  improve  resolution.  To  ensure  accuracy,  the  temperature  of  each  transducer  is 
closely  regulated  to  135  degrees  Fahrenheit  with  thermostatically  controlled  heaters.  Pitch  and 
yaw  angular  rate  sensors  are  installed  in  the  probe  in  close  proximity  to  the  gust  sensors.  The 
attitude  and  roll  rate  gyros  are  installed  at  the  base  of  the  boom. 
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Design  of  the  boom  on  which  a  gust  probe  is  to  be  mounted  requires  optimization  between  air 
How  and  vibration  characteristics.  The  probe  must  be  placed  far  enough  in  front  of  the  airplane 
to  minimize  the  fuselage  in^uences  on  airflow  .  On  the  other  hand,  the  effects  of  boom  flexure  on 
probe  instruments  must  be  considered.  To  achieve  this  end,  a  tapered  wall  thickness  was 
Lnco.porated  in  the  LO-LOCAT  Loom  design.  The  ratios  of  pressure  port  distances  from  the 
fuselage  to  effective  fuselage  diameter  are  1  62  and  1.16  for  the  a  and  fi  ports  and  for  the  static 
pressure  port*.  respectively.  The  gust  boom  and  probe  installation  is  shown  in  Figure  2.3. 

The  boom  originally  used  on  C-131B,  AFS3-7822,  was  shortened  and  a  mounting  was 
fabrictted  for  installation  on  the  T-33A  test  airplane.  Since  the  probes  used  on  the  C-131B 
airplanes  during  Phases  1  and  [1  were  not  calibrated  for  speeds  in  the  range  to  be  flown  during 
Phase  11!.  a  probe  head  previously  used  on  JB-S2H,  AF  61-023  (Reference  1.3)  was  installed  on 
the  T-3?-v  This  probe  head  had  been  previously  calibrated  both  statically  and  dynamically  in  the 
o-foot  transonic  wind  tunnel  at  Cornell  Aeronautical  Laboratory,  Buffalo,  New  York.  The  B-52 
probe  head  was  then  fitted  to  the  improved  probe  body  and  instrumentation  carriage  previously 
used  on  C-13  i  B,  AF  53-7795. 

A  solid  state  Narrow  Band  Frequency  Modulation  (NBF,M)  system  is  used  for  data  recording.  This 
system  is  the  same  as  that  used  during  Phases  1  and  II.  The  Inter-Range  Instrumentation  Group 
(I RIG)  Telemetry  Standards  were  used  as  a  guideline  in  the  design  of  the  system.  The  system  is 
compatible  with  remote  base  ground  stations  meeting  TRIG  standards. 

Excluding  transducers,  the  NBFM  system  consists  of  the  following  basic  components 

•  Signal  conditioners  to  control  transducer  excitation  and  bias  and  for  switching  in 
standanzauon  signals. 

•  Subcarrier  oscillators  (SCO)  to  provide  frequency  deviations  from  center  frequency  that 
are  proportional  to  the  conditioned  transducer  output  voltages. 

•  Mixer  amplifiers  which  combine  the  outputs  of  several  SCO’s  of  difterent  frequency 
bends  into  a  mixed  or  composite  signal. 

•  An  oscillator  to  provide  a  stable  reference  for  tape  flutter  compensation.  Transducer 
power  supplies  to  provide  stable  DC  voltages  for  transoucer  excitation  and  system 
correction  signals.  A  magnetic  tape  recorder  to  record  the  analog  outputs  of  the  mixer 
amplifiers 

•  A  standardization  control  panel  to  provide  the  required  system  standardization. 

•  A  control  system  installed  at  the  pilot’s  station  to  enable  the  data  gathering  sequences 
to  be  performed  with  a  minimum  of  pilot  effort. 

•  A  time  code  generator  to  provide  elapsed  time  on  the  magnetic  tape  for  time 
correlation. 

The  maximum  ereor  of  the  NBFM  data  recording  system  is  0.67  percent  of  full  scale  with  a 
confidence  of  95  percent  (Reference  2  1 )  for  the  following  recording  equipment  less  transducers: 
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•  Airborne  data  acquistion  System 

•>  FM  playback  and  discrimination  system 

•  Analog  to  digital  converter  system 

This  tatement  was  verified  by  a  test  program  of  1  3  flights  with  two  different  acquisition  systems 
covenng  1  RIG  Bands  2  through  1  1 . 

During  remote  operations,  a  playback  system  in  a  mobile  van  is  used  to  make  "quick-look'1  checks 
of  data  recorded  during  flight.  This  system  consists  of  a  tape  recorder/reproducer,  discriminators, 
time  decoder,  oscilloscope,  frequency  counter,  digital  voltmeter,  and  a  multichannel  oscillograph 
recorder  (Figure  2.4).  This  system  is  used  after  each  flight  to  establish  validity  of  the  recorded 
data,  perform  quality  assurance  checks,  and  to  determine  maximum  airplane  accelerations.  The 
equipment  installed  in  this  van  is  also  used  during  instrumentation  prefbght  and  maintenance 
checks.  A  command  radio  is  installed  in  the  van  to  provide  communication  with  the  test  airplane. 


Instrumentation  details  are  given  in  Appendix  1 
OPERATIONS 


Airplane  and  instrumentation  checkout  and  calibration  fights  began  on  9  August  1968.  These 
flights  consisted  of  tower  flybys  for  airspeed  calibra'ion,  flights  over  Cheney  Lake,  located  in 
south  central  Kansas,  for  radiometer  and  radar  altimeter  checks,  stabilized  and  maneuvering  flight 
in  smooth  air  for  gust  probe  checks,  and  angle  of  attack  and  sideslip  calibrations  using  the  confetti 
cloud  technique  With  the  exception  of  angle  of  attack  and  sideslip  calibrations,  these  checkout 
and  calibration  flights  were  all  completed  prior  to  start  of  turbulence  research  testing  which  began 
on  16  August  1968  Details  of  these  calibrations  are  given  in  Appendix  II. 

The  angle  of  attack  and  sideslip  calibration  flights  were  not  accomplished  until  4  September  1968 
due  to  delayed  availability  of  a  C-135  support  airplane  The  0135,  AF  55-3129,  and  crew  were 
provided  by  the  Air  Force,  Aeronatucial  Systems  Division,  located  at  Wright-Patterson  AFB, 
Ohio  Thu  airplane  arrived  at  Boeing-Wichita  on  3  September  1968,  and  the  confetti  hopper  was 
installed  in  the  cargo  area.  The  following  day  the  C-135  provided  two  support  fights  to  dispense 
confetti  particles  into  the  fight  path  of  the  T-33A  test  airplane. 

Turbulence  data  were  obtained  over  the  McConnell  AFB,  Kansas,  route  during  the  penod  16 
August  to  25  September  1968.  A  total  of  348  turbulence  samples  were  recorded  during  44  fights. 
Dunng  this  same  period.  66  samples  were  not  recorded  due  to  adverse  weather.  These  66  samples 
include  complete  fights  as  well  as  individual  legs  which  were  not  fown  due  to  fog  or  cloud  cover, 
Of  the  44  fights  fown,  21  were  fown  at  250  feet  absolute  altitude  and  23  at  750  feet.  No 
attempt  was  made  to  obtain  an  equal  number  of  fights  at  each  altitude  other  than  to  schedule 
alternate  fight  altitudes  each  succeeding  flight  day 

The  airplane  was  moved  to  Edwards  Air  Force  Base,  California,  and  data  were  obtained  dunng  the 
penod  8  October  1968  to  8  January  1969.  A  total  of  637  turbulence  samples  were  recorded 
during  100  fights  over  this  route.  During  Edwards'  testing,  255  samples  were  not  recorded 
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because  of  adverse  weather.  Complete  tests  as  well  as  individual  legs  were  not  flown.  Legs  2  and  7 
were  frequently  not  flown  due  to  smoke  and  haze  cover  and  leg  4  was  frequently  omitted  due  to 
fog.  Of  the  100  flights  accomplished,  49  were  flown  at  250  feet  absolute  altitude  and  51  at  750 
feet. 

One  flight  was  made  in  the  area  of  the  BREN  Tower  near  Lathrop  Wells,  Nevada,  on  8  January 
1969  (See  Section  VI).  The  purpose  of  this  flight  was  to  obtain  turbulence  data  for  comparison 
with  that  obtained  from  the  instrumerred  tower 

Turbulence  research  flights  over  the  Peterson  Field,  Colorado,  route  were  started  on  20  January 
1969.  A  safety-chase  airplane,  T-33A,  AF  58-0546,  was  bailed  to  Boeing  by  the  Air  Force  to 
support  the  turbulence  flights  over  the  high  mountains,  This  airplane  will  be  utilized  until 
completion  of  research  over  the  Peterson  route.  It  will  then  be  returned  to  the  Air  Force.  It  is 
anticipated  that  flying  over  the  Peterson  route  will  be  completed  on  2  May  1969. 

Upon  completion  of  flying  over  the  Peterson  route,  the  research  airplane  will  be  moved  to  Tinker 
AFB,  Oklahoma.  The  airplane  will  then  be  used  for  approximately  two  weeks  to  record 
turbulence  data  in  the  vicinity  of  thunderstorms.  This  flying  will  be  performed  in  support  of 
research  being  performed  by  the  National  Severe  Storms  Laboratories  at  Norman,  Oklahoma. 

The  airplane  will  be  moved  to  Griffiss  Air  Force  Base,  New  York,  in  the  latter  part  of  May  after 
completion  of  severe  storms  research  Turbulence  research  flying  will  be  performed  over  the 
Griffiss  route  until  completion  of  the  flight  penod  for  LO-LOCAT  Pnase  III,  which  will  be 
approximately  15  August  1969.  The  airplane  will  then  be  returned  to  Boeing-Wichita  for 
instrumentation  removal  and  airplane  restoration. 

Figures  2.5  and  2.6  show  the  cumulative  number  of  flights  and  samples  recorded  versus  calendar 
time.  As  of  3 1  March  1969,  90  percent  of  the  scheduled  flights  had  been  flown  and  77.5  percent 
of  the  scheduled  samples  had  been  obtained.  The  schedule  was  based  on  8  samples  per  flight  and 
no  allowance  was  made  for  legs  not  flown  due  to  clouds, smoke,  haze,  or  fog  cover.  The  percentage 
of  scheduled  samples  recorded  was  lower  than  the  percentage  of  scheduled  flights  accomplished. 

The  actual  number  of  flights  flown  have  been  less  than  those  scheduled  for  variops  reasons.  As 
discussed  previously,  several  flights  have  been  canceled  due  to  adverse  weather.  Also,  due  to  late 
availability  of  a  C-135  support  airplane,  the  angle  of  attack  and  sideslip  calibrations  were 
performed  after  turbulence  research  testing  had  begun.  Because  of  these  calibrations,  the  T-33A 
airplane  was  not  available  for  turbulence  research  flights  for  three  days, 

Unexpected  maintenance  problems  have  also  caused  the  airplane  to  be  grounded.  For  example, 
the  airplane  was  grounded  at  Edwards  Air  Force  Base  for  1 1  days  while  awaiting  a  replacement 
engine.  Four  flights  have  been  canceled  at  Peterson  Field  due  to  maintenance  problems  causing 
unavailability  of  the  safety-chase  airplane.  One  flight  was  not  flown  due  to  preparations  for  flyby 
of  the  BREN  Tower. 
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TEST  PROCEDURE 


r^vC:0~S'-;C'- 


Pnor  to  each  mission  a  complete  check  of  the  instrumentation  is  performed.  If  the 
instrumentation  is  acceptable  for  data  recording,  the  pilot  is  then  briefed  by  the  Boeing 
meteorologist  concerning  the  weather  conditions  over  the  route  The  pilot  decides  whether  or  not 
to  proceeed  w  ith  the  mission  based  on  the  weather  conditions  over  the  majority  of  the  route.  The 
airplane  is  not  flown  in  clouds,  rain,  etc.,  to  preclude  moisture  accumulation  in  the  gust  probe. 
The  pilot  considers  any  operational  or  (lying  safety  hazards  that  might  exist  and  has  final 
authority  10  ground  the  aircraft  if  he  deems  it  necessary  . 


The  routes  consist  of  eight  straight  legs,  each  of  which  is  approximately  30  nautical  miles  in 
length  The  McConnell,  Edwards,  and  Peterson  routes  are  shown  in  Figures  2.7  through  2.9,  On 
data  gathering  flights  the  routes  are  flown  beginning  with  the  starting  point  of  leg  1  and 
continuing  to  the  end  point  of  leg  8.  Each  leg  is  traversed  in  the  same  direction  on  every  flight  and 
every  effort  is  made  to  stay  exactly  on  course  while  flying  each  leg. 


Normally,  data  flights  are  scheduled  every  other  day,  This  is  varied  as  necessary  wher  weather 
conditions,  aircraft  maintenance  problems,  or  other  factors  interfere.  On  each  flying  clay  three 
missions  are  flown,  one  at  dawn,  one  at  mid-moming,  and  one  at  mid-aftemoon.  Twc  specific 
altitudes  are  used  for  the  gust  data  gathering  portion  of  the  flights.  These  altitudes  are  250  feet 
and  750  feet  above  the  terrain.  The  pilot  follows  the  terrain  contour,  as  closely  as  safety  allows, 
using  the  radar  altimeter  to  maintain  a  constant  absolute  altitude.  Only  one  altitude  is  used  for 
any  one  flying  day  and  it  is  alternated  on  successive  flight  dates. 

The  following  procedures  are  adhered  to  as  closely  as  possible  during  data  gathering  flights; 

I  Tb^  pilot  stabilizes  the  aircraft  at  1,000  feet  altitude  above  the  terrain,  and  flies  over 
tuc  starting  point  of  the  first  leg  of  the  course.  He  then  actuates  a  sequencing  switch  to 
initiate  the  automatic  data  recording  sequence. 

2.  Atmospheric  survey  data  are  recorded  for  10  seconds. 

3  An  instrumentation  standardization  cycle  is  then  initiated  automatically.  During  the 
standardization  cycle,  the  pilot  begins  a  descending  turn  to  an  altitude  of  100  feet 
above  the  terrain. 

4.  The  pilot  stabilizes  the  aircraft  at  an  altitude  of  100  feet  above  the  terrain  on  a  course 
that  will  pass  over  the  starting  point  of  the  first  leg. 

5.  After  the  standardization  cycle  is  completed,  atmospheric  survey  data  are  again 
recorded  for  10  seconds.  These  data  are  recorded,  as  near  as  possible,  in  the  same 
vertical  air  column  as  that  for  Step  2. 

6.  The  pilot  makes  a  climbing  turn  to  the  appropriate  data  gathering  altitude  and  aligns  the 
aircraft  with  the  Tint  leg  of  the  course. 

7  When  the  aircraft  passes  over  the  starting  point  of  the  leg,  gust  data  recording  is 
initiated  by  again  actuating  the  sequencing  switch.  A  turbulence  sample,  4  mir.utes  30 
seconds  in  length,  is  recorded.  During  this  period  the  airspeed  is  maintained  at 
approximately  360  knots. 

These  procedures  are  repeated  for  each  leg. 
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DATA  CATEGORIZATION 


All  data  are  categorized  according  to  terrain  features,  alt* tude ,  atmospheric  stability,  time-of-day 
season,  and  geographic  location 

The  categories  are  divided  into  precise  elements  as  indicated  in  the  following  list: 

Terrain  Features 

1.  High  mountain 

2.  Low  mountain 

3.  Desert 

4.  Flams 

5.  Swamp  (not  used) 

6.  Water 

Altitude 

1.  250  feet 

2.  750  feet 

Atmospheric  Stability 

1.  Very  stable  r  <  a 

2.  Stable  isr<5 

3.  Neutral  5  <  r  <  6 

4.  Unstable  6  i  r 


Time-of-Day 

TAKEOFF  TIMES 

-  LOCAL  TLME  * 

30  MINUTES 

Summer 

Spring-Fall 

Winter 

1.  Dawn 

Sunrise 

Sunrise 

Sunrise 

2.  Mid-morning 

0900 

0930 

1000 

3  Mid-afternoon 

1500 

1430 

1400 

Season 

1.  Spnng  -  March  21  to  June  20 

2.  Summer  -  June  21  to  September  20 

3.  Fall  -  September  21  to  December  20 

4.  Winter  -  December  21  to  March  20 

Geographic  Location 

1.  Edwards  Air  Force  Base,  California 

2.  Griffiss  Air  Force  Base,  New  York 

3.  Peterson  Field,  Colorado 

4.  McConnell  Air  Force  Base,  Kansas 
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Data  processing  is  accomplished  primarily  by  the  use  of  digital  computing  equipment.  The 
compensation  for  inadequate  frequency  response  of  angle  of  attack  and  sideslip  measurements  is 
accomplished  with  an  analog  computer.  Analog  compensation  networks  are  spot-checked  on  the 
digital  computer  for  accuracy  verification 

The  aata  processing  equipment  provides  precise  time  correlation  through  the  use  of  linear  phase 
filters,  sample  and  hold  features  to  assure  simultaneous  sampling  and  time  correlation  of  data 
during  analog-to-digital  conversion,  reduced  tape  playback  time,  and  automatic  plotting  and 
printing.  Automatic  methods  are  employed  that  precisely  correct  the  data  for  undesirable 
transducer  frequency  response  characteristics,  transducer  2ero  drift,  and  data  system  zero  drift 
and  sensitivity  changes.  Data  processing  details  are  given  in  Appendix  111. 

The  category  for  each  turbulence  sample  is  indicated  on  every  item  (listing  or  plot)  of  data  by  a 
code  number.  This  number  consists  of  six  digits  in  the  order  of  classification  as  given  in  the 
preceding  list.  The  number  preceding  each  item  is  the  code  number  for  that  item.  The  following 
example  illustrates  this  method: 

Code  No.  422324 

4  -  Terrain  Feature,  plains 
2  -  Altitude,  750  fee; 

2  -  Atmospheric  stability,  stable 

3  -  Time-of-Day,  afternoon 
2  •  Season,  summer 

4  -  Location,  McConnell  AFB 

The  categories  are  reduced  to  a  minimum  by  pooling  the  data  for  categories  which  show  like 
characteristics.  The  number  zero  is  used  to  indicate  that  all  data  under  a  given  heading  is 
pooled  An  example  to  illustrate  the  pooling  technique  follows: 

Code  No.  1 14000 

1  -  Terrain  feature,  high  mountain 
1  -  Altitude,  250  feet 
4  -  Atmospheric  stability,  unstable 
0  -  Time  of-Day,  all  pooled 
0  -  Season,  all  pooled 
0  -  Location,  ail  pooled 

All  category  numbers,  except  stability,  are  assigned  when  the  data  request  is  prepared.  The 
number  for  stability  is  assigned  during  data  processing  after  lapse  rate  has  been  computed.  A 
description  of  this  computation  is  given  in  Appendix  111. 


Gust  Probe 

19 

Signal  Cond  Rack 

Guit  Boom 

20 

Dual  Pwt  Supply 

SideiUp  Camera 

21 

Signal  Cond  Rack-Special 

Radiometer  Optics 

22 

Time  Code  Generator 

Stable  Platform  Gyro  ASN  SO 

23 

Circuit  Breaker  Panel 

Stable  Platform  Gyro  Amp-Pwr  Supply  ASN-50 

24 

Control  Pwr  Relays 

Stable  Platform  Gyro  Controller  ASN  -SO 

25 

Group  Selector 

Stable  Platform  Gyro  Compeniator  ASN-50 

26 

AutoCal  Network  Box 

Probe  Heater  Control 

27 

AN/APN  153  Doppler  Radar  RT  Unit 

Stable  Platform  Gyro  Roll  rate  Cutolf  Sw  ASN-50 

28 

AN/APN  153  Doppler  Radar  Antenna 

Rate  Gyro  Electronics 

29 

Radomc 

Roll  Rate  Gryo 

30 

Battery 

OAT  Probe 

31 

400  A  60  Cycle  Inverters 

Radiometer  Electronics 

32 

HG9050  Radar  Altimeter  RT  Unit 

Pilots  Inin  Panel 

33 

HG9050  Radar  Altimeter  Antenna 

14  Track  Tape  System 

34 

Rudder  Position 

SCO  Auembly 

35 

Elevator  Position 

Standardisation  Control 

36 

Aileron  Position- 

Figure  2.1  Instrumentation  Equipment  Locations 
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Figure  2.2  Gust  Probe  Instrumentation 
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Figure  2.6  Test  Progress  -  Number  of  Samples 


SECTION  III 


DATA  QUALITY 


A  detailed  consideration  of  the  qualitative  aspects  of  atmospheric  data  is  necessary  if  it  is  to  be 
applied  in  an  effective  manner  The  development  of  these  data  involves  the  application  of  many 
optimizations  and  engineering  judgments.  Individual  measurements  are  characterized  by  certain 
inaccuracies,  noise  content,  and  other  'imitations,  1  he  combination  of  these  variables  in  the  many 
calculations  involved  results  in  a  composite  of  characteristics.  Results  of  special  tests  and  analyses 
are  presented  to  show  characterisrics  of  LO-LOCAT  dat3  and  to  provide  insight  into  the  reasoning 
for  approaches  being  taken  during  the  progiam. 

ANGLE  OF  A  ITACK,  S  DESLIP  AND  AIRSPEED  SENSORS 


The  following  items  were  accomplished  in  order  to  accurately  determine  a  and  P  from  the 
probe  pressure  measurements 

•  1  he  probe  was  statically  and  dynamically  calibrated  in  the  wind  tunnel  *o  determine 
sensitivity  and  frequency  response. 

•  The  wind  tunnel  static  calibration  results  were  verified  and  fuselage  generated  pressure 
field  effects  were  checked  in  flight  by  the  confetti  cloud  calibration  technique. 

•  That  portion  of  the  gust  probe  overall  frequency  response  due  to  transduction  effects 
was  determined  by  subjecting  the  probe  pressure  system  to  step  inputs. 

The  probe  was  statically  and  dynamically  calibrated  in  1964  in  the  Cornell  Aeronautical 
Laboratory's  eight-foot  transonic  wind  tunnel.  The  static  wind  tunnel  calibration  was  performed 
to  calibrate  the  probe,  to  venfy  design  criteria,  and  to  determine  Mach  number  and  Reyno'd's 
number  effects  on  the  pressure  sensing  elements  over  wide  angle  displacements  and  large  speed 
ranges.  The  results  of  the  static  wind  tunnel  calibration  ate  given  in  Reference  1.3  and  Appendix 
11. 


A  dynamic  wind  tunnel  calibration  was  performed  to  determine  frequency  response  of  the  angles 
of  attack  (  a  )  and  sideslip  (  fi  ).  The  tests  were  conducted  over  a  Mach  number  range  from  0  50 
to  0  95  and  for  frequencies  from  2  to  1 2  cps.  It  was  found  that  measured  gust  velocity  frequency 
characteristics  were  aflected  by  probe  aerodynamics  and  transduction  (tubing  and  transducer) 
eftects.  These  results  art  also  given  in  Reference  1 .3  and  Appendix  II 

Ground  testing  was  accomplished  with  the  gust  probe  installed  on  the  airplane  (a  B-52H  airplane, 
Reference  1.3)  to  isolate  the  aerodynamic  effects  from  the  transduction  effects.  Results  from  the 
wind  tunnel  dynamic  calibration  provided  the  overall  frequency  response,  i.c .,  the  response  due  to 
combined  aerodynamics  and  transduction  effects.  Step  pressure  inputs  were  applied  to  the  probe 
pressure  ports  to  determine  the  frequency  response  characteristics  of  the  tubing  and  transducers. 
The  difference  between  the  frequency  response  functions  obtained  from  the  wind  tunnel  and  the 
ground  test  established  that  portion  of  the  overall  frequency  response  attributable  to  aerodynamic 
effects 
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Since  the  probe  head  being  used  during  LO-LOCAT  Phase  III  had  thus  been  previously  calibrated 
statically  and  dynamically,  it  was  not  necessary  to  repeat  this  procedure.  The  probe  shell,  which 
was  used  on  the  C-131B  aircraft  dunng  LO-LOCAT  Phases  1  and  II,  is  mated  to  this  calibrated 
head  and  is  identical  in  extenor  dimensions  to  the  one  used  in  the  wind  tunnel.  Therefore,  the 
aerodynamic  effects  previously  determined  from  the  wind  Pannel  and  ground  tests  are  applicable 
to  the  probe-shell  combinations  used  on  the  T-33A  airplane.  To  check  for  any  effects  on  a 
and  8  caused  by  the  pressure  field  generated  by  the  airplane  fuselage,  the  wind  tunnel  results 
were  checked  at  the  beginning  of  the  Phase  111  program  by  the  confetti  cloud  calibration 
technique  (Appendix  II).  Results  of  this  testing  verified  the  wind  tunrel  static  calibration  results. 
Figure  3.1  shores  the  agreement  of  the  results  from  these  two  methods  of  testing.  The  results 
indicate  that  a  and  8  measurements  are  not  affected  by  the  T-33A  fuselage  Therefore,  wind 
tunnel  derived  equations  of  a  and  8  are  being  used  in  processing  the  LO-LOCAT  Phase  Ill 
data.  These  equations  are  given  in  Appendix  II. 

The  frequency  response  characteristics  due  to  transduction  effects  were  determined  for 
LO-LOCAT  Phase  III  testing.  With  the  probe  installed  on  the  T-33A  airplane,  the  probe  pressure 
system  was  subjected  to  step  inputs.  The  transduction  effects  were  determined  from  the  resultant 
data  by  utilizing  Fourier  transform  techniques.  This  method  of  testing  and  analysis  is  detailed  in 
Appendix  II. 

The  transduction  effects  thus  determined  were  then  combined  with  the  known  aerodynamic 
effects  obtained  from  wind  tunnel  and  ground  testing.  The  addition  of  these  two  effects 
determined  the  overall  frequency  response  characteristics  of  the  probe.  This  overall  frequency 
response  was  found  to  be  the  same  as  that  for  the  Reference  1.3  testing.  The  gain  factor  is  shown 
plotted  versus  frequency  in  Figure  3.2.  The  reciprocal  of  this  gain  factor  is  used  as  the  frequency 
response  compensation  for  a  and  8  in  processing  the  Phase  111  LO-LOCAT  data.  The  value  of 
compensating  the  probe  pressures  for  frequency  response  characteristics  is  discussed  in  Reference 
1.2. 

In  addition  to  the  testing  and  analysis  done  for  the  probe  angles  of  attack  and  sideslip,  static 
pressure  position  error  was  also  investigated.  The  airspeed  systems  of  both  the  probe  and  airplane 
(pilot's  airspeed  system)  were  calibrated  dunng  flight  by  a  tower  The  techniques  of  this  airspeed 
calibration  method  are  discussed  in  Appendix  II. 

The  pressure  coefficients  for  the  probe  and  for  the  pilot's  airspeed  systems  are  plotted  versus 
indicated  Mach  number,  as  shown  in  Figure  3.3  It  was  determined  from  these  data  that  a  value  of 
0.045  for  the  probe  pressure  coefficient,  4  P,  /q.  ,  was  appropriate  for  the  T-33A  gust  probe. 
Further  substantiation  of  this  pressure  coeffcient  was  established  by  comparison  to  coefficients 
from  other  installations  of  this  same  probe.  The  comparison  is  shown  in  Figure  3.4  for  the  B-52, 
C-131B,  and  T-33A.  This  figure  illustrates  that  for  the  T-33A  installation,  a  value  of  0.04S  is  very 
close  to  the  expected  value.  Impact  pressure,  q  .  ,  and  true  static  pressure,  Ps,  are  then  calculated 
as  shown  in  Appendix  III. 

PROBh  MOTION  SFNSORS 


LO-LOCAT  instrumentation  and  computer  programs  are  designed  such  that  probe  attitude  angles 
may  be  obtained  from  either  attitude  gyro  data  (measured  armies)  or  by  integrating  the  output  of 
the  rate  gyros  (calculated  angles)  Also,  vertical  and  lateral  accelerations  may  be  obtained  either 
from  accelerometers  installed  in  the  probe  or  at  the  airplane  center  of  gravity  (e  g  ).  Since  these 
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parameters  are  interchanged  in  the  event  of  a  malfunction,  an  investigation  was  accomplished  to 
evaluate  their  equality.  This  evaluation  was  made  using  data  obtained  during  flight  in  smooth  air, 
turbulence,  and  simulated  turbulence  The  effects  of  the  turbulence  were  simulated  by 
maneuvering  the  airplane  in  smooth  air. 

Power  spectra  of  measured  and  calculated  angles  were  computed  for  each  of  these  samples.  They 
are  presented  in  Figures  3.5  through  3.12.  As  can  be  seen  from  these  data,  the  agreement 
between  the  two  measurements  is  very  good.  The  largest  difference  in  the  spectra  occurred  for 
angle  of  yaw  for  the  smooth  air  sample.  However,  the  power  level  is  low  and  the  difference 
in  the  standard  deviation  is  small,  as  can  be  seen  from  the  data  in  Table  3. 1 . 

TABLE  3.1 


COMPARISON  OF  MEASURED  AND  CALCULATED 
ATTITUDE  ANGLES 


Sample 

Type 

Variable 

Standard  Deviation  (degrees) 
Calculated  Measured 

Variable  Variable 

Smooth  Air 

Pitch 

0.28 

0.26 

Roll 

0.79 

0.74 

Yew 

0.37 

0.23 

Turbulence 

Pitch 

0.89 

0.79 

Roll 

2.60 

2.50 

L  .  . 

Yaw 

.  -  . 

0.95 

0.90 

Simulated 

Pitch 

0.96 

0.93 

Turbulence 

Roll 

• 

* 

Yaw 

0.9! 

_ 

0.84 

•Roll  rate  exceeded  the  limitations  of  the  recording  system  during  this  test, 
therefore,  roll  angle  comparisons  cannot  be  made 


Since  the  accelerometers  installed  at  the  probe  and  c.g.  arc  in  different  vibration  environments,  the 
effect  of  using  these  accelerometers  in  calculating  gust  velocn.es  was  evaluated  by  comparing  the  gust 
velocities  calculated  using  the  output  from  eacn  set  of  instruments. 

Power  spectra  of  gust  velocity  calculated  using  accelerations  measured  at  both  the  probe  and  e  g. 
were  calculated  for  a  turbulence  condition  and  are  presented  in  Figures  3, 13  through  3.15.  The  time 
series  standard  deviations  obtained  from  these  data  are  tabulated  in  Table  3.2. 
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TABLE  3.2 

COMPARISON  OF  GUST  VELOCITIES  COMPUTED  USING 
PROBE  AND  C.G.  ACCELEROMETERS 


Sample 

Type 

Component 

Component 

Standard  Deviation  (  ) 

Probe 

Acceleration 

C.G. 

Acceleration 

Turbulence 

u 

3.98 

3.99 

v 

3.89 

4.19 

w 

3.29 

3.13 

As  can  be  seen  from  Table  3.2  data,  there  is  little  effect  caused  by  the  interchange  in  the 
accelerometer  measurements  in  either  the  <rt  values  or  the  spectra  shape  It  was  determined 
that  calculated  and  measured  angles,  and  probe  and  c.g.  accelerations  may  be  interchanged  with 
no  harmful  effect  on  the  final  data. 

GUST  PROBE  AND  BOOM  NATURAL  FREQUENCY 

The  gust  probe-boom  assembly  was  designed  to  have  a  natural  frequency  greater  than  10  cps  since 
the  frequency  range  of  interest  for  LOLQCAT  is  0.04  to  10  cps. 

The  probe  and  boom  assemblies  used  during  Phases  1  and  II  of  the  LO-LOCAT  Program  had 
natural  frequencies  of  approximately  1 5.6  cps  (Reference  3.1).  Since  the  boom  used  during  Phases 
I  and  II  was  shortened,  for  the  U9e  on  the  T-33A,  it  was  expected  thar  the  natural  frequency  of 
the  Phase  III  assembly  would  be  greater  than  15.6  cps.  The  rigidity  of  the  boom  attachment  to  the 
airplane  has  a  large  influence  on  vibrational  characteristics;  therefore,  a  check  was  made  to 
determine  the  natural  frequency  of  the  assembly  for  the  Phase  III  Program.  Details  of  this  check 
and  associated  plots  of  the  vibration  characteristics  are  presented  in  Appendix  II. 

The  results  of  this  test  show  that  two  predominant  frequencies  are  present.  These  frequencies  are 
16.7  and  18.5  cps  in  both  the  vertical  and  lateral  dire'-  ..ons  (see  Appendix  II).  There  was  not  any 
aignillcant  vibration  between  the  frequencies  of  0.04  and  10  cps, 

CALCULATED  GUST  VELOCITY  ERRORS 


Longitudinal,  lateral,  and  vertical  gust  velocity  components  are  computed  using  the  equations 
given  in  Appendix  III.  Earth-referenced  probe  motion  is  sensed  in  terms  of  acceleration.  Probe 
velocities  are  computed  by  integrating  the  acceleration  signals.  A  small  bias  error  in  the  integrand 
will  produce  a  large  ramp  function  en or  in  the  resultant  integral  over  a  long  time  period.  This 
error  is  alleviated  by  establishing  a  zero  mean  for  the  integrand  prior  to  integration.  Results  of 
a  test  designed  to  demonstrate  the  effect  of  mean  removal  are  discussed  in  Reference  1 .2. 
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Maneuvers  were  performed  in  relatively  smooth  air  to  evaluate  the  effectiveness  of  motion 
instrumentation  in  the  removal  of  probe  motion  effects  in  computed  gust  velocities.  The 
maneuvering  used  to  simulate  airplane  motion  in  turbulence  was  comparable  to  the  gust  motion 
experienced  during  the  movement  of  an  actual  turbulence  sample  where  gust  velocities  reached  55 
fps,  and  their  average  »t  was  7.4  fps.  Maximum  excursions  of  probe  attitude  angles  and  rates 
of  change  of  attitude  for  the  two  cases  are  shown  in  Table  3.3. 


TABLE  3.3 


COMPARISON  OF  AIRPLANE  MOTION’  DURING 
ACTUAL  AMD  SIMULATED  TURBULENCE 


I 


Variable 

PEAK  TO  PEAK  VALUES 

Turbulence 

Sample 

Simulated 

Turbulence 

Condition 

Pitch  Angle-deg. 

16.7 

8.8 

Roll  Angle-deg. 

i  8 .  A 

16.7 

Yaw  Angle-deg. 

10.  1 

9.4 

Pitch  Rate-deg/sec 

29.0 

28.7 

Roll  Rate-deg/sec 

25.1 

30.8 

Yaw  Rate-deg/sec 

29.0 

16.7 

Plots  of  the  gust  velocity  data  calculated  from  the  simulated  turbulence  sample  are  presented  in 
Figure  3,16.  As  evidenced  by  this  data,  the  incremental  fluctuations  are  generally  ±2  fps  with  4.7 
fps  as  the  maximum.  Time  history  plots  of  gust  velocities  for  the  turbulence  sample  (Test  170, 
Leg  1)  are  presented  in  Section  V,  Figure  5.32. 

An  analysis  of  errors,  other  than  drift,  involved  in  the  gust  velocity  calculations  was  accomplished. 
Negligible  airplane  motion  and  small  angles  were  assumed  to  obtain  simplified  equations  for  the 
probable  error  of  the  three  gust  velocity  components: 


A  U  Or 

.6745 

4V 

A  v  oi 

.  6745 

l  (  V  A  3  ) 

A  v  a 

.6745 

[ (Via) 

(flAV)  +  (V  A* ) 

(orav)2  ♦  (vae)2 


(3.1) 

►  (  *  A  V  )  ]  (3  2) 

,  v*  1  1/2 

(»AV)  ]  (33) 
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error  in  the  measurement  of  airspeed  was  determined  to  be  3.2  fps  from  the  mean  square  of  the 
various  terms  in  the  differential  of  the  true  airspeed  equation  (Appendix  III)  assuming: 


P=  =  2:  in.  Hfc 
in.  Hg 
Ta  =  5^0  0R 
Kc  =  O.996 


A  p  ,  =  +  C  .  15  in.  iig 
Aq  .  =  +  U.C6  in.  Hg 
4Ta  =  t  2  °R 
dK.  =  i  0. 00- 


Gust  velocity  errors  were  then  determined  using  this  error  (  A  V  =  ±  3.2  fps)  and  the  following 

errors  assumed  for  angular  measurements: 

A  t  ar.d  Ad  =  1  0.5  Deg. 

A  fi  ar.d  Ac*  =  t  0.2C1  Deg. 

In  addition,  the  fol'owing  constant  values  were  used: 

a  and  g  =  6  Deg.  i/.  =  5  Dog. 

and 

V  =  620  fpo  q  =  1C  Deg. 

The  greatest  value  of  the  above  calculations  is  in  the  study  of  the  relative  contribution  of  various 
errors  in  the  computation.  Errors  deteimined  by  the  analysis  are  shown  in  Table  3.4. 


TABLE  3.4 


PROBABLE  GUST  VELOCITY  ERRORS 


Gust  Velocity 
Component 

Error  -  fps 

Cont  ribut  tons 

longitudinal 

t  3.0 

i 

A  V  -  iOOT 

Lateral 

i  3.9 

A  0  -  664 

A  $  -  21% 

A  V  -  7’. 

Vertical 

i  3.9 

A  0  -  66.. 

A  a  -  27% 

A  V  -  7'. 
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The  estimates  shown  in  Table  3.4  are  probably  pessimistic  due  to  the  fact  that  it  is  impossible  to 
account  for  the  cancelling  effects  of  errors  associated  with  individual  terms  in  the  gust  velocity 
equation.  In  fact,  the  maximum  gust  velocity  values  calculated  for  several  routine  turbulence 
samples,  where  the  air  happened  to  be  smooth,  have  been  considerably  less.  During  one  such 
sample,  for  example,  longitudinal  gust  velocity  varied  from  -2.7  to  +1.8  fps;  lateral  from  -0.9  to 
+  12  fps;  and  vertical  from  -1 .4  to  +1.2  fps. 

Based  on  these  data,  it  is  concluded  that  the  probable  error  in  LO-LOCAT  Phase  111  data  is  ^2 
fps  with  a  maximum  error  of  -4  fps  for  extreme  turbulence 

CALCULATED  GUST  VETOCITY  DRIFT  ANALYSIS 

An  evaluation  of  the  low  frequency  drift  in  unfiltered  calculated  gust  velocity  was  accomplished. 

The  test  consisted  of  performing  a  normal  instrumentation  prefhght  and  recording  data  for  two 
five  minute  intervals  spaced  two  hours  apart.  During  testing  the  airplane  was  in  a  hangar  and 
isolated  so  that  no  actual  physical  inputs  to  the  instrumentation  would  occur.  The  data  were  then 
processed  through  the  gust  velocity  equations.  Both  measured  (attitude  gyro  output)  and 
calculated  (integrated  rate  gyro  output)  probe  attitudes  were  used  in  the  calculations  to  evaluate 
the  drift  caused  by  each  set  of  measurements.  All  parameters  that  form  the  gust  velocity  equations 
were  used  as  recorded  except  for  q  jfairspeed),  F,  ^static  pressure),  and  T»  (temperature).  These 
parameters  were  set  equal  to  constants  that  would  yield  a  VT  (true  airspeed)  equal  to  a  value 
representative  of  the  inflight  airspeed. 

The  unfiltered  data  obtained  from  these  calculations  are  presented  in  Figures  3. 17  through  3.19 
The  letter  designations  on  the  figures  indicate  whether  calculated  (C)  or  measured  (M)  attitude 
angles  were  used.  The  subscripts  one  and  two  refer  to  the  first  and  second  recording  intervals.  The 
data  obtained  show  that  the  drift  increased  with  an  increase  in  instrumentation  operating  time 
when  using  calculated  angles.  The  drift  in  the  data  using  measured  angles  was  generally  more  linear 
than  the  data  obtained  using  calculated  angles. 

To  further  evaluate  the  effects  of  drift  and  determine  where  tts  greatest  influence  occurs,  the  gust 
velocity  equation  was  divided  into  three  parts;  probe  motion  terms,  indicated  gust  velocity  terms; 
and  constant  terms.  The  probe  motion  terms  define  the  changes  in  the  longitudinal,  lateral,  and 
vertical  probe  velocities  in  relation  to  the  earth  axes.  Their  computations  involve  the  integiation 
of  acceleration  information  obtained  from  probe  acceleration  measurements  and  rotation  to  the 
earth  axes  using  probe  attitude  measurements.  The  constant  terms  define  the  initial  conditions 
and  remove  the  earth  gravitational  effects  from  probe  vertical  acceleration  measurement.  They 
cannot,  therefore,  contnbute  to  the  drift.  The  indicated  gust  velocity  terms  define  the 
longitudinal,  lateral,  and  vertical  components  as  a  combination  of  probe  velocity  and  gust  velocity 
in  relation  to  the  earth  axes.  Their  computation  involves  attack,  sideslip,  and  probe  attitude  angles 
and  true  airspeed. 

The  drift  in  each  gust  velocity  component  was  found  to  be  almost  totally  caused  by  the  probe 
motion  terms  They  are  not  shown  because  they  did  not  vary  significantly  from  the  true  gust 
velocity  data  presented  in  Figures  3  17,  3.18,  and  3.19.  The  drift  in  the  indicated  gust  velocity 
terms.  Figures  3.10  and  3.2 1  did  not  exceed  10  fps  and' was  somewhat  cyclic  in  nature.  The  drift 
in  indicated  longitudinal  gust  velocity  is  not  presented  because  the  maximum  change  was  less  than 
one  fps. 
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The  high  drift  experienced  in  the  probe  motion  terms  was  due  to  very  low  amplitude  drifting  in 
the  instrumentation.  This  drift  varied  in  amplitude  in  a  periodic  manner.  Plots  of  integrated 
accelerometer  signals  and  calculated  and  measured  angles  are  presented  in  Figures  3.22  through 
3.26 

As  the  drift  experienced  in  the  individual  parameters  was  less  than  0.2  fps  for  the  integrated 
accelerometers  and  less  than  0.7  degree  for  the  calculated  angles,  it  was  not  readily  apparent  how 
thia  small  drift  could  build  up  to  high  values.  The  probe  motion  terms  were  further  segmented  as 
follows  and  each  segment  analyzed  separately: 


^  3.4  'i 


where: 


»u  •  Yi  ♦  7?  +  >3 

■Yj^  ■  -nK  co«  *  cos  | 

>2  •  n,  (  ■in  *  coa  *  -  coa  *  ain  $  sin  <t> ) 
>3  •  n,  (  aln  ♦  ain  *  ♦  coa  *  ain  g  cob  0 ) 


\ 

I 

i 

!  where: 


! 

1 


u. 


By  dt 


*v  ■  M  j.  *  1*2  *  **3 

■  n*  ain  ♦  coa  • 

**2  •*  nr  (  co*  ♦  001  4  *  ,ln  ♦  Bia  •  ,ln  *) 

■  n,(  coa  if  ale  $  -  ain  if  ain  f  coa  0 ) 


(3.5) 
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&  v  dt 


(3.6) 


where,  aw“  ti+t24,t3 

r,_  ■  n x  sin  § 
r  2  ■  -ny  sin  4  cos  0 
r3  -  n,  cos  4  cos  i 

The  data  selected  for  presentation  were  obtaine  1  using  calculated  angles  and  exhibited  the 
maximum  drift  amplitude  The  data  which  may  be  f  iund  in  Table  3.5  show  that  the  effect  of 
y  _  i  y  2  >  >  M2,  r  .  »nd  t  2  is  insignificant  and  that  7, ,  M  3  and  r3  contribute  almost  all 
of  the  drift. 


table  3.5 

COMPARISON  OF  DRIFT  IN  GUST  VELOCITY  EQUATION  TERMS 


Recording 

Interval 

No. 

— 

Gust 

Velocity 

Component 

Equation 

Variable 

Variable 

Mean 

(x) 

t 

g/  x  dt 

O 

2 

u 

?! 

-6.8S  x  10 

2.88  x  10  ~ 
6.58  x  10  ' 

-6.64  x  10  ^ 
2.20  x  10  j 
6.35  x  10 

v 

S} 

"3 

-2.21  x  10  “® 
5.20  x  10  ", 
2.38  x  10  "2 

-2.13  x  10  ~U 
5.02  x  10  , 
2.29  x  10  1 

„ 

T  1 

T  2 
r  3 

-1.75  x  10  ~l 
-7.43  x  10  ", 
-3.00  x  10  "3 

-1.69  x  10  "3 
-7.17  x  10  "3 
-2.90  x  10  1 

It  is  believed  that,  since  vertical  acceleration  is  common  to  all  three  equations,  it  is  the  primary 
source  of  drift.  This  is  caused  in  part  by  the  reduced  resolution  of  the  vertical  acceleration 
measurement  when  compared  to  lateral  or  longitudinal  acceleration.  The  resolution  of  the  reduced 
data  is  approximately  .005  g’s/digital  count  for  vertical  acceleration,  .002  g’s/digital  count  for 
lateral  acceleration,  and  .001  g’s/digital  count  for  longitudinal  acceleration.  It  should  be  noted 
that,  in  Figures  3.22  and  3.23,  integrated  vertical  acceleration  exhibited  the  maximum  amount  of 
drift.  It  should  also  be  noted  that  roll  angle  is  common  to  all  three  equations  (y,  u*„and  r?) ; 
m.  ,vever,  the  resolution  of  roll  angle  is  the  same  as  for  pitch  and  yaw. 

A  comparison  of  the  drift  in  u,  v,  and  w  revealed  that  u  and  v  generally  drifted  to  higher 
amplitudes  than  w.  The  w  drift  never  exceeded  30  fps  whereas  v  drifted  to  values  greater  than  1 25 
fps  and  u  drifted  to  values  in  excess  of  500  fps.  These  characteristics  have  been  consistent  in  all 
gust  velocity  calculations. 

The  drift  frequencies  observed  in  the  unfiltered  data  were  never  greater  than  .04  cps  and  were 
generally  much  lower.  Drift  frequencies  were  generally  .02  cps  or  less. 

Example  plots  illustrating  the  removal  of  drift  from  the  ground  test  data  are  presented  in  Figure 
3.27.  This  figure  contains  the  C  ?  data  of  Figures  3.17,  3.18,  and  3.19,  after  filtering  with  the 
LO-LOCAT  digital  filter.  Figures  3.28  and  3.29  contain  time  history  plots  of  ur.filtered  and 
filtered  data  obtained  in  smooth  air.  Figures  3.30  and  3.3 1  are  unfiltered  and  filtered  data  from  a 
turbulence  sample.  As  can  be  seen  from  these  data,  the  removal  of  the  undesirable  drift  in 
unfutered  gust  velocity  is  successfully  achieved  through  the  use  of  the  LO-LOCAT  high-pass 
digital  filter  Details  of  the  digital  filter  may  be  found  in  Apperdix  III. 

In  summary,  drift,  to  some  degree,  is  inevitable  in  the  present  state-of-the-art  processes  required  to 
compute  gust  velocity.  Low  frequency  (f  <  0.04  cps),  high  amplitude  drift  in  calculated  gust 
velocity  is  generated  by  small  amounts  of  drift  ir.  the  individual  recorded  parameters,  coupled 
with  the  mathematical  operations  required  to  compute  gust  velocities,  especially  integration  The 
primary  source  of  drift  is  attributed  to  the  vertical  acceleration  measurement.  The  variation  of  this 
drifting  with  time  is  unique  for  each  individual  gust  velocity  calculated,  Drift  in  gust  velocities 
computed  using  attitude  gyro  data  appears  to  be  more  linear  than  the  drift  which  occurs  when 
computations  are  based  on  attitude  angles  calculated  from  rate  gyro  information.  Successful  drift 
removal  may  be  attained  by  accurate  application  of  a  digital  filter. 
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NOISE  EVALUATION 


A  study  of  instrumentation  system  noise  levels  was  made  to  evaluate  the  effects  on  calculated  gust 
velocities  For  the  purposes  of  this  study,  noise  was  defined  as  all  spurious  oscillations,  at  high  and 
low  frequencies,  which  appear  in  the  data.  During  LO-LOCAT  Phases  I  and  II,  it  was  determined 
that  some  type  of  processing  technique  must  be  employed  to  remove  noise  from  the  static 
pressure,  OAT,  and  attitude  gyro  measurements  (Reference  1.2). 

A  100-point,  moving  time  averaging  technique  was  developed  to  reduce  the  noise  on  static 
pressure  and  temperature.  Evaluation  of  this  technique  was  accomplished  using  data  obtained 
during  flight  in  smooth  avr.  It  was  determined  that  by  smoothing  these  two  parameters  no 
distortion  occurred  in  the  longitudinal  spectra  at  low  frequencies  and  there  was  a  reduction  of 
noise  levels  in  the  high  frequency  end  of  the  spectra.  No  smoothing  effects  were  observed  in 
vertical  arH  lateral  gust  spectra.  As  a  result  of  this  test,  OAT  and  static  pressure  data  are  also 
being  smoothed  during  Phase  111  using  the  1  OO-point,  time-averaging  routine. 

During  Phases  1  and  II,  it  was  also  determined  that  high  frequency  noist  must  be  eliminated  from 
the  measured  attitude  angles.  This  was  accomplished  by  calculating  d<7  ,  i.e.,[fft  “cr]  for 
the  gust  velocity  data.  An  evaluation  of  dtr  givos  an  indication  of  the  amount  of  turbulence 
that  exists  between  frequencies  of  10  and  20  cps,  assuming  no  noise  This  occurs  because  o  is 
calculated  from  the  truncated  spectrum  from  0.04  to  10  cps  and  a  .  is  calculated  from  the 
time  series  Time  series  data,  considering  the  high-pass  digital  filter  and  the  low-pass  ground 
station  analog  filter,  contain  frequencies  ftomO.04  to  20  cps.  In  addition,  it  was  further  proved 
during  Phases  I  and  11  that  the  measured  attitude  angles  were  contributing  significant  noise  to  the 
gust  data  at  frequencies  less  than  lOcps.  This  was  proven  by  comparing  the  c  valur  for  gust 
velocity  spectra  calculated  both  with  measured  and  with  calculated  angles.  It  was  found  thai 
significant  noise  existed  in  the  attitude  measurement  at  frequencies  greater  than  five  cps.  A 
low-pass  filter  was  therefore  developed  to  remove  noise  above  four  cps.  The  filter  is  an  Ormsby 
low-pass  (See  Appendix  III j  with  a  cutoff  frequency  of  4  cps.  During  Phase  111,  the  filter  cutoff 
frequency  was  changed  to  three  cps  to  permit  removal  of  more  noise.  The  effect  of  this  filter  can 
be  seen  in  Figures  3.32  and  3.33  which  show  the  noise  level  on  measured  angles  both  before  and 
after  filtering.  Since  the  frequency  of  airplane  motion  is  one  cps  and  lower,  only  noise  is  removed 
by  this  filter  and  the  frequencies  in  the  airplane  motion  range  were  not  affected. 

The  complete  instrumentation  system  was  evaluated  using  ooth  ground  ar.d  flight  data  to 
determine  the  level  of  turbulence  below  which  LO-LOCAT  Phase  III  data  could  no  longer  be 
considered  valid  due  to  low  signal-to-noise  ratio 

The  ground  test  data  consisted  of  the  same  data  used  during  the  drift  analysis.  The  drift  data  in 
Figures  3.17  through  3.19  were  processed  through  the  high-pass  filter,  described  in  Appendix  Ill, 
and  power  spectra  and  standard  deviations  were  calculated  The  minimum  valid  (Tt  was  calculated 
assuming  that  10  decibels  (db)  is  the  minimum  signal-to-noise  power  ratio  required  for  valid  data. 
This  minimum  value, crt> , was  calculated  by  increasing  a.  by  a  factor  of  the  square  root  of  10  or: 


9 


(  10 


1/2 


(3.8) 


29 


Since  »t2  and  the  power  spectral  density  function  are  related  by 


a 


2 

t 


*  v k)  dk 


(3.9) 


increasing  9  2  by  ?.  factor  of  10  is  equivalent  to  increasing  the  power  level  of  tie  spectra  by 
10  db  end  integrating  the  new  function  to  obtain  the  minimum  a  2. 

The  result  of  this  test  is  presented  in  Table  3.6. 


TABLE  3.6 


MINIMUM  VALID  STANDARD  DEVIATION  VALUES 


Recording 

Interval 

Attitude 

Angles 

Used 

— 1 

Variable 

*t 

1 

Calculated 

u 

0.01 

0.03 

v 

0.02 

0.38 

w 

0,06 

0 .  i  9 

Measured 

u 

0.03 

0 .  u- 

V 

0.52 

1 . 64 

w 

0.43 

1.52 

2 

Ca lculated 

u 

0.02 

0.06 

V 

0.11 

0.35 

w 

0.07 

0.22 

Measured  ¥ 

u 

0.03 

0.09 

V 

0.10 

0.32 

V 

0.49 

1.55 

•Measured  angies  of  pitch  and  roll,  calculated  angle  of  yaw. 


The  lower  <7t  value  for  tin  lateral  measurement  during  the  second  interval  using  measured 
angles  is  probably  the  result  of  u  i  rg  calculated  instead  of  measured  yaw  angle.  Calculated  angles 
apparently  contribute  less  enor  than  measured  angles  to  the  standard  deviation  values.  This  is  a 
result  of  the  smoothing  effect  of  integration  on  high  frequenc.es  [f>l/(2rr)]when  using  rates,  and 
the  additional  high  frequency  noise  associated  with  the  attitude  gyros 

It  would  appear  from  the  above  table  that  there  are  some  cases  when  a  t  should  have  a  value  of 
at  least  1 .64  fps  to  be  assured  of  valid  turbulence  data  However,  it  must  be  realized  that  the  noise 
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vanes  slightly  from  time  to  time.  This  is  readily  seen  m  the  fact  that  at  values  as  low  as  0.34, 
0.28,  and  0.28  for  u,  v,  ar.d  w,  respectively,  have  been  observed  in  the  flight  data.  This  would 
yield  dtc  values  on  the  order  of  1 .07  fps  for  valid  data. 

The  shape  of  the  power  spectral  density  function  of  any  variable  will  be  degraded  if  the 
signal-to-noise  ratio  becomes  too  low.  It  will  also  be  degraded  if  the  time  (or  space)  function  is 
clearly  nonhomogeneous.  Ther'ore,  the  turbulence  data  which  had  irregularly  shaped  spectra 
Ci  e..  did  not  have  a  well  defined  inertial  subrange  with  a  -5/3  slope)  were  compared  with  the 
turbulence  data  having  a  well  defined  inertial  subrange  (reguiariy  shaped  spectra).  Figure  3  34 
contains  frequency  histograms  of  the  stai.oard  deviations  of  both  sets  of  data.  In  developing  the 
histograms,  only  turbulence  data  having  ox<  2.0  fps  were  used.  A  plot  of  the  ratio  of  the 
frequency  of  occurrence  of  irregularly  shaped  to  regularly  shaped  spectra  for  each  band  of<7tis 
presented  in  Figure  335  As  indicated  by  this  curve,  there  are  twice  as  many  irregularly  shaped 
spectra  at  values  of  approximately  1.2  fps.  At  ut  =  1 .4  fps  there  are  an  equal  number  of  both 
sets. 

Considering  the  ground  test  data,  inflight  standard  deviations,  and  power  spectra  data,  it  is 
believed  that  a  minimum  value  <rta=  1 .4  fps  will  yield  good  results  in  the  final  data  with  the  small 
number  of  noisy  samples  contained  in  the  final  data  having  little  effect  on  the  results 
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DIFFERENCE  BETWEEN  WIND  TUNNEL  CALIBRATION 
GAIN  FACTOR  -  db  AND  CONFETTI  CALIBRATION  -  Degrees 


O  Angle  of  Attack 
□  Angle  of  Sideslip 


2 


-2h 

I _ I - 1 _ I _ I _ I _ 1 

-6  -4  -2  0  2  4  6 

CONFETTI  CAUBFATIOH  ANGLES  OF  ATTACK  AIT)  SIDESLIP  =  Degrees 

iigure  3*1  Comparison  of  Wind  Tunnel  ar.d  "Confetti" 
Calibration  Results 


Figure  3.2  Dynamic  Response  of  Angles  of  Attack  and 
Sideslip  Measurements 
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GAIN  FACTOR  -  Output/input 


COEFFICIENT 


n1  ■  1  1  1  1  1  L  1  1 - l  i - i - 1 - 1 - j 

0  .2  >  .6  .8  1.0  1.2  1.4 

BOOM  SIATION/EFFECTI\E  FUSELAGE  DIAMETER 


F'gure  3-1*  Gust  Probe  Pressure  Coefficients  for 
'  JB-52E,  C-131B,  and  T -3 3A  Airplanes 
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Figure  3.11*  Lateral  Gust  Velocity  PSD  Calculated  Using 
Inputs  froci  Different  Acceleroneters 


ruguiKrr  •  ept 


Figur*  3*15  ^«rtlo»l  Ou»t  V«locity  PSD  Calculated  U*lng 
Input*  from  Different  Accelenometers 


Btorles  of  Simulated  Turbulence 
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Xj  -  Longitudinal  Velocity, 
Y^  -  Lateral  Velocity, 

Zl  -  Vertical  Velocity, 


X2  -  Longitudinal  Velocity, 
Yg  -  Lateral  Velocity, 

12  ~  Vortical  Velocity, 


TIME  -  Second* 


Pigu  re  3*23  Drift  in  Integrated  Acceleration  -Interval 
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Figure  3*30  Unfiltered  Guet  Voloclty  -  Turdul encu 
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Figure  3,31  Filtered  Gupt  Velocity  -  Turbulence 
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SECTION  IV 


METEOROLOGICAL  DATA  ANALYSIS 

Meteorological  parameters  are  being  recorded  to  determine  their  relationship  to  turbulence.  This 
effort  includes  development  of  forecasting  models,  correlation  of  observed  turbulence  and 
meteorological  data,  and  statistical  analysis  of  certain  basic  parameters.  Global  Weather  Center 
(GW O  forecasts  are  compared  to  Phase  Ill  gust  velocity  rms  values.  This  report  contains  Phase  111 
data  collected  at  the  McConnell  and  Edwards  .Air  Force  Base  locations  and  selected  Phases  1  and 
11  data. 

Three  types  of  meteorological  data  are  being  recorder  that  recorded  by  the  airplane  in  flight, 
that  recorded  by  rawinsonde,  and  that  derived  from  surface  observations.  The  meteorological  data 
recorded  by  the  airplane  are  primarily  obtained  during  the  surveys  performed  prior  to  the 
recording  of  each  4,5-minute  turbulence  sample.  During  these  surveys,  air  temp-mature,  ground 
surface  temperature,  wind  speed  and  direction,  absolute  and  pressure  altitude,  airplane  speed  and 
heading,  and  pilot  comments  concerning  the  state  of  the  atmosphere  are  recorded.  Utlng  the 
survey  data,  vertical  gradients  are  determined  to  define  lspse  rate,  wind  ahear,  and  atmospheric 
stability  parameten.  Standard  Weather  Bureau  and  Air  Force  rawinsonde  observations  are  u  tilled 
whenever  ipplicable;  however,  the  greatest  reliance  ia  placed  on  the  special  observations  provided 
at  each  base  by  the  6th  Weather  Squadron  (Mobile).  A  Boeing  meteorologist  at  each  baae  of 
opentioni  provides  weather  and  turbulence  forecuts  and  records  pertinent  data  for  subsequent 
analyse*.  Flight  weather  condition*  are  forecut  in  cooperation  with  the  Air  Force  base  weather 
personnel.  Hourly  surface  observations  are  received  by  teletype  from  adjacent  first  order  weather 
stations  and  special  surface  observation*  are  provided  by  cooperative  observers. 

A  summary  of  the  climatological  record  during  each  flight  period  Is  given  to  compare  the  average 
observed  condition!  during  the  period  to  the  put  record  (or  normal).  Racords  from  the  nearejt 
I’.S.  Waathar  Bureau  station  ware  used  for  this  comparison.  Thase  stations  are  located  at  the 
Wichita,  Kansu  Municipal  Airport  for  the  McConnell  rout*  and  at  the  Kem  County  Air  Terminal, 
Bakersfield,  California,  for  the  Edwards  routs.  The  Bakareflald  data  war*  not  available  for  this 
report.  MoConntU  dat?  presented  in  Fig.  4.1  show  the  departure  of  observed  date  from  the 
normal. 

Advert*  waathar  was  responsible  for  cancellttion  of  14%  of  thy  scheduled  turbulence  sampling  et 
th*  McConnell  route  and  30%  at  the  Edwards  route. 

Selected  meteorological  data  recorded  by  th*  airplane  instrumentation  art  Listed  in  the  Test  Log, 
Appendix  VI.  After  completion  of  the  program,  th*  rawinsonde  observations  and  the  specially 
prepared  large  Kale  synoptic  cherts  fer  each  Phase  III  flight  will  be  availsble  from  the  USAF 
Environmental  Technical  Application*  Center  (Data  Processing  Division)  in  Asheville,  North 
Carolina. 

WIND  VELOCITY 

Th*  relationihip  between  wind  velocity  and  gust  velocity  rms  values  Is  being  Investigated  (See 
Section  V),  In  conjunction  with  this  investigation,  the  variation  of  wind  speed  for  different 
geophysical  situations  is  being  studied  using  a  statistical  ipproach.  This  method  could  be  a 
valuable  aid  to  forecuting  wind  speeds  for  low  level  flights  Into  geographic  areas  where 
obtervationi  are  not  available. 
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The  north  and  east  components  of  airplane  measured  wind  velocities  are  computed  at  the  rate  ov 
100  times  per  second,  or  27,000  times  during  a  sample.  These  data  are  then  averaged  as  follows: 
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From  these  average  wind  components,  the  average  wind  during  the  sample  is  then  calculated.  The 
equations  for  these  calculations  are  shown  in  Appendix  111. 

The  data  are  grouped  into  bands  of  wind  speed,  each  band  having  a  width  of  5  fps.  The 
cumulative  probability  is  plotted  versus  the  lower  band  limit.  The  wind  analysis  includes  a  total  of 
741  samples.  Of  this  total,  232  sample*  were  recorded  at  McConnell  and  509  at  the  Edwards 
location.  For  each  category  presented,  a  minimum  of  30  wind  sampler  was  used  for  establithing  a 
cumulative  probability  curve  with  adequate  reliability.  Because  of  the  need  for  a  representative 
number  of  samples  in  each  category,  various  combinations  of  geophysical  situations  could  not  be 
analyzed.  However,  the  data  presented  here  indicate  the  wind  speed  -ends  u  they  vary  with  the 
most  pertinent  individual  category  components, 


During  LO-LOCAT  Phase*  land  I!  (Reference  1,2),  the  category  component*  which  »howed  the 
greatest  effect*  on  wind  cumulative  probability*  were  location  and  »e»*on,  Effects  of  altitude 
above  the  terrain  and  tlme-of-day  were  al*o  found  to  be  present,  but  were  dependent  upon  the 
type  of  terrain  over  which  the  data  were  obtained  and  the  geographic  location. 


For  the  Phase  III  data,  an  analyii*  by  both  lesion  and  location  wa»  not  made  because  the  twe 
category  component*  are  redundant.  At  McConnell  Air  Force  Base,  Ksnias,  approximately  90 
percent  of  the  dsta  were  obtained  during  the  summer  as.J  only  10  percent  during  the  fall.  At 
Edwins  Air  Force  Bsse,  California,  data  obtained  dunng  the  fall  accounted  for  approximately  80 
percent  of  the  total.  Approximately  20  percent  of  the  Edward*  data  were  obtained  daring  the 
winter.  The  cumulative  probability  of  wind  »peed  wu  not  analyzed,  therefore,  at  a  functiomof 
*ea»on.  The  wind  cumulative  probability  i*  preaented  a*  *  function  of  location  in  Figure  4.2.  A* 
can  be  teen,  the  probability  of  encountering  wind  of  a  given  ipeed  I*  greater  for  the  McConnell 
route  than  for  Edward*,  This  same  result  wu  al*o  found  during  Pha*e  1  and  II  of  the  program.  In 
fact,  during  Phase  I  and  II  the  highest  winds  were  recorded  at  McConnell  and  the  lowest  at 
Edwards. 
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Although  the  data  during  Phases  1  and  II  were  obtained  during  all  four  seasons  at  each  location, 
the  McConnell  data  agree  well  with  the  data  obtained  at  that  location  during  Phase  III  (Figure 
4.3).  The  Phase  III  winds  recorded  at  McConnell  during  the  summer  are  of  slightly  greater 
magnitude  than  those  recorded  during  Phases  I  and  II.  The  wind  speeds  recorded  at  Edwards 
during  Phase  III  are  somewhat  less  than  those  recorded  at  that  location  during  Phases  1  and  II 
(Figure  4.4). 

The  probability  of  encountering  a  given  wind  speed  was  higher  at  750  feet  absolute  altitude  than 
at  250  feet  for  data  recorded  over  the  plains  (Figure  4.5).  For  data  recorded  over  the  high 
mountain,  low  mountain,  and  desert  terrain,  there  is  very  little  difference  in  wind  speed 
cumulative  probabilities  at  750  and  250  feet.  These  results  are  shown  in  Figure  4.6,  4.7,  and  4.8. 
For  the  data  included  in  this  analysis,  all  of  the  plains  legs  are  located  at  the  McConnell  route  and 
all  high  mountain,  low  mountain,  and  desert  legs  are  located  Rt  the  Edwards  route. 

Tlme-of-dey  effects  on  wind  cumulative  probability  are  shown  in  Figure  4.9  through  4.13.  For  the 
data  recorded  over  the  McConnell  route  (Figure  4.9),  the  cumulative  probability  for  a  given  wind 
speed  is  greater  at  dawn  than  either  mid-morning  or  mid-aftemoon.  Winds  recorded  over  the 
Edwards  route  showed  a  different  tendency.  As  during  Phases  1  and  II,  these  wind  speeds 
increased  slightly  as  the  time-of-day  progressed  from  dawn  to  mid-moming  to  mid-aftemoon 
(Figure  4.10).  As  shown  in  Figures  4.11,  4.12,  and  4  13,  time-of-day  did  not  have  as  great  an 
effect  on  wind  over  any  of  the  terrains  of  the  Edwards  route  as  was  observed  for  the  plains  data  at 
McConnell. 

The  variation  of  wind  cumulative  probability  with  atmospheric  stability  is  different  for  each 
geographical  location.  At  Edwards  (Figure  4.14),  the  general  trend  is  for  the  wind  cumulative 
probability  to  Increase  for  a  given  wind  speed  at  the  atmosphere  become*  more  unstable.  This 
trend  doe*  not  appear  in  the  McConnell  data  (Figure  4.15).  For  these  data,  there  upDears  to  be 
very  little  change  In  the  expected  wind  speed  as  the  stability  change*. 

The  greatest  wind  velocitiei  obtained  during  this  part  of  the  progrem  were  ove'  *he  plains  Figure 
4.16  pnienta  an  overall  view  of  the  relative  wind  speed  cumulative  probebilKlet  at  a  furctlon  of 
the  terrain  over  which  the  data  were  obtained. 

Wind  speed  and  direction  are  listed  In  Appendix  VI  for  esch  gust  velocity  sample  The 
corresponding  geophysical  category  for  each  sample  is  tlto  listed. 

RICHARDSON  NUMBER 

The  following  equation  It  used  to  calculate  the  values  of  Richtrdson  number  presented  in  this 
report: 
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Probabilities  of  encountering  given  gust  velocity  rms  values  were  determined  for  each  of  20  bands 
(0.6  wide)  of  Richardson  number  from  -2.0  to  +10.0.  Approximately  11  percent  of  the 
Richardson  number  values  were  greater  than  +10,  and  2  percent  were  less  than  -2.  They  were  not 
included  for  this  analysis.  The  data  analyzed  are  presented  in  Figures  4.17  through  4.19.  The  gust 
velocity  rms  values  having  80  percent  probability  of  occurrence  were  plotted  versus  mid-band 
Richardson  number  (Figure  4.20).  This  figure  compares  Phase  I  and  II  to  the  Phase  III  data  and 
indicates  that  tne  probability  of  occurence  of  a  given  rms  level  is  greater  for  the  Phase  HI  data  at 
equivalent  Richardson  number  values  These  data  are  consistent  with  the  higher  rms  gust 
velocities  being  recorded  during  Phase  III  because  of  the  increased  airspeed  of  the  T-33A  airplane 
compared  to  the  C-13  IB  airplanes. 


The  data  presented  here  for  Richardson  number  statistics  must  be  considered  preliminary  because 
of  the  small  number  of  turbulence  samples,  and  their  distribution  in  the  Richardson  number 
bands.  The  number  of  data  points  in  each  Richardson  number  band  should  be  at  least  30  for  a 
good  statistical  sample.  For  this  interim  report,  onlj  the  three  bands  from  -0.2  toil  .6  contained 
over  30  data  points. 

SAMPLING  SITES 

Hourly  surface  observations  and  rawinsonde  soundings  provide  data  used  io  support  the 
LO-LOCAT  Phase  111  flights.  Special  rawinsonde  soundings  are  being  provided  by  the  USAF  Sixth 
Weather  Squadron  (Mobile).  These  special  soundings,  which  are  being  provided  for  each  flight, 
include  soundings  from  the  surface  to  500  mb  (400  mb  at  McConnell).  Additionally,  they  are 
providing  the  pressure,  temperature,  dew  point,  and  wind  velocity  at  levels  of  250,  750,  and  1 ,000 
feet  above  the  ground.  The  balloon  it  underinflated  so  cs  to  rise  at  the  rate  of  500  feet  per  minute, 
Azimuth  and  elevation  data  are  recorded  every  30  seconds.  Whenever  available,  standard 
rawinsonde  observations  from  nearby  fixed  installations  are  also  being  recorded. 

Cooperative  observers  in  the  vicinity  of  each  route  are  asked  to  provide  additional  observations. 
The  observers  are  provided  with  sling  psychrometerv  thermometers,  har.d-held  anemometers, 
obaerving  guideline*,  and  special  recording  forms.  They  provide  wet  and  dry-bulb  temperatures, 
surface  winds,  and  visual  observations  of  cloud  cover,  dust  devils,  and  precipitation. 

Meteorological  informstion  derived  from  debriefing  of  the  pilot  is  also  recorded. 

McConnell  Air  Force  Base,  Kansas  •  IS  August  1968  to  25  September  1968 


For  the  McConnell  route,  the  mobile  rawinsonde  was  located  at  Cedar  Vale,  Kansas,  which  la 
approximately  equidistant  from  ail  the  flight  legs  Since  ill  the  flight  legs  were  over  plains,  one 
sounding  was  considered  representative  of  a  given  flight.  Hourly  surface  observations  were 
obtained  from  Wichita,  Chanute,  and  McConnell  Air  Force  Baae,  Kansas,  and  Ponca  City, 
Oklahoma.  Weather  radar  reports  were  available  continuously  from  McConnell  Air  Force  Base  and 
hourly  from  the  Wichita  Weather  Bureau  Periodic  observations  were  received  from  Bartlesville, 
Oklahoma  Six  cooperative  weather  observer!  volunteered  to  participate  in  thi  program,  located  as 
shown  In  Figure  2.7.  Surface  wind  data  were  available  on  call  from  the  TITAN  II  misaile  sites  near 
legs  1.  2,  7,  and  8. 
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Edwards  Air  Force  Base,  Caiifomia  -  8  October  1968  to  8  January  1969 


The  mobile  rawinsonde  site  was  first  located  about  six  miles  north  of  Oxnard  Air  Force  Base, 
California,  on  8  October  1968  and  was  later  removed  to  Oxnard  Air  Force  Base  on  22  October 
1968.  Edwards  Air  Force  Base  also  provided  special  rawinsonde  soundings  for  each  flight  (prior  to 
23  October  1968,  only  wiresonde  observations  were  provided).  In  addition  to  these  special 
observations,  the  routinely  provided  rawinsonde  data  from  Vandenberg  Air  Force  Base,  Caiifomia, 
for  0400,  1000,  and  1600  PST  were  analyzed  to  assist  in  characterizing  the  upper  air  situation. 
The  terrain  in  the  overall  test  area  is  so  complex  that  it  would  have  taken  several  mobile 
rawinsonde  sets  to  obtain  representative  soundings  for  all  the  flight  legs.  However,  this  was  not 
feasible  for  budgetc*y  reasons.  Hourly  surface  observations  were  obtained  from  the  first  order 
weather  station  shown  on  Figure  2.8.  Only  one  cooperative  observer  (located  at  Juncal  Dam) 
participated  in  the  program. 

forecasting  techniques 


Subjective  Forecasting  Techniques 


! 
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The  Second  Weather  Squadron  at  Offutt  Air  Force  Base,  Nebraska.  Global  Weather  Center, 
provided  special  low-level  turbulence  forecasts  for  the  Phase  III  LO-LOCAT  flights.  These 
forecasts  were  based  on  forecast  values  of  the  vertical  wind  shears,  winu  speeds  from  the  surface 
to  1,000  feet,  pressure  changes,  and  a  value  for  terrain  type.  Forecasts  were  prepared  for  each 
type  of  terrain  and  specifically  for  the  dawn,  mid-moming,  and  mid-afternoon  flights.  The 
forecast  consisted  of  a  nondimensional  number  corresponding  to  a  specific  degree  of  turbulence  as 
follows: 


GWC  Numerical  Forecast 
Less  than  80 
80  to  119 

Graater  than  or  equal  to  120 


Expected  Turbulence 
No  Turbulence 
Light  Turbulence 
Moderate  Turbulence 


An  evaluation  of  the  GWC  turbulence  forecast  wu  performed  by  computing  the  correlation  of 
GWC  numbers  to  observed  gust  velocity  rm»  values.  The  correlation  wsi  computed  for  the  general 
case  (all  values)  by  using  regression  snalysis  with  inputs  of  either  linear  (L)  or  squared  (S)  terms. 
The  resulting  correlation  coefficients  are  listed  in  Table  4.1,  a  value  of  zero  Indicating  no 
correlation  and  a  value  of  1 .0  perfect  correlation.  The  mean  error  between  the  actual  gust  velocity 
rm»  and  that  predicted  from  the  regreialon  equation  is  listed  in  Table  4  1.  Also  listed  art  the 
errors  for  regression  equations  of  rms  values  greater  than  a  selected  threshold  value,  which  is  of 
interest  in  defining  how  well  the  forecast  predicts  only  the  higher  rms  values.  The  best  forecasts 
•re  shown  to  have  been  made  for  the  plains  terrtin  with  errors  of  as  much  as  42  percent  in  some 
other  categories, 

Another  forecasting  tool,  the  Showalter  stability  index,  was  compiled  for  each  flight  and  the  tame 
type  of  regression  snalysis  performed  as  was  done  for  the  GWC  forecasts  This  index  is  an 
Indication  of  either  a  atable  or  unitable  atmospheric  condition  between  the  850  to  500  mb  level. 
Tne  reiulti  are  shown  in  Table  4.2  and  indicate  that  the  Showalter  stability  index  is  poorly 
correlated  with  the  gus'  velocity  rms  values. 


58 


For  comparison,  Table  4.3  lists  the  results  of  regression  analysis  equations  to  evaluate  Richardson 
number  as  a  forecast  tool.  The  errors  are  quite  large  and  the  correlation  coefficients  are  near  zero. 


The  evaluations  represented  by  Table  4.1  through  4.3  emphasize  the  difficulty  being  experienced 
in  the  forecasts  of  low  level  turbulence.  These  results  indicate  the  importance  of  the  need  for 
improved  methods  and  procedures  in  the  performance  of  low  level  turbulence  forecasting,  and  in 
evaluating  forecasting  techniques. 


TABLE  A . 1 


EVALUATION  OF  GW C  FORECASTS  (PHASE  III) 


Vertical  ms 

Gusts 

Terrain 

Input 
Terms 
(L  or  S) 

Cr  ..relation 
Coef  f iclents 

Overall 

Error 

(X) 

Error  Above 
Threshold 
(X) 

rms 

Threshold 

(fps) 

High  Mountains 

L 

-0.0505 

40.00 

34.18 

mm 

Low  Mountains 

S 

0.1981 

40.82 

46.36 

4.0 

Desert 

L 

0.3039 

47.36 

47.36 

2.5 

Plains 

L 

0.6285 

31.00 

20.00 

4.0 

Lataral  raa 

Guata 

Terrain 

Input 

Tama 

(L  or  S) 

Correlation 

Coefficient! 

Overall 

Error 

(X) 

Error  Above 
Thraahold 
(X) 

me 

Thraahold 

(fpa) 

High  Mountain* 

S 

0.0447 

30.21 

19.07 

mm 

Low  Mountain* 

L 

0.2974 

30.63 

16.81 

Daaert 

L 

0,4581 

30.13 

29.79 

mm 

Plain* 

L 

0.6087 

27.26 

14.77 

4.0 
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TABLE  4.2 


EVALUATION  OF  SHOWALTER  INDEX  (PHASE  III) 


Vertical  rms 

Gusts 

Terrain 

Input 
Terms 
(L  or  S) 

Correlation 

Coefficients 

Overall 

Error 

(%) 

Error  Above 
Threshold 
(%) 

rms 

Threshold 

(fps) 

High  Mountains 

S 

-0.0865 

39.80 

35.14 

4 . 5 

Low  Mountains 

L 

0.1372 

43.42 

45.27 

4.0 

Desert 

L 

-0.0428 

57.93 

57.93 

2.5 

Plains 

S 

-0.2964 

41.06 

32.43 

4.0 

Lateral  rms 

Gusts 

Terrain 

■a 

Correlation 

Coefficients 

Overall 

Er  ror 
(X) 

Error  Above 
Threshold 

(X) 

rms 

Threshold 

(fps) 

High  Mountains 

s 

-0.0657 

30  11 

35.14 

4.5 

Low  Mountains 

L 

0. 2023 

31.02 

49.15 

4.0 

Dasart 

L 

-0.0663 

31.98 

31.98 

2.5 

Plains 

S 

-0.3384 

34.33 

36.27 

4.0 
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'.ABLE  4.3 


EVALUATION  OF  RICHARDSON  NUMBERS  (PHASE  III) 


Vertical  rns  Gusts 

Terrain 

Input 

Terms 

(L  or  S) 

Correlation 
Coef  f icients 

Overall 

Error 

(2) 

Error  Above 
Threshold 
(2) 

rm a 

Threshold 

(fps) 

High  Mountains 

L 

-0.1790 

33.73 

44.25 

4.5 

Low  Mountains 

L 

-0.0540 

44.08 

46.02 

4.0 

Desert 

L 

-0. 1508 

56.97 

34.29 

2.5 

Plains 

S 

-0.3344 

44.59 

36. 0< 

4.0 

Lateral  rms  Gucts 

Terrain 

Input 
Terms 
(L  or  S) 

Correlation 

Coefficients 

Overall 

Error 

(2) 

Error  Above 
Threshold 
(2) 

rms 

Threshold 

(fps) 

High  Mountains 

L 

-0. 1683 

26.39 

34.62 

4.5 

Low  Mountains 

L 

-0,0726 

31.38 

37.96 

4.0 

Desert 

L 

-0.1206 

31.89 

31.39 

2.5 

Plains 

L 

-0.3292 

36.49 

34.81 

4.0 

McConnell  Air  Force  Base,  Kansas 


Turbulence  forecasting  af  McConnell  Air  Force  Base  was  accomplished  using  guidelines  of  the  Air 
Weather  Service  Manual  55-8,  Volume  l  .The  forecasts  were  primarily  based  upon  the  intensity  of 
the  surface  winds,  the  vertical  wind  shear  from  the  surface  to  1,000  feet  above  the  ground,  and 
the  type  of  convective  clouds  observed  or  forecast  for  the  period  of  a  given  flight.  For  wind  shears 
of  3  to  5  knots  per  1 ,000  feet,  light  turbulence  was  forecut;  6  to  9,  moderate;  and  10  or  higher, 
severe.  For  surface  winds  greater  than  IS  knots  but  leu  than  25,  light  turbulence  was  forecast  If 
the  air  was  neutrally  stable  or  unstable.  Similarly,  moderate  turbulence  wu  forecast  for  winds 
greater  than  25  knots.  When  cumuliform  clouds  were  forecast,  the  following  turbulence  waa 
forecast  in  that  vicinty:  fairweather  cumulus,  light  turbulence;  thunderstorm*  or  towering 
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cumulus,  moderate;  mature  or  rapid'y  growing  thunderstorms,  *cvere;  and  severe  thunderstorms, 
extreme  turbulence.  The  forecast  turbulence  level  compared  to  the  pilot's  estimate  of  the 
turbulence  (Tables  4.4  through  4.7)  indicates  that  of  the  methods  used,  wind  shear  (Table  4.4) 
was  the  test  indicator.  The  tendency  to  underestimate  at  lower  turbulence  levels,  and  ovemimate 
at  higner  levels,  is  evident. 


TABLE  4 


TURBULENCE  FORECAST  USING  WIND  SHEAR  ( SURFACE- 1 000 ' ) 


FORECAST 

OBSERVED  BY  PILOT 

NONE 

LIGHT 

MODERATE 

SEVERE 

NONE 

1 

12 

4 

0 

LIGHT 

4 

2 

0 

0 

MODERATE 

2 

4 

1 

0 

SEVERE 

3 

4 

0 

0 

TABLE  4.5 

TURBULENCE  FORECAST  USING  WIND  SPEED 


FORECAST 

OBSERVED  BY  PILOT 

NONE 

LIGHT 

MODERATE 

NONE 

14 

19 

4 

LIGHT 

0 

6 

1 

MODERATE 

0 

0 

0 
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TABLE  4.6 


TURBULENCE  FORECAST  USING  CONVECTIVE  CLOUD  TYPE 


FORECAST 

OBSERVED  BY  PILOT 

LIGHT 

MODERATE 

LIGHT 

3 

1 

MODERATE 

3 

0 

TABLE  4.7 


TURBULENCE  FORECAST  CONSIDERING  MAXIMUM 
VALUE  FROM  WINS  SHEAR,  SURFACE  VJIND,  OR 
CONVECTIVE  CLOUD  TYPE. 


FORECAST 

OBSERVED  BY  PILOT 

NONE 

LIGHT 

MODERATE 

SEVERE 

NONE 

1 

12 

13 

0 

LIGHT 

4 

4 

1 

0 

MODERATE 

5 

6 

1 

0 

SEVERE 

3 

4 

0 

0 

Edwards  Air  Force  Base,  California 

Turbulence  occurrences  from  leg  to  leg  for  a  given  flight  varied  considerably.  During  the  period 
from  8  October  1968  to  8  January!  969,  atmospheric  conditions  in  the  area  were  characterized  as 
stable,  except  for  brief  occurrences  of  post-frontal  instability  ;  however,  showers  and  low  clouds 
over  tne  mountains  usually  prevented  flights  during  these  times. 

The  most  useful  subjective  forecasting  aid  for  turbulence  over  the  flight  legs  was  wind  direction 
and  speed  in  relation  to  the  terrain  covered  by  each  leg.  Particularly  useful  in  determining 
low-level  wind  conditions  from  the  surface  to  about  2,000  feet  above  the  terrain  were  upper-air 
soundinp  from  Vandenberg  Air  Force  Base,  Oxnard  Air  Force  Base,  and  Edwards  Air  Force  Base, 
and  surface  observations  from  Sandberg,  California'.  Sandberg  was  the  only  observing  station 
which  reported  wind  velocities  in  the  mountains, 
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The  turbulence  occurrences  appeared  to  be  highly  related  to  terrain.  The  pilot  observed  that  the 
most  severe  encounters  occurred  near  the  higher  peaks  of  legs  3  and  5.  For  this  reason,  the 
subjective  techniques  are  categorized  by  flight  leg. 

c  Legs  1  and  8:  The  pilot  reported  that  any  low  altitude  wind  speed  over  10  knots 
normally  produced  at  least  light  turbulence  over  these  legs.  A  low-altitude  wind 
direction  from  west-southwest  through  north-northwest  or  from  east-northeast  through 
south-southeast  with  speeds  in  excess  of  25  knots  resulted  in  moderate  or  greater 
turbulence  On  several  occasions,  turbulence  was  encountered  over  the  desert  to  the 
southeast  of  the  mountains  in  an  apparent  lee  wave  effect. 

•  Legs  2  and  7:  Turbulence  was  reported  on  only  a  few  fill  'its  over  these  legs.  The  low 
altitude  flight  path  above  the  valley  floor  was  apparent!"  somewhat  sheltered,  even 
from  possible  effect!  of  strong  northerly  winds  down  the  valley. 

•  Legi  3  and  5;  The  most  severe  turbulence  encounters  occurred  on  these  legs.  From  a 
forecut  standpoint,  any  low  altitude  wind  speed  over  10  knots  resulted  in  reports  of  at 
leut  light  turbulence  by  the  pilot  over  one  of  the  legs.  A  low  altitude  wind  direction 
from  northwest  through  northeast  or  from  southwest  through  southeast,  accompanied 
by  speeds  In  excels  of  25  knots,  resulted  in  moderate  or  greater  turbulence.  An 
Interesting  feature  is  that  the  reported  turbulence  occurred  in  a  different  position 
relative  to  the  terrain,  depending  on  the  wind  direction.  The  position  of  the  hardest 
strikes  was  generally  to  the  lee  of  a  m^jor  terrain  feature. 

•  Leg  6  :  Turbulence  v>as  reported  by  the  pilot  on  this  leg  leu  often  than  on  legs  3  and  5, 
In  general,  low  altitude  wind  speeds  of  20  knots  or  more  produced  at  least  light 
turbulence  on  this  leg. 

•  L«8  4  :  Turbulence  was  rarely  reported  on  this  leg  over  the  coastal  water.  Occasions  on 
which  light  turbulence  wu  reported  over  this  leg  were  asaoclated  most  often  with  a 
strong  offshore  preiaure  gradient  and  strong  Santa  Ana  wind  channeled  down  the  Santa 
Gar*  River  Valley. 

Major  emphasis  was  placed  on  the  low  altitude  wind  from  the  surface  to  about  2,000  feet  above 
the  aurface,  aince  this  appeared  to  be  a  key  factor  related  to  reports  of  turbulence.  From  a 
forecut  standpoint,  low  altitude  wind  In  both  direction  and  speed  must  be  forecut  over  the  flight 
area  from  available  soundings  and  obiervations  prior  to  the  turbulence  forecast.  Understanding  of 
synoptic  conditions  and  resulting  wind  conditions  is  useful  in  this  regtrd,  but  this  was  not 
attempted  since  It  involves,  enentially,  a  case  study  spproach.  Given  an  accurate  low  altitude  wind 
velocity  forecut  for  the  seasonal  period  of  Phase  III  LO-LOCAT  flights,  one  can,  through  terrain 
considerations,  prepare  a  subjective  low  altitude  turbulence  forecut, 

Objective  Forecuts _ 


The  objective  forecut  models  being  developed  are  based  on  mu’tiple  regression  analysis  (Section 
V  and  Appendix  III).  Because  most  models  developed  use  the  gust  velocity  rms  as  the  dependent 
variable,  an  analysis  of  maximum  guit  velocities  wu  conducted  to  determine  if  a  turbulence  index 
based  on  msximum  gust  velocity  could  prove  more  significant.  It  was  found  that  the  maximum 
gust  velocities  occurred  frequently  over  the  same  legs  of  a  particular  course.  This  indicstes  that 


62 


some  local  terrain  feature  or  combination  of  terrain  feature  and  wind  pattern  was  responsible  for 
the  maximum  gusts.  For  instance,  legs  7  and  8  at  Edwards  AFB  (Phases  I  and  II  routes)  were 
subject  to  highly  localized  lee  wave  turbulence  under  certain  stability  and  wind  conditions, 
whereas  the  remaining  legs  would  have  light  or  no  turbulence  Therefore,  only  the  gust  velocity 
rms  values  are  being  evaluated  to  develop  forecast  models.  It  was  also  determined  that  the  lower 
rms  values  are  present  so  often  that  they  constitute  a  form  of  “noise”,  and  as  such  correlate 
poorly  with  the  forecast  data.  In  a  case  study  of  twelve  of  the  highest  rms  values  for  each  baae, 
for  example,  it  was  found  that  variables  such  as  vertical  wind  shear  and  Richardson  number  were 
better  indicators  of  turbulence  than  when  all  of  the  data  were  used  in  this  regression  analyst. 
Therefore,  it  was  decided  to  eliminate  the  lower  gust  velocity  rms  vaiuea  from  the  Phases  I  and  II 
data.  Analysis  indicated  that  2  40  fps  was  a  suitable  cutoff  point  for  all  but  the  water  leg.  The 
small  amount  of  data  precluded  the  possibility  of  eliminating  the  lower  gust  velocity  rms  values 
for  this  leg.  All  the  Phase  Ill  data  will  be  retained  since  there  are  fewer  data  than  was  the  case 
with  Phases  1  and  II. 

The  input  variables  used  in  the  regression  analysis  are  tabulated  in  Appendix  III. 

CASE  STUDIES,  PHASES  1  AND  11 

Case  Studies  of  Phases  I  and  II  data  have  been  completed  for  the  twelve  highest  vertical  gust 
velocity  rms  values  from  each  of  the  four  bases.  Although  rms  gust  velocities  have  quite  a  limited 
range,  the  maximum  vertical  gust  velocities  for  these  flights  range  from  about  16  to  50  fps,  with 
an  average  of  about  27  fps.  No  consistent  synoptic  pattern  could  be  discerned. 

The  most  striking  features  noted  were  the  vertical  wind  shears.  The  low-mountain  legs  had  an 
average  wind  shear  of  3.29  knota  per  100  feet,  for  high  mountains  and  plains,  the  values  are  1 .94 
and  1 .53  knota  per  100  feet,  respectively. 

The  next  most  consistent  feature  is  the  low  value  of  the  Richardson  numbers.  For  the  low 
mountain  legs,  the  Richardson  number  (R)  never  exceeded  0.015.  The  R  over  the  high  mountain 
lega  are  small  negative  values,  with  the  exception  of  one  case  where  It  ia  plus  0/9.  Over  the  plains 
there  is  less  consistency:  there  Is  a  high  value  of  2,42,  three  values  between  one  and  two,  one  value 
vetween  one-half  and  one,  and  seventeen  values  lust  than  one-half. 

'he  mountain  legs  have  higher  vertical  gust  velocity  average  rms  values  than  do  the  plain*  lep 
(high  mountains  6.59  fps,  plains  4.94  fps).  The  lowest  gust  velocity  rms  value  over  the  mountains 
it  also  higher  than  the  average  for  the  plains.  These  same  relationships  also  hold  true  for  the 
maximum  vertical  gust  velocities  (average,  high  mountains  33.49  fps;  low  mountains  32.31  fps; 
pli  ins  2 1 .8 1  fps). 

T hire  is  some  tendency  for  these  high  vertical  gusts  to  occur  during  the  spring  season.  Over  the 
p'ains,  this  is  true  for  16  out  of  26  cases;  over  the  low  mountains,  it  is  true  for  6  out  of  14  cases. 

A  ciumal  trend  appears  only  foi  the  plains  lep:  mid  morning  and  afternoon  flights  account  for 
2*  Dut  of  the  26  cases.  This  could  be  the  result  of  convective  activity,  especially  since  the  sine  of 
the  solar  elevation  is  greater  than  0,55  for  these  24  cases.  However,  the  high  wind  shears  recorded 
for  '.hese  cases  and  the  fact  thai  unstable  lapse  rates  are  in  the  minority  indicate  that  convective 
activity  is  not  responsible  for  the  diurnal  trend, 
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It  is  probably  a  coinciceme  that  the  ratio  of  the  average  wind  speeds  for  the  mountains  and  plains 
is  fairly  close  (1.40  versus  1  33)  to  the  ratio  of  the  respective  average  gust  velocity  rms  values. 
Over  the  plains,  the  range  of  v/ind  speeds  (3  to  35  knots)  combined  with  the  low  average  value  (15 
knots)  suggest  that,  fcr  tnese  flight  samples  at  least,  wind  speed  by  itself  is  not  causally  related  to 
high  rms  gusts,  Over  the  mountains,  the  average  wind  speeds  (21  knots)  is  submarginal  for 
producing  full  scale  lie  waves;  however,  this  is  high  enough  to  produce  mechanical  turbulence. 
The  higher  winds  (up  tc  30  knots)  could  of  course  produce  lee  waves  given  the  proper  stability 
profile  and  terrain  orientation. 

An  example  case  study  follows.  This  case  was  chosen  because  five  out  of  a  total  of  twelve  large 
gust  velocity  rms  values  occurred  on  the  same  day.  In  addition,  a  maximum  gust  velocity  of  50 
fps  was  recorded. 

Edwards  Air  Force  Eax  (13  December  1967, 

General  Synoptic  Situation 

An  extensive  and  intense  high  pressure  area  was  centered  over  Bntish  Columbia  with  a  ridge  line 
extending  southward  into  California.  A  low  pressure  area  was  located  south  of  the  flight  route  as 
shown  in  the  surface  analysis  for  1500  GMT  in  Figure  4.21 .  This  synoptic  pattern  is  frequently 
associated  with  Santa  Ana  winds,  (relatively  strong,  warm  and  dry'  off-shore  winds  observed  along 
the  southern  California  coast).  A  fast-moving  cold  front  had  passed  through  the  urea  the  previous 
night,  bunging  very  dry  air.  The  cold  front  then  became  stationary  for  several  hours  near  the 
Mexican  border  (see  Figure  4.21).  Throughout  the  flight  route  area  (outlined  by  dotted  lines)  the 
surface  pressure  gradient  gradually  increased  and  therefore,  the  winds  were  expected  to  increase. 
The  surface  winds,  which  varied  in  direction  from  northwest  to  northeast,  were  considerably 
Influenced  by  topography  The  temperature  pattern  isotherms  of  5*  C  are  represented  by  the 
dashed  lines  and  was  similarly  strongly  influenced  by  topography.  Figure  4.22  shows  the  surface 
conditions  at  1800  GMT. 

The  synoptic  patterns  at  850  mb  for  0000  and  1200  GMT  are  shown  in  Figure  4.23  and  4.24 
where  the  pressure-height  gradient  is  not  unusual.  The  flow  as  represented  from  vertical  soundings 
U  northeasterly  in  contrast  to  the  more  easterly  flow  suggested  by  the  pressure  gradient.  At  this 
height  (about  2500  feet  above  lowest  ground),  the  effects  of  terrain  are  still  present.  Low  pressure 
at  both  the  850  and  700  mb  "evels  was  centered  over  San  Diego, 

Dawn  Flight 

Turbulence  was  encountered  over  legs  6,  7  and  8,  wuh  maximum  vertical  gusts  of  50  fps  over  legs 
7  and  8.  Neutrally  stable  air  and  scattered  clouds  existed  over  the  route.  Winds  at  flight  altitude 
were  20  knots  over  leg  6  and  28  knots  over  leg  7,  Legs  6  and  7  lie  over  the  Santa  Clara  River 
valley  directly  in  the  lee  of  the  coastal  mountains.  Without  consulting  a  radiosonde  sounding,  it 
cannot  be  known  that  lee  waves  were  present  in  the  valley.  However,  the  flow  at  850  mb, 
corresponding  approximately  to  the  heights  of  tb"  mountains,  was  about  35  knots.  This  is  more 
than  adequate  for  the  generation  of  lee  waves  if  a  stable  layer  is  present  at  higher  levels.  This 
could  explain  the  difference  of  20  fps  in  the  maximum  gust  velocity  over  legs  6  and  7.  Leg  S  lies 
partly  In  the  same  valley  so  that  lee  waves  could  also  be  experienced  there.  Leg  8  also  crosses  part 
of  the  mountain  range  making  it  vulnerable  to  lee  effects  from  some  of  the  higher  isolated  peaks. 
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The  win^  shears  were  1.4  knots  per  100  feet  over  leg  6  and  1.3  over  leg  7.  The  Richardson 
numbers  were  small  and  negative 

Mid  Morning  Flight 

Turbulence  was  encountered  over  legs  7  and  8,  with  maximum  vertical  gusts  of  50  fps  on  leg  7. 
Unstable  lapse  rates  were  recorded.  Winds  at  flight  altitude  were  22  knots  over  leg  7  and  28  knots 
over  leg  8.  Wind  shears  were  2.6  knots  per  100  feet  over  leg  7  and  2.3  over  leg  8.  The  Richardson 
numbers  were  small  and  negative. 

CASE  STUDIES,  PHASE  111 


Edw  ards  Air  Force  Base  ( 1 3  November  1 968) 


General  Synoptic  Situation 

Figures  4.25  and  4.26  show  the  surface  synoptic  analysis  for  1500  and  1800  GMT,  respectively. 
Through  this  three-hour  time  period  in  which  the  flights  were  conducted,  very  little  variation  was 
noted  in  the  synoptic  pattern.  A  large  high  pressure  system  was  centered  in  the  eastern  Pacific 
Ocean  at  37*  N  and  135"  W  from  which  a  broad  ndge  extends  eastward  into  Oregon  and  towards 
Colorado. A  minor  ridge  was  observed  to  develop  through  central  California  and  south  towards 
Yuma.  Arizona.  The  effects  of  this  ndge  were  noticed  in  a  generally  high  pressure  gradient 
through  the  flight  area.  Winds  through  this  area  were  mostly  from  the  north  to  east  and  increased 
through  the  time  period  from  8  to  10-15  knots.  The  temperature  pattern  shown  by  the  dashed 
isotherms  was  generally  influenced  by  the  topographical  variation  in  surface  heights. 

Aloft  at  850  mb  for  0000  and  1200  GMT  as  indicated  in  Figures  4.27  and  4.28  a  low  pressure 
center  occurred  over  the  southern  California  border.  This  was  associated  with  a  surface  low 
pressure  cell  and  frontal  system  which  extended  through  Arizona.  Winds  at  850  mb  through  the 
flight  area  were  north  to  northeast  During  the  12-hour  period  presented,  little  change  was 
observed  in  the  pattern  of  pressure-heights  through  the  flight  area.  At  500  mb,  a  broad  trough 
having  a  north-south  axis  along  the  coastal  states  intensified  during  the  12-hour  period  examined. 
Such  development  is  generally  reflected  in  a  slow-down  of  movement  of  surface  systems. 

Mid-Morning  Flight 

Stable  and  relatively  dry  atmospheric  conditions  prevailed  over  the  flight  legs  in  Test  No.  87,  The 
pilot  reported  moderate  to  servere  turbulence  over  legs  3  and  5  with  the  hardest  strikes  occurring 
to  the  lee  of  major  terrain  features,  The  north  winds,  in  some  local  cases  exceeding  25  knots  at 
flight  altitudes,  are  approximately  crossing  the  coastal  mountains.  Considerable  variation  was 
observed  between  the  wind  directions  and  speeds  observed  over  individual  legs,  and  it  was 
concluded  that  this  variation  is  closely  correlated  with  the  resulting  turbulence.  Legs  3  and  5  have 
the  largest  wind  speed  components  along  or  against  the  wind,  respectively,  which  when  combined 
with  the  cross-ridge  effects  contributes  to  the  resulting  rms  gust  velocities  observed  of  1 1 .46  and 
8.50  feet  per  second.  This  compares  with  rms  gusts  less  than  3.5  feet  per  second  in  the  other  legs 
of  this  flight.  It  is  interesting  to  note  that  flight  along  leg  4  resuited  in  the  lowest  rms  gust  velocity 
of  this  test  (1.24  fps).  This  leg  is  essentially  parallel  to  and  in  the  lee  of  the  coastal  ridge  but  out 
over  the  water.  It  must  be  concluded  that  the  terrain  perturbations  of  the  winds  do  not  extend 
i.or  travel  far  downwind  under  this  synoptic  pattern. 
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OBSERVED  MINUS 


O  Phases  I  and  II  Data  (All  Seasons) 

□  Phase  III  Data 

o phases  I  and  II  Data  (Fall  Season  Only) 
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Figure  lf.U  Wind  Speed  Cumulative  Probability  at  Edwards 
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igure  t.li  Variation  of  'Wind  Speed  Cumulative  Probability 
over  High  Mountains  with  Time-of-Day 
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figure  4.13  Variation  of  Wind  Speed  Cumulative  Probability 
ove:  Desert  with  Time-of-Day 
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CUMULATIVE  PROBABILITY 
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Figure  4.17  Longitudinal  Gust  Velocity  RMS  Cumulative  Probability 
Associated  vith  Richardson  Number 
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SECTION  V 


PRELIMINARY  TURBULENCE  ANALYSIS 

Preliminary  turbulence  analysis  was  accomplished  on  data  obtained  dunng  the  first  half  of  the 
Low-Level  Critical  Air  Turbulence  Phase  III  program.  Approximately  72  hours  of  low-level 
(0-1000  feet)  turbulence  and  associated  meteorological  data  were  recorded  from  16  August  , 
1968,  through  8  January,  1  969.  The  data  were  obtained  during  flight  over  routes  near  McConnell 
Air  Force  Base.  Kansas  and  Edwards  Arr  Force  Base,  California  Also  included  are  time  histones 
of  some  of  the  iarger  magnitude  gust  velocities  recorded  over  the  Peterson  Field,  Colorado  route. 

Dunng  the  LO-LOCAT  Phase  III  program,  gust  velocities  are  not  being  computed  for  data  samples 
where  the  level  of  turbulence  is  low,  which  reduces  computer  effort  on  data  which  would 
eventually  be  discarded  due  to  a  low  signal-to-noise  ratio  The  standard  deviations  of  the  angles 
of  attack  and  sideslip  differential  pressures,  as  sensed  by  the  gust  probe  pressure  ports,  were  the 
criteria  used  in  determining  which  samples  fall  into  this  category.  Data  were  not  processed 
when<k<0.07  in.  Hg  or  c^<0.05  in.  Hg  For  this  situation,  the  standard  deviations  of  gust  velocities 
were  nearly  always  less  than  1 .0  fps,  The  values  of  cra  and  are  determined  from  the  calibration 
program  (see  Appendix  III).  When  these  small  values  are  obtained,  then  all  of  the  data  processing, 
including  plotting  that  involves  gust  velocities  is  halted  early  in  the  data  processing  sequence. 
Thus,  this  criteria  results  in  a  savings,  while  assuring  processing  of  the  meaningful  turbulence  data. 
The  manner  in  which  this  low  turbulence  level  data  are  included  in  the  various  statistical  analyses 
Is  explained  in  each  appropriate  section, 

Three  different  types  of  counting  techniques  are  being  used  to  evaluate  gust  velocity  data  These 
techniques  are  the  peak  count,  amplitude  count,  and  level  crossings  count.  A  description  will  be 
given  of  each  of  these  techniques  hollowing  the  descriptions,  results  will  be  presented  and 
comparisons  made. 

PRIMARY  PEAK  COUNT 

In  the  peak  count  procedure  each  maximum  positive  or  negative  excursion  between  adjacent 
crossings  of  the  mean  is  defined  as  a  primary  peak  As  shown  in  Figure  5,1,  secondary  peaks  are 
not  considered.  The  number  of  primary'  peaks  which  fall  in  each  2  fps  wide  amplitude  band  are 
counted. 

In  the  peak  count  procedure,  high  frequency,  low'  amplitude  data  (or  noise)  occu.  as  primary' 
peaks  in  the  lowest  band  Figure  5.2  shows  a  peak  count  distribution  and  illustrates  how  it  was 
extrapolated  to  eliminate  the  undesirable  peaks  in  the  lowest  band.  The  shaded  points  indicate 
the  cumulative  number  of  peaks  in  the  lowest  band  as  they  ware  originally  counted.  The 
unshaded  points  are  the  extrapolated  values.  This  curvilinear  extrapolation  is  performed  based  on 
the  shape  of  'he  distribution  curve  This  is  done  by  fitting,  by  computer,  a  curve  through  the  data 
pomts  and  extrapolating  it  to  zero  gust  velocity.  The  new  value  of  N [>,  and  f;  obtained  by  this 
extrapolation  is  then  used  in  computation  of  the  standard  deviation  on 
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(5.1) 


Also,  Npe  is  u  1  for  calculation  of  the  cumulative  probability . 


AMPLITUDE  COUNT 


The  amplitude  count  technique  of  amplitude  analysis  was  performed  on  computed  values  of  the 
gust  velocity  components  Each  component  was  computed  at  100  samples  per  second.  In  this 
procedure,  amplitude  bands  2  fps  wide  are  placed  on  either  s*de  of  the  mean.  The  number  of 
samples  which  fall  in  each  band  is  determined,  then  corresponding  positive  and  negative  bands  are 
added  together.  This  procedure  is  illustrated  in  Figure  5.3.  These  data  are  used  to  compute 
cumulative  distributions,  cumulative  probabilities,  and  probability  density  distributions  The  time 
senes  standard  deviation  (?t)  is  computed  using  the  following  equation 

r  * 

"t  "  [  T  E  <  *•  '  *  J  (5  2 


LEVEL  CROSSINGS  COUNT 


Another  method  of  amplitude  analysis  utilizes  the  level  crossing  technique.  I  evels  of  gu-it  velocity 
are  established  at  2  fps  intervals  from  the  zero  mean.  Only  crossings  with  a  posi'.h  slope  srt 
counted.  Figure  5.4  illustrates  the  level  crossings  technique.  Crossings  of  correspond iru;  posiU'e 
and  negative  levels  are  combined  and  zero  'evel  crossings  are  doubled.  The  number  of  leval 
crossings  per  mile  of  the  level  x(t)  =o  is  computed  for  each  level  crossed  by: 


N 


(3.3) 


A  level  crossings  standard  deviation  (<?,)  is  calculated  using  the  following  equation: 


9  L  ■ 


(5  4) 


COMPARISON  OF  COUNTING  TECHNIQUES 

The  three  counting  techniques  used  to  analyze  gust  velocity  time  series  are  compared  In  the 
following  paragraphs.  Only  data  obtained  during  Phase  111  o-  *,  the  McConne'l  and  Edwards  route* 
were  considered.  For  these  data,  individual  turbulence  tamples  were  originally  claaaified  into 
categories  as  discussed  in  Section  II.  Terrain,  altitude,  and  atomospheric  stability  were  found  to  be 
the  primary  categories  and  all  samples  were  grouped  into  them  for  this  analyaia, 
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It  should  be  noted  that  the  effects  of  removing  samples  where  turbulence  level  was  low  has  not 
been  taken  into  account  for  these  data  presented  in  this  report.  A  computer  program  is  being 
modified  to  incorporate  these  samples  for  further  analysis. 

Cumulative  distributions  of  peaks  per  mile  and  kvel  crossings  per  mile  were  calculated  for 
comparison.  Figures  5.5,  5.6,  and  5.7  show  comparisons  of  the  occurrences  per  mile  of  peak 
count  end  level  crossings  for  tne  three  gust  velocity  components  of  the  all-data  category 
(McConnell  and  Edwards).  These  figures  show  that  the  occurrences  per  mile  for  level  crossing  arc 
always  higher  than  for  peaks,  except  at  the  highest  gust  velocities.  If  the  peak  count  data  had  not 
been  extrapolated,  these  distribution  curves  would  have  agreed  at  zero  gust  velocity.  The  reason 
for  extrapolating  is  explained  in  the  peak  count  discussion. 

Probabilities  of  encountering  a  value  equal  to  or  greater  than  a  given  value  of  gust  velocity  for 
the  three  counting  techniques  are  shown  for  the  three  gust  velocity  components  in  Figures  5.8, 
5.9,  and  5.10.  The  category  represented  in  these  figures  includes  all  McConnell  and  Edwards  data. 
The  relationship  of  the  probability  distributions  of  these  three  counting  techniques  varies 
somewhat  from  one  gust  component  to  another.  However,  in  general,  the  amplitude  count 
distributions  are  lowest  and  the  level  crossing  distributions  are  the  highest.  The  probability 
distributions  for  the  peak  count  technique  are  consistently  lower  than  the  level  crossing  technique 
except  at  the  higher  gust  velocities  where  they  are  higher. 

The  amplitude  and  peak  count  distributions  given  in  the  LO-LOCAT  Phases  I  and  11  final  report 
were  more  nearly  a  straight  line  than  those  given  here.  The  Phase  III  distributions  now  being 
calculated  are  curving  up  and  showing  higher  probabilities  of  encountering  given  gust  velocities  at 
the  larger  gust  velocities.  The  reason  for  this  is  believed  to  be  due  to  the  increased  speed  of  the 
Phase  III  test  aircraft  which  permitted  measurement  of  the  long  ;r  wavelength  higher  amplitude 
gusts.  These  gusts  have  the  effect  of  making  the  distribution  cur;  turn  up  slightly  at  the  higher 
gust  velocity  values. 

Since  the  relationship  between  peak  count,  amplitude  count,  and  level  crossings  has  already  been 
shown,  only  peak  count  data  will  be  shown  to  illustrate  the  effect  of  different  geophysical 
conditions.  Figures  5.11,  5.12,  and  5.13  show  the  terrain  effects  on  peak  count  distributions.  The 
relationship  of  these  distributions  are  similar  to  those  for  LO-LOCAT  Phases  1  and  11  data  except 
that  the  difference  between  high  and  low  mountain  data  is  more  pronounced.  These  plots  show 
that  the  number  of  peaks  per  mile  exceeding  a  given  value  is,  in  general,  larger  for  high  mountain 
data  and  decreases  for  low  mountain,  plains,  desert,  and  water,  in  that  order. 

Figures  5.14,  5.15,  and  5.16  show  the  altitude  effects  on  peak  count  distributions.  These  figures 
indicate  a  higher  frequency  of  occurrence  of  gust  velocity  peaks  at  250  feet  altitude  than  at  750 
feet.  Here  again  the  relationship  of  these  distributions  agrees  well  with  those  for  Phases  1  and  II. 

Figures  5.17,  5.18,  and  5.19  show  the  effect  of  atmospheric  stability  on  peak  count  distributions. 
The*  plot*  show  that  in  general  very  stable  data  have  the  lowest  distribution,  then  stable  and 
neutral,  in  that  order.  The  distribution  of  unstaole  data  has  a  much  steeper  slope  than  the 
distributions  of  the  other  stabilities.  Unstable  data  did  not  show  this  characteristic  for  Phases  land 
II  data.  It  is  possible  that  this  relationship  will  change  when  the  Peterson  Field  high  mountain  data 
are  included. 


Plots  showing  the  relationship  of  the  three  components  of  guit  velocity  for  all  three  counting 
techniques  are  given  ir  Appendix  V. 

GAUSSIAN  DISTRIBUTION  CHECK 


1 

1 


Amplitude  count  d  ;ta  obtained  for  each  turbulence  sample  were  checked  for  normality  using  the 
chi-squared  *ood  iess-of-fit  test.  This  test  compares  the  actual  distribution  with  a  normal 
distribution. 

A  chi-squared  test  was  performed  on  each  four  and  one-half  minute  turbulence  sample  using  28 
bends  (Nv)  with  widths  of  constant  probability.  The  degrees  of  freedom  for  each  sample  were 
(Nb-1).  Two  hundred  gust  velocity  amplitude  values  were  used  ir.  the  calculation: 

2  8 

x2  -  <iu  2  **  *  200  (5,5) 

b«l 


A  cumulative  probability  distribution  of  the  chi-square  values  for  all  McConnell  and  Edwards  data 

is  shown  in  Figure  5.20.  Levels  of  significance  based  on  the  number  of  degrees  of  freedom  for  / 

these  data  are  shown.  Inspection  of  this  fi;  ire  shows  that  approximately  83  percent  of  the 

vertical,  84  percent  of  the  longitudinal  and  88  percent  of  the  lateral  5ust  velocity  samples  were 

accepted  as  Gaussian  at  the  0  02  level  of  significance.  Approximately  10  percent  more  samples 

were  accepted  dunng  Phases  I  and  II.  This  is  attributed  to  the  fact  that  turbulence  samples  arc 

being  taken  over  a  60  percent  longer  distance  dunng  Phase  III. 

RUN  TEST  OF  STATION AR1TY 

To  thoroughly  analyze  clear  air  turbulence  data,  the  proof  of  nationality  is  required.  When  a 
single  time  hsitory  record  is  refereed  to  as  being  stationary,  thus  generally  means  that  statistical 
properties  computed  over  short  time  intervals  do  not  vary  more  than  would  be  expected  due  to 
normal  sampling  variatons.  Only  a  limited  analysis  of  stationarity  was  accomplished  prior  to  this 
interim  report,  however,  results  thus  far  have  shown  reasonable  agreement  with  homogeneity 
checks.  Run  test  results  will  be  presented  in  the  final  report. 

COMPARISON  OF  STANDARD  DEVI  AT  IONS 

Gust  velocity  standard  deviations  from  amplitude  count,  yeak  count,  and  level  crossings  were 
calculated  as  shown  by  equations  5  1,  5.2,  and  5.4,  Figures  5.21,  5.22,  and  5.23  show 
comparisons  of  the  cumulative  distributions  of  these  standard  deviations  for  the  three  gust 
velocity  components.  AU  McConnell  and  Edwards  data  are  included.  Time  series  standard 
deviations  show  a  lower  probability  of  occurrence  than  peak  count  or  level  crossing  stanaard 
deviations. 

CHARACTERISTIC  FREQUENCY 

It  should  be  noted  that  the  characteristic  frequency,  N0,  3S  refereed  to  in  this  report  is  not  the 
actual  characteristic  frequency  of  the  turbulence  but  rs  equivalent  to  the  intercep*  of  the  peak 
count  cumulative  distribution  function.  Measurement  of  the  actual  characteristic  frequency  of 
the  turbulence  would  require  calculation  out  to  sery  high  frequencies  (around  30.000  cps). 
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Vglue*  of  the  characteristic  frequency  were  obtained  by  two  different  method*.  In  the  first 
method,  N&iin  cycles  per  foot-'  was  obtained  from  PSD  data  using  equation  5  6. 


N0  ■  /  <r 


(5.6) 


In  the  auvond  method,  Nr  (peaks  per  mile)  was  determined  by  extrapolation  of  the  peak  count 
distribution  curve  as  previously  di^cribed  under  peak  count.  Values  of  were  determined  by  this 
method  for  each  turbulence  sample. 


In  order  to  compare  cliaracteristic  frequencies  calculated  by  both  methods,  N0 values  were 
converted  to  occurrences  pet  mile  and  averaged.  Average  values  for  each  method  are  given  in 
Table  5.1. 


TABLE  5.1 

COMPARISON  OF  CHARACTERISTIC  FREQUENCIES  CALCULATED  FROM 
POWER  SPECTRA  AND  PEAK  COUNT  DATA 


N0  (occurenccs/mile) 

Longitudinal 

Lateral 

Vertical 

26.39 

31.73 

34.72 

Np  (peaWmile) 

23.30 

30.81 

38.36 

These  averages  represent  McConw !'  and  Edwards  data.  Power  spectra  (and  hence,  N  0  values)  were 
calculated  for  approximately  50  percent  of  the  data  processed.  In  addition,  characteristic 
frequency  values  were  not  calculated  for  irregularly  shaped  spectra  (i.e.  non-homogeneous 
turbulence  or  'ow  slgnal-to-noisc  ratio).  Extraneous  NF  values  which  deviated  more  than  four 
standard  deviations  from  the  mean  were  also  removed  before  averaging. 

Characteristic  frequency  for  the  longitudinal  component  is  the  smallest  and  vertical  is  the  largest. 
This  is  true  for  both  N0and  Np.  Comparison  of  corresponding  N0and  NP  values  shows 
reasonable  agreement  between  the  two  methods  of  computation. 

ENSEMBLE  AVERAGING 


Atmospheric  turbulence,  over  long  periods  of  time,  must  be  considered  to  be  nonstationary, 
Stationarity  Implies  that  statistical  properties  computed  over  short  time  invervals  do  not  vary 
more  than  would  be  expected  due  to  normal  sampling  errors.  The  fact  that  turbulence  is 
nonatationary  is  seen  by  comparing  mean  square  values  computed  from  turbulence  samples  for 
different  days.  Also,  in  mountainous  regions  where  the  terrain  is  not  homogeneous  and  the  wind 
currenta  vary  considerably,  mean  square  values  for  adjacent  4. 5-minute  turbulence  samples  have 
been  found  to  vary  from  3  to  75  (fps)2 

Nonstationary  data  are  segmented  into  lengths  which  are  stationary  and  ensemble  averaged.  The 
turbulence  data  are  automatically  segmented  by  the  sampling  routine  of  recording e-.iy  270 seconds 
of  continuous  data.  These  data  segments,  for  the  most  part,  are  stationary  (i.e.  homogeneous); 
however,  only  samples  which  have  been  proven  to  be  self-stationary  are  useu  in  the  ensemble 
averaging  technique. 


98 


Ensemble  averaging  was  accomplished  by  taking  an  ensemble  (group)  of  stationary  record-1  and 
averaging  the  ensemble  at  corresponding  times.  That  is,  individual  points  of  the  ensemble  average 
time  functional)  is  calculated  by; 


fti)  *  ~7  2  *a(tl)  (5.7) 

k  ■  1 

Ensemble  averaging  was  accomplished  for  a  total  of  rventy-thret  4.5-minub  turbulence  samples 
for  tnree  different  categories.  Power  spectra  and  associated  data  (dis:us*ed  later  in  this  report)  for 
category  424000,  consisting  cf  12  samples,  are  presented  in  Figure  5.24.  Data  for  the  o»her 
categories  are  not  presented  due  to  the  limited  number  of  samples  and,  hence,  low  statistical 
confidence. 

Due  to  the  limited  quantity  of  data  analyzed  thus  far,  no  definite  conclusions  can  be  expressed. 
However  certain  consistent  observations  have  be'.n  noted.  These  are:  (1)  As  the  number  of 
ensemble  samples  increases,  the  standard  deviation  of  the  ensemble  average  time  function 
decreases,  (2)  the  scale  lengths  of  the  ensemble  average  time  function  calculated  from  the  power 
spectra  are  approximately  equal  to  the  average  of  the  scale  lengths  for  the  ensemble  records,  (3) 
the  ensemble  average  time  function  is  both  stationary  and  normally  distributed,  and  (4)  the 
ensemble  average  data  are  isotropic,  homogeneous,  incoherent,  and  in  good  agreement  with  von 
Karman  mathematical  spectra. 

SEVERE  AND  EXTREME  TURBULENCE  ENCOUNTERS 

Time  histories  of  longitudinal,  lateral,  and  vertical  gust  velocities  are  being  plotted  for  turbulence 
samples  during  which  the  pilot  reported  severe  or  extreme  turbulence  and  where  the  value  of  one 
or  more  of  the  gust  velocities  exceeded  50  feet  per  second. 

The  above  criteria  were  met  on  one  turbulence  sample  from  the  Edwards  route  and  on  no  samples 
from  the  McConnell  route.  Eighteen  turbulence  encounters  meeting  c  exceeding  the  above 
criteria  occurred  during  flights  at  Peterson  Field  during  the  time  period  from  20  January  1969  to 
15  Feoruary  1969. 

The  time  histories  of  gust  velocities  for  the  single  Edwards  encounter  rre  presented  in  Figure  5.25 
and  those  for  the  Peterson  occurrences  are  shown  in  Figures  5.26  through  5.43.  A  data  Wock 
included  on  each  figure  identifies  the  test  and  leg  during  which  the  severe  turbulence  was 
encountered,  Also  tabulated  in  the  data  block  are  the  maximum  and  minimum  gust  velocities, 
gusto^s,  average  wind  direction  and  velocity,  average  ambient  air  temperatur,  average  true 
airspeed,  average  radar  altitude,  and  category  number. 

The  vertical  gust  velocity  component  in  Figure  5.25  should  be  particularly  noted  because  of  the 
long  wavelengths  of  approximately  1 4,300  feet  and  the  high  magnitudes  of  74.5  fpa  and  -65  fps 
of  the  turbulence.  The  long  wavelength  component  appears  at  a  frequency  of  approximately  0.04 
cps  in  the  time  history  This  same  long  wavelength  turbulence,  if  encountered  during  the  Phases  I 
and  II  testing,  would  heve  appeared  as  0  02  cps  data  due  to  the  lower  speed  at  which  the  C-131 


aircraft  were  down.  Phases  1  and  I!  gust  velocity  measurements  in  this  frequency  regime  would 
have  been  attenuated  approximately  80  percent  due  to  the  low  frequency  data  filtering  techniques 
utilized  to  alleviate  the  effects  of  instrumentation  drift.  This  is  a  cogent  example  of  the  fact  that 
there  is  turbulence  at  low  altitudes  that  contains  wavelengths  longer  than  7000  feet,  and  this  is  the 
reason  why,  in  comparing  Phase  Ilf  results  with  Phases!  and  II  for  equivalent  categories,  that  (1)  the 
probability  distributions  of  peak  gust  velocities  arc  higher,  (2)  the  probability  distributions  for  rms 
gust  velocities  are  higher,  (3)  integral  scale  lengths  are  longer,  (4)  Taylor’s  microscales  are  longer,  and 
(5)  viscous  dissipation  rates  are  less. 

The  striking  similarity  between  Figure  5.33,  5.37,  and  5.41  should  also  be  noted.  These 
turbulence  samples  were  gathered  over  leg  2  of  the  Peterson  route  on  three  separate  flights 
conducted  on  the  same  day.  Leg  2  starts  approximately  one  mile  south-vest  of  Elk  Creek  Acres 
and  proceeds  in  a  southwesterly  direction  over  Tarryall  Reservoir  toward  Antero  Junction.  The 
terrain  profile  of  this  particular  leg  is  shown  in  Figure  5.44.  The  elapsed  time  between  the  first 
turbulence  sample  gathered  during  the  dawn  flight  (Figure  5.33)  and  the  last  turbulence  sample 
gathered  during  the  mid-aftemooh  flight  (Figure  5.41)  is  approximately  seven  hours. 

GUST  VELOCITY  RMS  VALUES 


Another  method  of  investigating  turbulence  magnitude  is  to  analyze  the  standard  deviation  of  the 
gust  velocities.  The  standard  deviation  (<rt)  of  each  of  the  three  gust  velocity  components  was 
calculated  hy  using  equation  5.2. 

The  values  of  x^re  taken  from  the  time  histories  of  gust  velocities  which  are  computed  at  the  rate 
of  100  samples  per  second  giving  27,000  values  for  each  four  and  one-half  minute  turbulence 
sample.  Therefore,  for  this  particular  calculation,  N  =  27,000.  The  mean  gust  velocity  value  is 
indicated  by  x. 

Since  the  mean  value,  x ,  is  made  essentially  equal  to  zero  by  high-pass  filtering,  the  standard 
deviation  is  actually  equal  to  the  root-mean-square  (rms)  value  as  shown  in  equation  5.8.  The 
two  terms  are  used  interchangeably  in  this  report. 


a 


t 


=  rms  with  x  a  0 


(5.8) 


The  gust  velocity  rm.  .alues  are  being  analyzed  in  a  number  of  ways.  Primarily  these  involve 
analysis  of  rms  magnitude  on  a  statistical  basis,  a  study  of  the  gust  velocity  rms  values  as  a 
function  of  wind  speed,  and  investigation  of  various  turbulence  forecasting  models  (discussed  in 
Section  IV)  using  A  velocity  rms  values  and  regression  techniques.  Since  the  gust  velocity  rms 
values  are  being  investigated  within  the  realm  of  the  categorization  technique  described  in  Section 
II,  then  any  relationship  with  these  geophysical  categories  also  becomes  apparent. 
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Gust  Velocity  RMS  Statistical  Analysis 


Probabilities  of  gust  velocity  rms  values  being  a  given  value  or  greater  are  computed.  The  rms 
values  are  grouped  into  bands  O.S  fps  m  width  3nd  the  cumulative  probabilities  ca'culsted.  The 
cumulative  probabilities  determine  1  for  all  data  recorded  thus  far  during  Phase  II!  ai*  presented 
in  Figure  5.45  throi  gh  5.47.  For  comparative  purposes,  the  data  obtained  during  Phases  I  and  II 
of  the  program  are  also  presented.  As  can  be  seen,  the  Phase  lil  gust  velocity  rms  vahies  are 
greater  in  magnitude  than  the  corresponding  seasonal-location  Phases  I  and  II  data.  This  is 
attributed  to  the  fact  'hat  *he  longer  turbulence  wavelengths  are  being  measured  more  accurately. 

A  statistical  analysis  of  rms  values  by  geophysical  category  is  being  accomplished.  In  order  to  assure  a 
good  statistical  sampling,  a  minimum  of  30  data  samples  within  the  category  being  analyzed  war 
established  as  the  criteria  for  a  valid  analysis  of  this  type.  Therefore,  all  of  the  various  combinations 
of  geophysical  categories  being  used  during  this  program  cannot  be  analyzed.  There  an  not  enough 
data  for  such  ;n  extensive  investigation.  However,  those  categories  which  have  the  greatest  effect  on 
gust  velocity  rms  magnitudes  are  being  analyzed. 

The  results  analyzed  for  this  report  are  shown  in  Figures  S.48  through  5.58.  Season  of  the  year 
has  not  been  analyzed  separately  in  any  of  these  plots.  As  mentioned  previously .  it  is  redundant 
with  location  for  the  Phase  III  testing 

As  can  be  seer,  in  these  Figures,  the  results  are  similar  for  all  three  gust  velocity  components.  A 
discussion  of  general  effects  related  to  individual  category  components  is  presentee  below, 

Figure  5.48  shows  the  effects  of  terrain  over  which  the  data  were  obtained  Turbulence  over  the 
rougher  terrains  produced  the  highest  gust  velocity  rms  values.  Data  recorded  over  the  water  leg  show 
the  lowest  turbulence  magnitudes. 

The  data  were  recorded  at  the  two  nominal  test  altitudes  of  250  and  750  feet  above  the  terrain. 
Resultant  gust  velocity  rms  probabilities  are  shown  in  Figure  5.49  and,  as  can  be  seen,  the  rms  values 
are  greater  at  the  lower  altitude.  As  shown  in  Appendix  VI,  the  average  radar  altitude  did  not  always 
correspond  with  the  250  and  750  font  category  altitudes  ar  far  as  exact  magnitude  is  concerned.  For 
the  greater  part,  this  occurs  for  flight  over  the  rougher  terrains,  bmall  crevices  and  canyons  into  which 
the  airplane  cannot  maneuver  to  maintain  absolute  contour  flight  increase  the  average  radai  altitude. 

Effects  of  atmospheric  stability  on  gust  velocity  rms  values  are  shown  in  Figure  5.50.  The  lowest  rms 
values  were  recorded  when  ‘he  atmospheric  sm’  :iity  wss  very  stable.  The  cumulative  probability 
increases  with  decreasing  atmospheric  stability. 

As  the  time-cf-day  during  which  the  data  were  recorded  vanes  iiom  dawn  to  mid-morning  or  to 
mid-afternoon,  the  cumulative  probability  of  encountering  ar  mis  gust  of  a  given  magnitude 
increases  This  is  shown  in  Figure  5.5  1 
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Figure  5.52  shows  the  gust  velocity  rms  cumulative  probability  curves  fcr  data  from  each  of  the  two 
locations  associated  with  this  report. 

During  Phases  I  and  II  all  combinations  of  terrain,  altitude,  and  stability  were  analyzed.  This  involves 
40  different  category  combinations,  For  the  Phase  III  data,  however,  only  six  of  these  combinations 
contained  30  or  more  data  samples  Therefore,  for  this  interim  report  the  effects  of  terrain  and 
altitude  were  determined  by  studying  the  combinations  of  these  two  category  components. 
Cumulative  probability  plots  for  the  combinations  are  shown  in  Figures  5  53  through  5.58.  A 
summary  of  these  plots  is  shown  in  Figure  5.59  where  the  gust  velocity  rms  mean  values  are  plotted  as 
a  nmction  of  terrain  and  altitude. 

These  data  show  how  the  turbulence  level  increases  as  altitude  decreases  for  all  types  of  terrain  .  Also, 
the  turbulence  level  is  seen  to  increase  with  the  rougher  terrains.  The  one  exception  to  this  is  with 
regard  to  the  plains  and  desert  terrain.  Although  the  desert  terrain  is  rougher  than  for  the  plains, 
larger  gust  velocity  rms  magnitudes  were  recorded  over  the  plains. 

It  has  been  found  that  the  effects  of  noise  are  significant  at  the  lower  gust  velocity  rms  magnitudes. 
According  to  theory  fornormal  distributions: 

<  -  *1  ♦  (5.9) 


Where:  a  t  =  standard  deviation  of  the  turbulence 

=  standard  deviation  of  the  noise 
<r  -,  =  standard  deviation  of  the  rejorded  data 

If  #  is  assumed  constant  for  a  given  category,  this  value  becomes  a  greater  and  greater  percentage 
of#yas*tbecumes  smaller  and  smaller.  Therefore,  at  lower  turbulence  levels  the  distribution  is 
distorted  by  the  high  percentage  of  noise  content  in  the  data  where  the  noise  has  a  greater  effect  on 
the  distribution  *\andocs  the  turbulence. 

Because  of  the  low  signal-to-noise  ratio  criteria  and  the  results  from  the  noise  analysis  discussed  in 
Section  ill,  the  cumulative  probability  distribution  plots  in  Figures  5.45  through  5.58  do  not  extend 
below  a  gust  velocity  rms  value  of  1 .5  fps.  However,  so  that  no  error  will  occur  through  evaluation  of 
these  data,  the  rms  values  below  1 .5  feet  per  second  have  been  used  in  computing  the  probabilities  for 
the  data  above  1.5  fps.  Those  samples  that  were  not  processed  due  to  the  low  turbulence  levels 
have  been  assumed  to  have  gust  rms  values  below  1.5  fps  and  arc  included  in  the  cumulative 
probability  distribution  calculations  by  being  added  into  the  total  number  of  data  samples  within 
the  category  being  investigated. 

Addition  of  these  lew  turbulence  level  conditions  does,  of  course,  have  an  effect  on  th~  mean  values 
of  the  data.  These  mean  vaiues  were  corrected  by  asaunnng  that  1/3  of  the  samples  added  are 
contained  in  each  of  the  three  bands  below  1 .5  fps  (i.e.,  1/3  in  the  0  to  0.5  fps  bsnd,  1/3  in  the  0.5  to 
1 .0  fps  band,  and  1  /3  in  the  1 .0  to  1 .5  fps  band).  The  average  value  of  the  added  samples  within 
each  of  the  bands  is  then  assumed  to  be  the  middle  value  of  the  band.  The  values  plotted  in  Figure 
5.59  were  corrected  in  this  manner. 


The  distribution  shape,  or  form,  of  all  gust  velocity  rms  values  recorded  during  Phases  1  and  II  of 
the  LO-LOCAT  program  was  determined  to  be  a  sum  of  two  normal  distributions  (see  Reference 
1.2).  Distributions  of  the  data  obtained  so  far  during  the  Phase  III  program  and  those  from  Phases 
1  and  II,  are  shown  together  in  Figures  5.60,  5.61,  and  5.62  for  comparison.  In  these  figures,  the 
cumulative  probabilties  are  plotted  versus  a  standardized  variable  t,  where  t  is  computed  as 
follows: 


t 


Ou*t  Velocity  rms.p 


(5.10) 


In  this  manner,  the  effects  of  turbulence  intensity,  as  represented  by  the  mean  (m)  and  the  dispersion 
of  the  distribution  <et)  are  removed  and  the  data  are  compared  on  the  basis  of  distribution  form  only. 

The  rms  values  which  were  not  calculated  because  of  low  turbulence  levels  were  accounted  for  in 
the  calculation  of^ando^by  estimating  the  distribution  of  these  values  within  the  gust  velocity 
rms  bands  below  1.5  fps.  These  estimations  were  made  by  extrapolating  the  cumulative 
probability  curves  to  zero(Figures  5.63,  5.64,  and  5.65).The  number  of  samples  within  each  of  the 
lower  bands  can  then  be  determined.  The  mean  value  of  the  data  can  then  be  calculated  by 
assuming  that  the  mean  value  of  the  samples  within  each  of  these  lower  three  bands  (0.0  -  0.5,  0.5 
to  1 .0,  and  1 .0  to  1 .5  fps)  is  at  the  mid-point  value  of  each  respective  band.  The  dispersion  of  the 
distributions  is  corrected  by  noting  that: 


(5.11) 


(5.12) 


where  the  subscript  1  represents  the  set  of  estimated  rms  values  and  the  subscript  2  represents  the 
calculated  set  of  rms  values.  It  can  be  shown  that  the  dispersion  of  the  distribution^)  for  the 
combination  of  both  sets  of  data  is: 


^  -2  1  1"  Ni*i  +  n2  Me  i  i 

h  ••J'l  ^xa'-J  ( 


1/2 


(5.13) 


After  these  ecirections  were  made,  the  values  of  the  standardized  variable,  t,  were  calculated, 
and  Figure  5  60,  5.61,  and  5.62  were  prepared.  The  distribution  shapes  of  the  different  Phases 
match  well  at  the  higher  cumulative  probability  values  but  seem  to  vary  at  some  of  the  lower 
orobabuities  of  the  higher  gust  velocity  rms  magnitudes. 
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Gust  Velocity  RMS  Correlation  with  Measured  Phenomena 


Gmt  Velocity  RMS  Versus  Inflight  Measured  Wind  Speed 

The  effects  of  wind  velocity  on  gust  velocity  rms  values  are  being  investigated  during  Phase  111 
testing.  Because  of  the  mechanical  generation  of  turbulence  caused  by  wind  blowing  across  the 
terrain,  the  relationship  between  gust  velocity  rms  values  and  recorded  inflight  wind  velocity  has 
been  investigated  with  relation  to  type  of  terrain  and  altitude  above  the  terrain  associated  with  the 
data. 


The  gust  velocity  nns  values  were  plotted  versus  the  wind  velocity  and  the  first-oider  least  square  line 
which  best  fit  *he  data  was  computed.  The  dope  of  the  least-square  line  indicates  the  amount  of 
correlation  between  rms  gusts  and  wind  velocity.  As  the  slope  increases,  a  greater  correlation  is 
indicated.  Calculation  of  the  y-intercept  then  allows  calculation  of  the  equation  of  the  line  using  the 
point-slope  method. 


where 


*  y0  ♦  ew 

m  =  slope  of  line 


(5.14) 


y.  =  y  intercept 

The  standard  deviation  of  the  data  about  the  least-square  line  is  calculated  as  follows: 


where: 

x  =  value  of  gust  velocity  rms  at  a  given  wind  speed 
=  ',alue  of  least  square  line  at  the  given  wind  speed 


The  standard  deviation  of  the  data  about  the  line  indicates  the  degree  of  data  scatter  The  results  are 
shown  <n  Figures  5.66  through  5.73.  Figure  5.74  summarizes  the  wind-turbulence  relationship  by 
showing  the  slopes  of  the  least-square  lines  plotted  as  a  function  of  terrain  and  altitude. 

Certain  trends  apparent  in  these  data  are  as  follows: 

•  The  turbulence  level  increases  at  a  greater  rale  with  increasing  wind  velocity  at  the 
250-foot  altitude  over  high  mountains  than  for  any  of  the  other  terrain-altitude 
combinations  studied.  At  750  feet,  the  greatest  correlation  occurs  for  data  obtained  over 
low  mountains.  The  correlation  between  turbulence  and  wind  speed  is  less  over  plains 
and  desert  terrains.  At  750  feet,  over  desert  terrain  there  is  a  decrease  in  turbulence 
intensity  with  an  increase  in  wind  speed. 
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•  The  effect  of  altitude  varies  depending  upon  the  type  terrain  involved.  O'er  hi^h 
mountains,  the  correlation  between  turbulence  level  and  wind  velocity  is  considerably 
greater  at  2  SO  feet  than  at  750  feet.  Over  the  plain?,  the  effect  of  wind  speed  on  mis  values 
is  small  and  approximately  the  same  regardless  of  altitude.  The  gust  velocity  rms  values 
increase  at  a  greater  rate  with  increasing  wind  speed  for  data  obtained  at  750  feet  than  for 
250  feet  over  low  mountains  This  is  the  only  type  of  terrain  over  which  this  phenomena 
occurs.  These  trends  appear  to  be  the  same  for  all  three  components  of  turbulence. 

Gust  Velocity  RMS  Correlation  with  Surface  Wind  Measurements 


An  attempt  was  made  to  define  a  relationship  between  routine  meteorological  data  from  established 
U.SW.B.  stations  an.  the  rms  of  observed  gust  velocities.  The  following  assumptions  and 
simplifications  were  made. 

•  The  stations  were  sufficiently  close  to  the  flight  routes. 

•  The  wind  recorded  at  the  3-hour  time  period  nearest  the  flight  time  was  representative  of 
the  wind  at  flight  time 

A  gust  factor  was  computed  (similar  to  reference  5.12)  by  dividing  the  maximum  wind  speed  by 
the  average  wLnd  speed  for  one  day  as  follows' 


MAX  WIND  FOR  24  HOURS 
AVG  WIND  FOR  24  HOURS 


Daily  gust  factors  are  plotted  versus  time  i  i  Figure  5.7c  Fcr  the  McConnell  route,  the  mean  value 
is  approximately  1 .9.  with  extreme  values  as  high  as  7.C. 

The  3-hour  winds  for  each  test  were  multiplied  by  the  gust  factor  and  compared  to  the  maximum 
rms  value  recorded  as  shown  in  Figure  5.76.  It  appears  that  there  is  a  relationship  of  increasing 
gust  velocity  rms  with  increases  in  the  product  of  gust  factor  and  wind  speed. 

Gust  Velocity  PMS  Regression  Models 

Correlation  of  the  Phase  III  gust  velocity  rms  values  with  some  of  the  regression  models  from 
other  turbulence  studies  is  currently  being  investigated.  Results  of  some  of  these  studies  using  the 
LO  LOCAT  Phases  1  and  II  data  have  been  prepared  and  are  presented  herein.  These  analyses  are 
for  the  vertical  component  of  turbulence  only. 

The  regression  models  presented  in  Tables  5.2  through  5.5  have  been  tested  with  the  I  O-LOCAT 
Phases  I  and  1!  data  These  tables  show  the  average  percent  difference  between  the  actual  recorded 
gust  velocity  rms  value  and  the  value  predicted  by  the  mod"I.  They  were  d.termined  by 
computing  the  percentage  difference  between  each  sample  and  the  value  pr  dieted  by  the 
applicable  regression  model,  and  then  averaging  the  absolute  values  for  all  these  percentage 
differences 
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'  It  should  be  pointed  out  that  the  Phases  I  and  II  data  utilized  in  this  analysis  do  not  contain  any 

i  gust  velocity  rms  values  below  2.40  fps.  As  explained  in  the  discussion  of  Forecasting  Techniques 

■  ir  Section  IV,  rms  values  below  this  magnitude  have  shown  poor  correlation  with  the  forecast 

data.  Thus,  these  lower  gust  velocity  rms  values  have  been  excluded  from  the  data  included  in  t.u; 
regression  model  study  Exclusion  of  these  data  decreased  the  percentage  difference  values. 

In  Reference  S.I,  regression  models  are  presented  for  flight  over  each  cf  three  different  types  of 
j  terrain:  smooth,  low  mountains,  and  high  mountains.  In  these  equations,  gust  velocity  rms  value  is 

expressed  as  a  linear  function  of  mean  wind  speed.  The  equations  are  based  on  the  B-66B 
t  Low-Level  Gust  Study,  The  results  using  these  models  are  shown  in  Table  5.2. 
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TABLE  5.2 

REGRESSION  MODELS  BASED  ON  B-66B  AIRPLANE  LOW  LEVEL  GUST  STUDY  DATA 


Terrain  Type 

Model 

Percent 

Different 

High  Mountains 

»t,_  -  5.91  +  0. 1W 

111.038 

Lou  Mountains 

*tv "  4.2/  +  0.07  W 

61.028 

Deser c 

-  2.15  +  0.12  W 

39.016 

Plait,'5 

*t.,-  2.15  +  0.12  W 

86.174 

As  can  be  men,  the  percentage  differences  are  quite  la-ge. 

Table  5.3  presents  some  of  the  regression  models  derived  from  data  obtained  during  the  B-52G/H 
Flight  Load  History  Program  (Reference  5.2)  from  May  1963  to  May  1964.  Not  all  the  models 
were  tested  with  the  Phases  I  and  II  data  because  the  height  difference  between  the  700  and  850 
mb  levels,  which  occurs  as  a  variable  in  some  of  the  models,  was  not  measured,  In  addition,  some 
of  the  models  tested  are  not  presented  in  the  table  because  of  redundancy.  Although  more  than 
1200  hour,  of  usable  data  were  obtained,  the  B-52  program  was  limited  by  the  lack  of 
meteorological  instrumentation  and  by  the  fact  that  accelerations  were  used  to  estimate  thn  gusts. 
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TABLE  5.3 


REGRESSION  MODELS  BASED  ON  DATA  OBTAINED  DURING  B-52G/H  AIRPLANE 
FLIGHT  LOAD  HISTORY  PROGRAM 


Terrain  Type 

Model 

Percent 

Difference 

Plains 

3.258  +  0.830expf(S~- 12 . 3) 2 
L  13 

+  1.758  x  10-4  (W> 2  +  3.167 

33.091 

High  Mountains 

3.316  +  1.070  (S;e)2  -1.264  x  10'4  li 

+  7.372  x  10-3  W 

►— 

XT' 

UJ 

00 

Low  Mountains 

■  Saine  as  High  Mountains 

16.926 

High  Mountains 

*t.  “  3.345  +  0 . 735exp(~(6  12 . 3)  3 

_3  L  8  J 

+  4.875  x  10  J  W 

-1.241  x  10“4  h 

17.327 

Low  Mountains 

*tw“  Same  as  High  Mountains 

19.102 

High  Mountains 

«t  -  3.195  +  1.720  S  , 

■ 

14.570 

Low  Mountains 

■  Same  as  High  Mountains 

20.951 

The  equations  presented  are  in  terms  of  absolute  altitude  (H),  in  feet,  true  solar  time-of-day  (ST), 
in  hours,  and  sine  of  the  solar  elevation  (Sss  ). 

Table  5.4  presents  some  derived  regression  models,  baaed  on  Phases  1  and  II  data,  that  are 
analogous  to  the  Reference  5.1  models  previously  shown,  in  that  the  only  independent  variable  is 
the  mean  wind  speed. 


TABLE  5.4 

REGRESSION  MODELS  BASED  ON  PHASES  I  AND  II  DATA 


Terrain  Type 

Model 

Percent 

Difference 

High  Mountains 

17.324 

Low  Mountains 

Desert 

3.4857  +  0-004203  W 

13.257 

Plains 

*tv-  2.9269  +  0.04086  W 

15.311 

f 

I 
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The  models  of  Table  5.4  indicate  an  acceptably  low  percentage  error  but  they  should  be  used 
with  caution.  The  equations  indicate  that  gust  velocity  rms  values  can  be  considered  constant 
except  for  a  small  correction.  This  is  not  a  bad  assumption  most  of  the  time,  but  these  models  fail 
in  predicting  high  gust  values.  Since  the  regression  coefficients  are  quite  small,  very  strong  winds 
would  be  required  to  predict  even  moderate  rms  values. 

In  conjunction  with  this  reasoning,  it  must  be  pointed  out  that  wind  speed  and  terrain  are 
obviously  not  the  only  factors  which  determine  gust  velocity  magnitude.  In  Reference  i.2, 
altitude  above  the  terrain  and  atmospheric  stability  were  shown  to  be  vtry  important  influences 
on  the  magnitude  of  gust  velocity  rms  values  at  low  altitudes.  Other  factors  influencing 
turbulence  level  include  time-of-day,  season  of  the  year,  geographic  location,  wind  gradients, 
temperature  gradients,  synoptic  weather  patterns,  and  other  variables.  Thus,  any  analysis  of 
turbulence  levels  as  2  function  of  wind  only  will  not,  under  most  circumstances,  provide  the  exact 
values  of  gust  velocity  rms  recorded,  although  trends  can  be  established  from  the  relationships.  It 
should  also  be  pointed  out  that  in  Reference  1.2,  the  slopes  and  intercepts  of  the  gust  velocity 
rma-wind  linear  relationship?  for  Phases  I  and  II  data  are  presented  as  a  function  of  terrain, 
absolute  altitude,  and  atmospheric  stability.  All  of  the  gust  velocity  rms  values  recorded  are 
included  in  that  analysis,  (i.e.,  even  those  rms  values<2.4  fps). 

The  equations  in  Table  5.3  present  gust  velocity  rms  values  as  a  function  of  more  than  just  terrain  and 
wind  speed.  As  previously  discussed,  absolute  altitude  is  included  as  well  as  solar  time-of-day  and  sine 
of  solar  elevation.  These  latter  two  variables  do  relate  to  season  of  the  year  and  atmospheric  stability. 
The  relationship  between  time-of-day  and  stability  for  the  Phases  I  and  I!  data  is  presented  in 
Reference  1.2. 

The  Phases  I  and  II  data  have  also  been  analyzed  with  respect  to  the  combination  of  wind  speed  and 
sine  of  solar  elevation.  Some  preliminary  regression  models  for  vertical  gust  velocity  rms  magnitudes 
are  presented  in  Table  5.5. 


TABLE  5.5 


REGRESSION  MODELS  BASED  ON  PHASES  I  AND  II  DATA 


TT9 *rl -if 


DERIVED  EQUIVALENT  GUST  VELOCITY 

An  approximate  value  of  derived  equivalent  gust  velochy  was  calculated  for  each  turbulence 
sample.  These  values  were  used  to  establish  the  approximate  turbulence  level  encountered  and  as  an 
immediate  guideline  in  determining  if  gust  encounters  were  of  a  magnitude  to  warrant  an  airplane 
structural  inspection  Maximum  incremental  acceleration  from  a  one  g  mean  was  determined  from 
“quick-look”  strip  charts  of  vertical  acceleration  at  the  airplane  center  of  gravity.  Computation  of 
was  accomplished  using  the  equation. 


C  rt 

i.  p  o  S'/.  Kg  , 


Allies  (feet/second) 


(5.17) 


where  Kg,  the  gust  alleviation  factor,  is  expressed  as: 


.88* 

5.3+* 


(5.18) 


and  m  .  the  dimensionless  airplane  mass  ratio  is  defined  by: 

2  W 


C,  Pc g  S 

^  ■ 


(5.19) 


Using  the  equation  of  state,  p  (air  density)  may  be  expressed  as: 


^(4r)-— 4r(^) 


(5.20) 


The  values  of  U4e equation  constants  for  the  T-33A  airplane  are  defined  below" 
c  =  mean  aerodynamic  chord  =  6.72  feet 
g  =  gravitational  constant  =  32.2  feet,  sec^ 


*0  =standard  day  air  densuy  =  .0023  78  slug'feet 
S  =  wing  area  =  238  feet” 
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Constant  values  were  determined  for  the  equation  variables  (CL,  P,/T  ,  W  and  Ve)  to  facilitate 
calculations  in  the  field.  This  was  accomplished  by  using  an  average  gross  weight  fW)  for  the 
airplane  over  the  test  route,  an  equivalent  airspeed  (Ve)  approximated  by  target  test  airspeed,  and  the 
slope  of  the  lift  coefficient  versus  angle  of  attack  curve  (CL  )  for  a  Mach  number  equivalent  to  the 
target  airspeed  and  the  average  altitude.  The  ratio  of  static  pressure  to  ambient  temperature  (F^/Ta) 
was  estimated  from  data  obtained  during  tne  LO-LOCAT  Phac*s  1  and  11  Program.  Values  for 
P,/T»  were  calculated  by  averaging  the  Phases  1  and  II  data  over  a  corresponding  seasonal  time  period 
for  each  similar  route  leg.  The  values  for  these  variables  are  listed  below: 


C La  ~  5.95/radian 

W  =  12,500  lbs. 

Ve  =591.15  feet/sec, 

Fj/T**  .0534 

R 

McConnell  (all  legs) 

P.O;=  .0500 

in.  hg. 

_rR 

Edwards  (legs  1,3,5,  and  8) 

P, /I.  “.0550 

in.  hg. 
°R 

Edwards  (legs  2  and  4) 

P./V  .0528 

in.  hg. 

Edwards  (legs  6  and  7) 

Substitution  of  the  above  constants  in  Equation(5.1 7)  gives  the  folio »ving  equations  for  calculations 


ofU4#fr°m  "quick-look”  strip  charts  of  : 

U4e  =16.31Ar^ro  (McConnell,  all  legs) 

(5.21) 

Ude  =16.18AU,C0  (Edwards,  legs  1 ,3,5,  and  8) 

(5.22) 

Ude  =  16.37Anzcq  (Edwards,  legs  2  and  4) 

(5,23) 

Ude  =  16.29Aj^ca  (Edwards,  legs  6  and  7) 

(5.24) 

The  values  of  U4.are  tabula’ea  in  Appendix  VI  for  each  turbulence  sample.  The  maximum  vertical 
gust  velocity  for  each  turbulence  sample  is  also  tabulated. 


Values  of  and  maximum  vertical  gust  velocity  were  categorized  with  respect  to  terrain  and 
cumulative  probability  calcui  ted  for  each  p arameter.  Tiese  data  are  shown  in  Figures  5.77  through 
5.81  for  the  high  mountain,  low  mountain,  plains,  desert,  and  water  terrain  categories,  respectively. 
The  cumulative  probability  of  the  maximum  calculated  vertical  gust  velocity  and  for  the 
combined  McConnell  and  Edwards  routes  were  also  computed  and  the  two  cumulative  probability 
distributions  are  presented  in  Figure  5.82. 

These  distributions  indicate  that  consistently  underestimated  tne  maximum  vertical  guvt 
velocities  for  the  higher  gust  velocities  and  overestimated  them  at  the  lower  gust  velocities. 
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TURBULENCE  SPECTRA 


Power  spectra  were  calculated  for  429  out  of  985  four  and  one-half  minute  turbulence  samples. 
These  samples  were  evaluated  for  instrumentation  anomalies,  low  signal-to-noise,  high  coherency, 
and  ronhomogeneous  turbulence.  Approximately  150  samples  not  falling  in  these  categories  were 
analyzed  for  conformance  to  isotropic  theory,  comparison  to  mathematical  spectra  models,  and 
scale  lengths. 

Instrumentation  Anomalies 


An  example  of  power  spectra,  with  suspected  instrumentation  problems,  is  presented  in  Figure 
5.83.  Here  it  is  quie  obvious  that  the  spectra  are  suspect  due  to  the  large  hump  in  the  lateral 
spectra,  which  occur;  at  the  airplane  short  period  mode,  and  due  to  the  irregular  shape  of  the 
vertical  spectra.  Experience  has  shown  that  when  spectra  such  as  that  shown  in  Figure  5.83  are 
obtained,  it  can  ultimately  be  proven  to  be  caused  by  instrumentation  problems. 


Low  Turbulence  Level 

An  example  plot  of  power  spectra  where  the  turbulence  level  was  low  is  presented  in  Figure  5.84. 
Turbulence  data  samples  obtained  during  Phase  Ill  containing  low  signal-to-noise  ratios  ty  pically 
have<rtvalues  less  than  1 .5  fps,  irregularly  shaped  spectra  (i.e.,  spectra  which  do  not  have  a  definite 
-5/3  slope  and  in  which  the  power  turns  upward  at  the  higher  frequencies),  and  uncommonly  large 
integral  scale  lengths 

Coherency 


The  statistical  independence  of  the  gust  components  u,  v,  and  w,  is  evaluated  using  the  coherency 
function.  The  coherency  function  between  u  and  v,  and  v  and  w,  is  calculated,  e  g. 


,00 


[*uv00]2 

*u(k)  fvW 


(5.25) 


Two  parameters  are  said  to  be  completely  dependent  (coherent)  if  their  coherency  function  is  equal 
to  1.0  and  completely  independent  (not  coherent)  if  the  coherency  is  equal  to  zero.  Due  to 
statistical  variations  in  the  recorded  data,  electronic  noise,  etc.,  a  coherency  of  zero  is  impossible  to 
achieve.  The  gust  velocity  components  are  therefore  considered  to  be  independent  if  the  coherency 
function  approaches  zero. 

The  coherency  function  is  used  as  an  aid  in  defining  those  turbulence  samples  which  are  invalid. 
Electronic  noise  and  instrumentation  malfunctions  such  as  sensitivity  changes  will  increase 
coherency  since  the  same  measurements  are  used  in  calculating  the  three  gust  velocity  components. 
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The  average  coherency  function  is  presented  in  Figure  5.85.  Also  presented  is  a  band 
representative  of  the  scatter  in  the  data.  The  average  coherency  function  was  calculated  by 
determining  the  value  of  the  coherency  function  of  each  sample  at  spatial  frequencies  of  7  x  10'  , 
1  x  10'\  1.5  x  10*\  2  x  Ur\  3  x  10‘\  4  x  10‘\  6  x  10'\  1  x  KT3,  1.5  x  10'\  2  x  10'3,4  x 
1  O'3 ,  7  x  10’3,  1  x  10"*,  and  1.5  x  10**  cpf.  The  values  obtained  at  each  of  these  frequencies 
were  then  averaged  and  the  standard  deviation  calculated.  The  average  line  in  Figure  5.a5  was 
obtained  by  plotting  the  average  value  obtained  at  each  of  the  above  spatial  frequencies.  The 
scatter  jind  was  obtained  by  plotting  the  mean  plus  and  minus  the  standard  deviation  of  the  data 
at  each  point. 
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Homogeneity 


To  accurately  define  valid  scale  lengths  and  spectra  shape  only  turbulence  samples  considered  to  be 
homogeneous  can  be  used. 

Homogeneity  is  the  property  which  in  space  coordinates  corresponds  to  the  stationarity  property 
in  t  me  coordinates.  The  entire  analyses  of  time  histories  of  gust  velocities  as  recorded  by  m 
airplane  are  based  upon  the  validity  of  Taylor's  hypothesis.  Assuming  Taylor’s  hypothesis  to  be 
valid,  then  the  stationarity  of  the  time  history  is  indicated  by  the  spatial  homogeneity  of  tire 
turbulence  sample.  Detailed  discussions  of  Taylor's  hypothesis  may  be  found  in  References  5.3 
and5. 11 .  An  indication  of  turbulence  homogeneity  was  obtained  by  dividing  the  spectrum  of  the 
middie  third  of  a  data  sample  by  the  spectrum  of  the  entire  sample.  Complete  homogeneity  is 
indicated  when  the  specna  ratio  averages  approximately  1 .0. 

Example  spectra  and  homogeneity  plo<  'or  data  considered  to  be  nonhomogeneous  are  presented  in 
Figures  5.86  and  5  87.  Similar  plots  for  homogeneous  data  are  presented  in  Figures  5.88  and  5.89. 
Time  history  plots  of  the  r.onhomogeneoux  and  homogeneous  data  are  presented  in  Figures  5.33  and 
5.34. 

Data  showing  average  homogeneitv  for  the  turbulence  data  considered  to  be  homogeneous  is 
presented  in  Figure  5.90.  The  data  plotted  consists  of  a  mean  line  and  the  scatter  band.  The 
.iter  band  was  developed  in  the  same  manner  as  for  tne  coherency  function. 

Isotropy 

A  turbulent  field  which  is  isotropic  is  one  in  which  there  is  no  preferred  direction  and  therefore  all 
average  functions  describing  this  field  must  remain  unchanged  regardless  of  the  rotation  of  the 
coordinate  system. 

Theoretical  considerations  and  experimental  evidence  have  shown  that  most  fine  or  small  scale 
structure  of  actual  anisotropic  turbulence  is  nearly  isotropic  (local  isotropy).  The  fourand  one-half 
minute  data  samples  obtained  during  this  program  are  considered  to  be  “local”  with  respect  to  large 
scale  turbulence  and  have  been  analyzed  for  isotropic  characteristics 

According  to  References  5.3.  5.4,  and  5.11.  the  isotropic  characteristics  will  follow 
♦v(k)/<J>v  (k)  "1  and.  ♦u  (k) /♦v(k)  will  vary  from  2.0  to  0.75  with  frequency  and  scale 
length.  The  characteristics  of  the  exp  erimenial  data  were  established  by  dividing  the  unnormalized 
power  spectra  of  u  and  w  by  the  power  spectrum  of  v. 
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In  this  manner,  if  the  turbulence  is  isotropic,  the  w/v  spectra  ratio  will  be  equal  to  1 .0  and  the  u/v 
ratio  will  approach  a  form  similar  to  that  in  Figure  5.91.  The  data  in  Figure  5.91  illustrate  how 
the  u/v  ratio  vanes  with  scale  length  and  was  calculated  by  taking  the  ratio  of  the  longitudinal  to 
lateral  von  Karman  equations  for  turbulence  spectra. 

Data  showing  the  average  isotropic  characteristics  and  the  scatter  band  for  Phase  III  homogeneous 
turbulence  samples  are  presented  in  Figure  5.92 

Normalized  Spectra 

The  power  spectra  for  each  condition  were  normalized  by  dividing  the  power  magnitude  by*t2 
Since  turbulence  spectra  shift  up  and  down  in  power  with  a  respective  increase  or  decrease  in*t2, 
normalizing  the  data  in  this  manner  permits  a  comparison  of  the  spectra  shape  without  having 
shifts  due  to  amplitude  changes.  The  power  spectra  were  also  normalized  by  e^L,  that  is, the  power 
magnitude  was  divided  by  the  oM.  product.  By  normalizing  the  data  in  this  manner  and  plotting 
the  data  versus  kL,  variations  due  to  bothetand  L  are  removed. 

Power  spectra  normalized  by<rt‘and  a? L  re  presented  in  Figures  5.93  and  5  94,  respectively.  The 
data  in  Figure  5.93  represent  the  average  of  turbulence  spectra  for  the  McConnell  and  Edwards 
routes.  The  data  normalized  by*2L  in  Figure  5.94  are  based  on  an  average  of  52  turbulence 
samples  except  for  the  very  low  frequencies  where  the  number  of  samples  decrease  to  a  minimum 
of  5.  The  turbulence  spectra  agree  very  well  with  the  -5/3  logarithmic  relationship  expressed  by 
Kolmogorov’s  theory.  In  general,  the  shape  of  the  spectra  agree  very  well  with  the  von  Karman 
mathematical  spectra  except  at  the  lowest  frequencies  This  is  more  readily  seen  in  Figure  5.95 
which  is  a  plot  of  the  ratio  of  the  expcnmental  and  von  Karman  spectra.  As  can  be  seen,  the  value 
of  the  PhD  is  grea*er  for  the  experimental  data  than  for  the  von  Karman  relationship  at  low  f. 
This  difference  was  also  observed  in  the  Phases  1  and  II  data.  Only  moderate  confidence  can  be 
placed  on  the  data  at  the  lowest  kL  values  due  to  the  low  number  of  values  averaged.  It  is 
anticipated  that  flight  at  higher  speeds  which  has  been  performed  over  leg  8  of  the  Peterson  route 
will  produce  additional  data  in  the  lower  k  L  region 

Experimental/Mathematical  Spectra  Comparison 

Use  of  the  power  spectral  density  approach  in  aircraft  design  is  predicated  upon  the  availability  of 
a  mathematical  model  of  the  PSD  curve.  Many  mathematical  expressions  have  been  proposed  to 
represent  the  atmospheric  turbulence  spectra.  The  expressions  most  widely  used  in  the  aircraft 
industry  at  the  present  time  include  expressions  suggested  by  Theodore  von  Karman  (Reference 
5.5)  and  by  H  L.  Dryden  (Reference  5  6)  These  expressions  vary  according  to  the  gust  velocity 
component  being  represented  and  are  as  follows. 

von  Karman  Vertical  and  Lateral  - 


von  K  arm  an  Longitudinal  - 
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Dryden  Vertical  and  Lateral 
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Dryden  Longitudinal  - 
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(5.28) 


(5.29) 


As  suggested  in  Reference  5.7,  the  equations  may  be  represented  b;  the  fallowing  general  forms 
Vertical  and  Lateral  - 


♦v  (k)  2L[l  +  8  a2  k2  (n  ♦  l)] 

•x.  "  (l  +  l*»2aak2)r-*3/s 


Longitudinal  - 
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( l+Ur2a2k2)r'^/2 


(5.30) 


(5.31) 


where  a  is  related  to  the  scale  of  turbulence,  L  as  follows: 

L[(a-1)  ■  ] 


5.32) 


Thus,  depending  upon  the  value  of  n  chosen,  various  expressions  for  the  power  spectral  density 
may  be  derived.  If  n=l/2  these  general  equations  produce  the  Dryden  expressions  The  von 
Karman  equations  are  the  special  form  of  the  general  equations  when  n  =1/3 . 

These  expressions  are  being  compared  to  the  experimentally  determined  spectra  from  the 
LO-LOCAT  program  in  order  to  investigate  the  validity  of  using  the  von  Karman  or  Dryden 
expressions  in  design.  T  hese  comparisons  are  made  by  dividing  each  experimental  spectrum  by  the 
mathematical  expression  in  question  and  plotting  this  ratio  versus  spatial  frequency  A  ratio  of  1.0 
indicates  perfect  correlation  of  the  experimental  and  mathematically  defined  spectra.  These  plots 
are  shown  in  Appendix  IV 

Average  values  of  these  experimental  to  mathematical  ratios  have  been  computed  and  plotted. 
This  is  done  by  determining  the  value  of  the  ratios  for  each  turbulence  sample  at  certain  values  of 
spatial  frequency.  These  ratios  are  then  averaged  at  each  of  the  chosen  frequencies. 
Approximately  150  turbulence  samples  determined  to  be  homogeneous  were  used  in  calculating 
th:  averages.  The  results  are  shown  in  Figure  5.96. 
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Figure  5.96  also  shows  the  good  agreement  of  the  von  K  arm  an  expressions  and  experimentally 
determined  spectra.  This  fact  is  supported  by  examining  the  standard  deviations  of  the  ratio*  for 
the  individual  sample  about  the  average  values.  1  hese  data  are  shown  in  Table  5.6. 

TABLE  5.6 


EXPERIMENTAL /VON  KARMAN  DATA 


Spatial 

Frequency 

(cpf) 

Lateral 

Average  Std.  Dev. 

Vertical 

Average  Std.  Dev. 

0.0005 

m 

0.18 

in 

0.34 

m 

mm 

0.0007 

m&m 

0.18 

■sen 

0.34 

m i 

0.00010 

1.10 

0.14 

1.46 

0.29 

1.39 

fSfs 

0.00015 

1.24 

0. 14 

1.40 

0.24 

1.33 

0.0002 

1.15 

0.18 

1.13 

0.20 

1.08 

0.20 

0.0003 

1.05 

0.20 

0.88 

0.20 

0.88 

0.23 

0.0004 

0.99 

0.20 

0.78 

0.22 

0.81 

0.22 

0.0006 

0.92 

0.17 

0.70 

0.20 

0.79 

0.18 

0.0010 

0.94 

0.17 

0.86 

0.21 

0.88 

0.20 

0.0015 

1.00 

0.  18 

0.94 

0.17 

0.96 

0.20 

0.002 

0.99 

0.  18 

0.99 

0.19 

1.08 

0.22 

0.004 

1.06 

0.23 

1.12 

0.21 

1.09 

0.21 

0.007 

1.06 

0.19 

1.14 

0.23 

1.13 

0.24 

0.010 

1.05 

0.21 

1.17 

0.22 

1.12 

0.20 

0.015 

1.01 

0.18 

1 .  10 

0.24 

1.18 

0.29 

The  Dryden  expressions  do  not  agree  well  with  the  experimental  data,  especially  at  the  higher 
frequencies.  The  basic  disagreement  is  in  the  inertial  subrange  where  the  Dryden  expressions  have 
a  -2  logarithmic  slope  rather  than  tne  -5/3  slope  of  the  data. 

The  results  discussed  here  are  similar  to  those  obtained  during  LO-LOCAT  Phases  1  and  II.  During 
that  program,  expressions  suggested  by  U.  O  Lappe  (Reference  5.8)  and  Lumley  and  Panofsky 
(Reference  5.11)  were  investigated.  The  Lumley  and  Panofsky  expression,  which  is  for  the 
longitudinal  component  only,  agreed  well  with  the  experimental  data.  The  Lappe  expression, 
however,  showed  the  same  type  of  disagreement  as  that  discussed  above  for  the  Dr/den 
expressions.  During  the  latter  half  of  the  Phase  111  program,  the  Lumley  and  Panofsky  expression 
for  the  longitudinal,  and  an  expression  suggested  by  Busch  and  Panofsky  for  the  vertical  and 
lateral  components,  are  also  being  studied. 

Turbulence  Scale  Length 


The  atmospheric  turbulence  integral  scale  length,  L,  is  an  indication  of  the  average  eddy  sue 
associated  with  the  longitudinal  gust  velocity  component.  The  value  of  the  scale  length  must  be 
known  in  order  to  determine  the  power  spectra!  density  from  either  the  von  Karman  or  Dryden 
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mathematical  expressions  (  Equations  (5.26)  through  (5.29)  ). 


Scale  lengths  can  be  computed  from  either  the  autocorrelation  function  or  the  power  spectrum 
(Reference  5.4).  Scale  lengths  are  computed  from  power  spectra  ui  the  following  manner.  A 
standard  deviation,  .covering  only  the  high  frequency  range  of  the  spectrum  is  defined,  e.. 


‘kJ 

♦  he)  dk 

k2 


1/2 


(5.33) 


For  calculations  of  L  using  the  von  Karman  expressions,  equations  5.26  and  5.27  are  simplified  by 
noting  that,  in  the  inertial  subrange,  377.5  (Lk)“and  70.78  tLkT»  1.  The  expression 
for  <t>  (k j ,  from  either  Equation  5.26  or  5.27  is  inserted  in  Equation  5.33  and  the  integration  is 
carried  out  and  solved  for  scale  length,  yielding  Equations  5.34  and  5.35,  respectively. 


I 

f 

i 

\ 

e 
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Scale  lengths  to  be  used  in  the  Dryden  expressions  are  determined  in  the  same  manner.  Equations 
5.28  and  5.29  are  simplified  by  noting  that(2rLk)  >>  1  in  the  inertial  subrange.  The 
resulting  equations  obtained  by  inserting  these  simplified  equations  into  Equation  5.33  and 
integrating  are : 


(5  36) 


(5.37) 
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Equations  5.34  and  5.36  are  expressions  for  longitudinal  scale  length  as  computed  from  either  the 
vertical  or  lateral  gust  velocity  components  while  Equations  5.35  and  5.37  represent  longitudinal 
scale  length  as  computed  from  the  longitudinal  component  of  gust  velocity. 

Scale  lengths  may  also  be  calculated  by  using  the  area  under  the  gust  velocity  autocorrelation 
curve  (Reference  5.3)  where: 


L 


Gs 


[R(r)/R(0)]dr 


(5.38) 


The  autocovariance  function  R(r)  is  normalized  by  dividing  it  by  R(0),  wnere: 


R(0) 


4>(k)  dk 


(5.39) 


During  Phases  I  and  II  of  the  LO-LOCAT  Program,  scale  lengths  were  calculated  from  the 
autocorrelation  function  for  selected  data.  -The  results  of  this  method  were  compared  to  scale 
lengths  obtained  for  the  same  data  from  Equations  5.34  and  5.35.  The  results  agreed  very  well, 
lending  credence  to  the  assumptions  made  in  deriving  Equations  5.34  and  5.35. 

Since  the  von  Kannan  expressions  give  an  accurate  representation  of  the  experimental  spectra,  the 
scale  lengths  from  these  expressions  are  being  analyzed  statistically.  The  cumulative  probability 
distributions  of  these  data  are  shown  in  Figure  5.97. 

As  mentioned  above,  the  scale  lengths  tend  to  increase  in  magnitude  as  altitude  increases.  The 
scale  length  cumulative  probability  is  plotted  in  Figures  5.98,  5.99,  and  5.100  to  illustrate  this 
point.  As  can  be  seen  in  these  figures,  the  scale  length  computed  from  the  vertical  component  is 
more  sensitive  to  altitude  changes. 

During  Phases  I  and  II  of  the  LO-LOCAT  Program,  scale  length  was  found  to  be  a  function  of 
more  than  just  height  above  the  ground.  At  the  250-foot  altitude,  scale  length  was  found  to 
increase  with  increasing  terrain  roughness.  At  the  750-foot  altitude,  atmospheric  stability 
appeared  to  have  an  effect  on  scale  length.  As  the  atmosphere  became  more  unstable,  scale  length 
showed  a  tendency  to  increase  in  magnitude.  These  trends  were  more  evident  for  scale  length 
computed  from  the  vertical  component  than  from  the  other  two  components. 

The  Phase  III  data  are  being  investigated  for  the  same  types  of  trends.  Figure  5.101  shows 
variations  of  mean  scale  length  values  with  terrain  and  altitude.  Scale  lengths  at  the  250-foot 
altitude  tend  to  increase  with  increasing  terrain  roughness,  while  those  recorded  at  750  feet  are 
somewhat  inconsistent  with  relationship  to  terrain. 
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The  mean  values  of  scale  length  are  plotted  -s  a  function  of  atmospheric  stability  for  each  of  the 
test  altitudes  in  Figure  5.102.  No  consistent  relationships  are  apparent  in  these  particular  data. 

The  average  scale  length  values  recorded  thus  far  during  the  Phase  III  program  are  664,  572,  and 
462  feet  for  the  longitudinal,  lateral,  and  vertical  components,  respectively,  Average 
corresponding  values  for  all  data  recorded  during  Phases  I  and  II  were  387,  409,  and  367  feet.  As 
mentioned  in  the  gust  velocity  tms  discussion  in  this  section,  the  values  of  at  being  recorded 
during  Phase  III  are  greater  in  magnitude  than  the  corresponding  values  recorded  during  Phases  I 
and  II.  This  is  due  to  the  longer  wavelength  turbulence  being  measured  during  Phase  III.  These 
larger  gust  velocity  rms  values  increase  the  scale  length  values  appreciably  as  can  be  seen  in 
Equations  (5.34)  and  (5,35) 

The  scale  lengths  analyzed  in  this  report  are  for  homogeneous  turbulence  conditons  only.  This  is, 
of  course,  a  preliminary  analysis  of  these  data.  At  the  conclusion  of  Phase  III,  scale  lengths  from 
all  three  phases  of  the  LO-LOCAT  Program  wit]  be  analyzed  to  provide  values  based  upon  large 
statistical  samples. 

According  to  Reference  5.9,  the  scale  lengths  for  clear  iir  turbulene  below  an  altitude  of  2,500 
feet  for  use  in  von  Karman  equation;  are  to  be  calcualted  as: 


rv  V 


(5.40) 


lk." 


Lk.-  184H1/3 


(5.41) 


A  comparison  between  the  scale  lengths  recommended  by  these  equations  and  those  obtained  so 
far  during  the  LO-LOCAT  Phase  III  Program  are  shown  in  Table  5.7, 


TABLE  5.7 


m-mcai  scale  ms  m  m.  kashas  a?fcssigns 

TO  THOSE  RECOMMENDED  BY  MIL-)’-008785A  (USAF) 


Absolute 
Altitude  ~Ft . 

Reference  5.9 
Recommended 

Scale  Length  ~  Ft . 

Average  LO-LOCAT 
Phase  III 

Sceln  Length  ~  Ft. 

L 

L, 

L 

L 

L 

L 

K  j 

Kv 

X  , 

K  % 

2  50 

1159 

1159 

250 

657 

539 

373 

750 

1670 

1670 

750 

670 

611 

556 

■■ 
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Although  the  Dryden  expressions  do  not  fit  the  LO-LOCAT  data  as  well  as  the  von  Karman 
expressions,  the  Dryden  expressions  are  used  in  aircraft  design.  Reference  5.9  ttates  that  uae  of 
the  Dryden  expressions  is  permissable  when  it  is  not  feasible  to  use  the  von  Karman  forma. 
According  to  this  reference,  the  scales  of  turbulence  to  be  used  in  the  Dryden  equations  for  flight 
below  1 ,000  feet  are: 


L-  -  H 


Ld„*  LDi-  ;ooH 


1/3 


(5.42) 

(5.43) 


Table  5.8  shows  how  the  Phase  III  data  are  comparing  with  values  calculated  from  Equations  5.42 
and  5.43. 


TABLE  5.8 
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LO-LOCAT  SCALE  LENGTHS  FOR  DRtDEN  EXPRESSION:  COMPARED 
TO  THOSE  RECOMMENDED  BY  MIL-F-008785A  (USAF) 


ABSOLUTE 
ALTITUDE ~ FT 

REFERENCE  5.9 
RECOMMENDED 

SCALE  LENGTH -FT. 

AVERAGE  LO-LOCAT 
PHASE  III 

SCALE  LENGTH  ~  FT. 

1-  A 

L-, 

L  D„ 

Ln, 

250 

630 

630 

250 

47C 

468 

366 

750 

909 

909 

750 

479 

512 

478 

As  can  be  seen,  the  Phase  Hi  scale  lengths  included  in  this  analysis  do  not  agree  well  with  those 
recommended  in  Reference  5.9.  In  geneitl,  scale  lengths  calculated  from  Equations  5.40  and  5.41 
overestimate  the  Phase  111  experimentally  determined  scale  lengths.  The  only  exceptions  to  this 
cverestlmution  occur  for  the  longitudinal  scale  length  computed  from  the  vertical  component  for 
both  the  von  Karman  and  the  Dryden  expressions.  There  i:  some  increase  In  scale  'ength  with 
increasing  attitude,  as  predicted  by  Equations  5,40  through  5.43. 


According  to  Lappe  in  Reference  5.1,  type  of  terrain  may  be  taken  into  coniideration  by  using 
the  following  formula  for  calculating  L: 


H  Lh 


(5.44) 
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where  l^and  l^are  determined  as  follows : 


Terrain  Class 

ho 

Smooth 

2/3 

135 

Low  Mountains 

1/2 

300 

High  Mountains 

CD 

675 

Values  from  Equation  5.44  are  compared  to  Phase  III  scale  lengths  obtained  over  the  McConnell 
and  Edwards  routes  in  Table  5.9. 


TABLE  5.9 

LO-LOCAT  SCALE  LENGTHS  COMPARED  TO  THOSE 
RECOMMENDED  BY  AFFDL-TR-67-122 


Absolute 
Altitude  ~Ft 

Type  of 
Terrain 

L  FROM 
EQUATION 
5.44 
~  FT 

Average  LO-LOCAT 
Phase  III 

Scale  Length  ~  Ft 

V 

B 

LK„ 

250 

Smooth 

402 

621 

472 

282 

250 

Lo.  Mount. 

425 

674 

543 

426 

250 

Hi  Mount. 

706 

819 

728 

567 

750 

Smooth 

635 

635 

581 

550 

750 

Lo  Mount. 

675 

722 

688 

506 

_ — _ 

Hi  Mount. 

769 

742 

628 

627 

It  should  be  pointed  out  that  L,  as  calculated  in  Reference  5.1,  is  designed  for  use  in  the 
mathematical  spectra  expression  suggested  therein,  not  for  use  in  the  von  Karman  equations.  The 
Reference  5.1  mathematical  expression  was  evaluated  during  Phases  I  and  II.  Those  experimental 
data  did  not  correlate  well  with  that  expression.  However,  the  fact  that  terrain  is  considered  in  the 
calculation  of  L  is  in  harmony  with  the  results  obtained  for  the  250-foot  data  during  the 
LO-LOCAT  Program.  The  results  from  Equation  5.44  agree  better  with  the  experimentally 
determined  data  than  did  results  from  Equations  5.40  and  5.41,  where  only  altitude  is  considered. 
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VISCOUS  DISSIPATION  RATE 


Dissipation  of  kinetic  energy  present  in  the  atmosphere  plays  an  important  role  in  a  complete 
atmospheric  model.  In  the  past,  several  investigators,  using  the  dissipation  of  smoke  cloudi  (or 
puffs),  have  estimated  the  viscous  dissipation.  The  numerical  results  of  these  investigations  vary 
considerably  (Reference  1.2,  Figure  4 .75).  Following  the  method  shown  in  Reference  1.2, 
dissipation  rates  were  calculated  using  the  following  equation: 


.761  /  L*. 


(5.45) 


The  small  number  of  samples  (277)  were  not  categorized  into  stability  classes.  F  igure  5 . 1 03  shows 
the  cumulative  probability  of  encountering  a  value  equal  to  or  greater  than  s.  given  dissipation 
rate  for  each  altitude.  These  data  indicate  that  the  dissipation  rate  decreases  as  the  altitude 
increase*  from  250  feet  to  750  feet,  as  expected.  The  mean  for  these  data  is  somewhat  less  than 
that  for  Phases  1  and  II. 

MICROSCALES 

It  has  been  suggested  by  Kolmogorov  that  turbulence  spectra  should  generally  be  a  function  of 
the  kinematic  viscosity  and  the  rate  at  which  the  energy  is  dissipated.  From  a  dimensional  analysis 
approach,  a  length  can  be  formulated  as  follows  (Reference  5.10): 


(5.46) 


This  length  represents  the  eddy  size  at  which  the  viscous  dUaipation  forces  become  very  strong. 
The  value  of  n  calculated  for  Phase  III  of  LO-LOCAT  la  approximately  one  mm,  wh,ch  agrees 
well  with  other  estimates  for  dlaalpaticn  in  the  loweit  layers  of  the  atmoiphere  (Reference  5,11). 
Figure  5.104  shows  the  distribution  of  n  for  each  teat  altitude.  Increases  in  thi  scale  length  are 
suggested  to  be  a  function  of  altitude  from  this  figure. 

The  Taylor  miscroscale  is  also  a  measure  of  the  smallest  eddies  responsible  for  viscous  dissipation, 
and  is  given  by  the  following  equation: 


X  ■>  4.1+309 


[@fe)I 


(5.47) 
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This  value  L  the  intercept  of  the  parabola  drawn  to  coincide  with  the  vertex  of  the 
autocorrelation  function.  This  length  lies  oetween  the  Kolmogorov  microscale  and  the  integral 
scale  as  reported  by  Lumley  and  Panofsky  (Reference  5  1 1).  Figure  5.105  shows  the  properties 
of  the  Taylor  microscale  for  each  test  altitude. 

Stability  does  not  seem  to  influence  this  length  to  the  extent  which  was  observed  for  y  The 
microscales  being  measured  during  Phase  III  of  LO-LOCAT  are  somewhat  larger  than  those 
measured  during  Phases  I  and  II.  This  is  in  line  with  the  fact  that  Phase  III  integral  scale  lengths 
too  are  somewhat  larger.  The  larger  scale  lengths  are  also  attributed  to  the  fact  that  the  long 
turbulence  wave  lengths  are  being  measured  bettei  during  Phase  III. 
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Figure  5.1  leek  Count  Example 
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Figure  5.1+  Level  Crossing 'Count  Exanple 
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Cumulative  Statute 
Miles  Represented 


0  -  Peak  Count  246 76. 9 

□  -  Level  Crossings  24676.9 


Terrain  -  All 

Altitude  -  290  and  7S0  Feet 
Stability  -  All 
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Figure  9»P  Comparison  of  Peak  and  Level  Croaaing  Count  Distributions  - 
Longitudinal  Gust  Velocity 


RRENCES  PER  STATUTE  MILE 


Cumulative  Statute 
Miles  Represented 


Fib'ui'v. 


Comparison  of  Peak  and  Level  Crossing  Count  Distributions 
Lateral  Gust  Velocity 
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Cumulative  Statute 
Miles  Represented 
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Figure  5,l6  Altitude  Effects  on  Vertical  Gust  Velocity  Peak: 
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Cumulative  Statute 
Miles  Represented 


Figure;  5. 


£1  Comparison  of  Longitudinal  Guat  Velocity  Standard 
Deviation*  Calculated  by  Different  Technique* 
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Tifur*  5.2)*  Turbulence  Spectra  Data  Calculated  trom 
Inaeeble  Avarace  Tie#  Function 
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Figure  5.25  Ti»e  El  story  of  Goat  Velocity 
Turbulence  tocounter 
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Figure  ^.26  Time  History  of  Gust  Velocity  Components  -  Severe 
Turbulence  Encounter 
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T’lrbuleni'e  Encounter 


Figure  3 *29  Ti*e  History  of  Gust  Velocity  Components  -  Severe 
Turoulence  iii  counter 


Figure  5.30  Time  History  of  Gust  Velocity  Cc^ooentB  -  Severe 
Turbulence  to  counter 
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Figure  3.32  Ti»*e  History  of  Gust  Velocity  components  -  Severe 
Turbulence  Encounter 
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5.  V*  Time  History  of  Gust  Velocity  Components  -  Severe 
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5.33  Ttae  History  of  Gust  Velocity  Oowponeuts  -  Severe 
Turbulence  Encounter 


5.36  Time  History  of  Gust  Velocity  Components  -  Severe 


Figure  37  Time  History  of  Gust  Velocity  Components  -  Severe 
Turbulence  Encounter 


Fl«ure  5.38  Tine  History  of  Gust  Velocity  Components  -  Severe 
Turbulence  Encounter 


5.39  Time  History  of  Gust  Velocity  Components  -  Severe 
Turbulence  tn counter 
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Figure  5»**0  Tine  History  of  Gust  Velocity  Components  -  Severe 
Turbulence  Encounter 


rH 

-3“ 

• 

m 

lA 

> 

CD 

• 

i-1 

-3 

rH 

I 

3^-3 

Oj 

• 

3 

i 

i 

c 

t 

K 

j  fl 

1  3 

«b 

a 

c 

« 

3 

a 

V. 

>< 

* 

> 

« 

a. 

fc* 

c 

■H 

B 

> 

« 

CL 

%- 

K 

1 

> 

1 

• 

£i 

c 

B 

> 

1 

CD 

» 

I 

CD 

i 

O 

• 

LA 

s 

ft 

• 

eo 

-it 

h- 

• 

• 

vO 

j 

u 

• 

} 

f 

u 

< 

♦-* 

c 

V 

S 

i 

« 

i 

a 

to 

u 

< 

V 

e 

♦-> 

i, 

•J 

«-> 

3 

u 

1 

a. 

&  £  £ 

Ji  5 

OJ 

f- 

l_ 

$ 

i 

vO 

dl 

OJ 

m 

* 

* 

rH 

2i  2 

2  2 

ii  i 

i 

l 

V 

u 

« 

i 

I 

£ 

l 

t 

■o 

B 

«-» 

U 

K 

P 

y 

1 

« 

e 

& 

8 

ri 

4» 

> 

| 

3 

e  t  5  $  q  s  s~  y  q 

—»  •  uimtA  uno  •<)  •  juiwt»a  uno 

163 


Turbulence  Kr counter 


172  Ambient.  Air  Temp  (**■-)  5^.9  u  w  (fP,)  kB  .6 


Figure 5 . 42  Tine  History  of  Gust  Velocity  Components  -  Severe 
Turbulence  Encounter 
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Hgure  5-**3  Time  History  of  Gust  Velocity  Components  -  Severe 
Turbulence  encounter 
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Figure  5*UL  Terrain  Profile  cf  Lag  2  -  Peteraon  Field  Route 
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Figure  5.1+5  Longitudinal  Gust  Velocity  AMS  Onuiative 
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Figure  5.46  Gust  Velocity  RMS  Cumulative  Probability  Associated 
vith  Terrain 
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Figure  5.50  Gust  Velocity  IMS  Cumulative  Probability  Associated 
vith  Atmospheric  Stability 
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Figure  5«?1  Gust,  Velocity  RMS  Guiiulatii'e  Probability  Associated  / 

with  Time -of -Day  ‘I'l1 
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Figure  5.54  Effects  of  Terrain  or.  Cucrolative  Probability  of 
Gust  Velocity  R4S  Values  Obtained  at  750  Feet 
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5«56  Effects  of  Altitude  ever  Low  Mountains  or.  Gust 
Velocity  RMS  Cumulative  Probability 
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Figure  5.59  Variation  of  Dust  Velocity  RMS  Mean  Values 
■with  Teii'-in  and  Altitude 
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igure  5*73  Gust  Velocity  RMS  Values  versus  Wind  Speed  at 
75C  Feet  Over  Plains 


awn  of  oust  moc 


$  -  Flight  Days 


LATERAL  i  GUST  VELOCITY  RMS  -  fps 

LrCOUNGLL 
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Surface  Winds 
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Cumulative  Probability  of  Ude  Figure  b.78  Cumulative  Probability  of 
finri  Haiiu  Vertical  Gust  and  Maximum  Vertical  Gust 
Velocity  -  High  ^fountains  Velocity  -  Low  Mountains 
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Figure  5-95  Average  Quit  Velocity  Coherency  Cheracterietlce 
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Figure  5.89 
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Figure  5.101  Variation  of  von  Kaiman  Longitudinal  Scale  Length 
vith  Altitude  and  Terrain 
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Figure  5.103  Viscous  Dissipation  Rat®  Statistics  for  750  Foot  oM 
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SECTION  VI 


BREN  TOWER  FLYBY 

The  test  airplane  was  flown  by  AEC's  BREN  tower  near  Jackass  Flats,  Nevada  (Figure  6.1)  on  8 
January  1969.  The  purpose  of  this  flyby  was  to  compare  turbulence  data  obtained  from  the 
tower  with  the  turbulence  data  obtained  from  the  airplane.  BREN  tower  is  located  in  a  restricted 
area  identified  as  area  number  R-4808.  The  tower  is  1520  feet  high  and  is  instrumented  to  record 
temperature  and  wind  at  various  levels.  At  the  time  of  this  flyby,  only  the  565  foot  level 
instrumentation  was  operational.  The  terrain  surrounding  the  tower  is  rolling  desert  and  hills,  with 
mountains  located  on  the  east  side.  The  Shoshone  mountains  are  directly  north  of  the  tower. 

Eight  passes,  each  approximately  nine  miles  in  length,  were  made  on  alternate  60  degree  and  240 
degree  magnetic  headings.  The  wind  was  reported  to  be  northerly  at  27  fps.  The  pilot  reported 
that  the  turbulence  was  light  during  the  passes.  Subsequent  data  processing  has  revealed  a 
maximum  gust  velocity  during  the  test  to  be  approximately  31  fps,  with  the  average  of  the 
maximum  values  for  all  8  passes  being  less  than  15  fps.  Table  6  1  lists  the  average  meteorological 
values  for  two  of  tne  passes. 

Atmospheric  survey  data  were  obtained  prior  tc  each  flyby  flown  on  a  60  degree  heading.  These 
surveys  were  conducted  by  flying  the  airplane  over  a  fixed  point  at  100  feet  and  1000  feet  above 
the  terrain.  Survey  dat  were  used  to  determine  the  vertical  gradients  which  were  then  employed 
to  describe  the  state  of  the  atmosphere  for  the  following  flyby. 

The  first  and  seventh  pas  'i  were  choaen  for  spectral  analysis  because  of  the  good  homogeneity 
and  isotropy  characteristics  of  the  turbulence.  Time  histories  of  the  gust  velocity  components  for 
each  of  ,heee  flybys  are  shown  in  Figures  6.2  and  6.3. 

The  power  spectral  densities  are  shown  in  Figures  6.4  and  6.5.  The  von  Karman  scale  lengths  were 
approximately  600  feet  The  homogeneity,  isotropy,  and  coherency  characteristics  are  also  shown 
in  Figures  6.4  and  6.5. 

The  flybys  selected  for  PSD  analysis  are  representative  low-level  turbulence  tamples  'jid  should 
provide  a  good  bails  for  evaluation  of  the  tower  data  when  they  become  available. 

TABLE  6.1 
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SECTION  VII 


SUMMARY 


Conclusions  regarding  low-level  atmospheric  turbulence  would  be  untimely  at  this  interim  point  in 
the  LC-lOCAT  Phase  III  program  A  number  of  observations  are  made,  however,  to  summarize 
the  accomplishments  and  results  achieved  to  date. 

The  most  distinguishing  characteristic  of  the  Phase  III  data  when  compared  to  data  obtained 
during  Phases  I  and  II  is  the  incidence,  on  the  average,  of  higher  rms  values,  longer  scales  of 
turbulence,  longer  Trylor  microscales,  higher  gust  velocity  magnitudes,  and  lower  viscous 
dissipation  rates.  This  is  due  to  the  fact  that  longer  wavelength  turbulence  is  being  sensed  during 
Phase  III  Long  wavelength  (14,000  feet)  turbulence  being  recorded  during  Phase  III  at  0.04  cps 
would  have  appeared  at  0.02  cps  during  Phases  I  and  II,  since  the  C131B  airplanes  were  flown  at 
one-half  the  speed  of  the  T-33A.  This  0.02  cps  data  was  below  the  useful  range  of  the 
instrumentation  causing  gust  velocity  rms  values  to  be  reduced  and  hence  the  scale  length, 
Taylor's  microscale,  etc.  Phase  III  data  are  illustrating  the  fact  that  low  altitude  turbulence 
wavelengths  longer  than  7000  feet  are  significant.  In  addition,  the  following  observations  were 
notea: 

4  Comparison  of  Phase  111  time  series  statistical  distributions  with  those  from  Phases  1 
and  II  shows  that  the  Phase  111  attributions  tend  to  increase  or  level  off  at  the  higher 
gust  velocities.  This  is  also  attributed  to  the  use  of  the  higher  speed  Phase  III  test 
airplane  which  permits  measurement  of  longer  wavelength  turbulence. 

•  Increased  gust  velocity  rms  values  are  associated  with  decreased  altitude  and  increased 
terrain  roughness.  The  only  exception  occurs  for  data  recorded  over  the  plains.  These 
data  exhibit  rms  values  of  greater  magnitude  than  those  obtained  over  desert,  even 
though  the  desert  terrain  is  rougher  than  tiie  plains. 

•  Surface  winds  recorded  near  the  flight  routes  can  be  related  to  the  rms  gust  velocity  by 
a  scale  factor  based  on  the  ratio  of  maximum  to  average  wind  speeds. 

•  There  Is  greater  correlation  between  gust  velocity  rms  values  and  recorded  inflight  wind 
speed  over  the  rougher  types  of  terrain. 

•  The  experimental  spectra  show  good  correlation  with  spectra  computed  from  the  von 
Kerman  expressions  It  was  noted  that  experimental  spectra  estimates  tended  to  be 
jomewhat  higher  thin  vow  Karnar.  estimates  a'  very  low  values  of  k L  The  Dryden 
expressions  do  not  give  as  good  a  representation  of  the  data.  This  Is  particularly  true  in 
the  inertial  subrange  where  these  eypreuaions  have  a  -2  logarithmic  slope  rather  than  the 
-5/3  slope  of  the  dc.tr.. 

•  The  avrruge  longitudinal  integral  scale  lengin  of  turbulence,  bs  calculated  from  the  vor. 
Karman  expresaions.  is  approxL*. lately  5''0  feet,  beak  length*  are  higher  for  Phase  III 
data  than  for  data  recorded  during  Phases  I  and  II 
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•  The  probability  of  increased  gust  velocity  peak  values  with  increasing  terrain  roughness, 
decreased  altitude,  and  decreased  atmospheric  stability  (with  the  exception  of  the 
unstable  category  )  is  indicated  by  the  data. 

•  Gust  velocity  peak  count  probability  distributions  were  found  to  compare  reasonably 
well  with  corresponding  amplitude  count  and  level  crossing  probability  distributions. 

•  Calculated  values  of  derived  equivalent  vertical  gust  velocity  (U,„  )  from  the  T-33.A 
airplane,  overestimate  the  maximum  measured  true  vertical  gust  velocity  component  at 
the  low  true  gust  velocity  values  and  underestimate  at  the  higher  true  gust  velocity- 
values 

•  The  maximum  discrete  gust  velocity  encountered  as  of  15  February  1969  was  74.5  fps 
for  the  vertical  component  and  occurred  over  high  mountainous  terrain  on  leg  3  of  the 
Edwards  route  The  maximum  gust  velocity  rms  was  13.1  fps  for  the  lateral  com¬ 
ponent  and  occurred  on  leg  6  of  the  Peterson  route 

•  Characteristic  frequencies  calculated  from  both. peak  count  and  power  spectra  data 
compared  favorably. 

t  As  in  the  past,  drift  existed  in  all  computed,  unfiltered  gust  velocities.  The  drift  was  not 
characterized  by  constant  amplitude  or  frequency,  but  varied  in  amplitude  and 
frequency  from  sample  to  sample 

0  Approximately  86  percent  of  the  4-1  2  minute  turbulence  samples  were  determined  to 
be  normally  distributed  based  on  gust  velocity  amplitude  statistics. 

0  LO-LOCAT  Phase  III  signal-to-noise  ratio  is  at  an  acceptable  level  when  >  1.4 

fps  Standard  deviations  of  gust  velocity  (data  +  noise  computed  from  data  obtained  in 
smooth  air)  as  low  as  .28  fps  were  obtained 

0  Probable  errors  in  the  calculation  of  LO-LOCAT  gust  velocities  are  t  2  fps  with 
maximum  ercors  of  ±4  fps 

0  Winds  recorded  at  McConnell  during  Phase  111  compare  well  with  those  recorded  at 
McConnell  during  Phases  I  and  II  Winds  recorded  at  Edwards  are  slightly  less  for  Phase 
III  than  for  Phases  I  and  I!.  The  highest  winds  recorded  so  far  were  obtained  over  the 
plains  legs  at  McConnell 

0  Climatalogieal  variance  during  the  first  two  flight  periods  (at  the  McConnell  and 
Edwards  routes)  was  approximately  normal,  with  the  exception  of  higher  than  normal 
precipitation  amounts 

0  Viscous  dissipation  rates  calculated  for  Phase  III  of  the  LO-LOCAT  program  are  lower 
than  obtained  during  Phases  I  and  11  A  comparison  of  Taylor  microscales  indicates  an 
increase,  but  near  the  expected  value  Kolmogorov  microscales  are  approximately  one 
millimeter 
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•  Data  recorded  and  processed  from  the  BREN  tower  flyby  provides  an  opportunity  for 
evaluation  of  data  from  this  tower  based  on  airplane  measured  gust  velocities. 
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Colorado  route.  Approximately  72  hours  of  low-level  (0-1000  feet)  turbulence  and  associated 
meteorological  data  recorded  from  August  16,  1968,  through  January  8,  1969,  are  presentecl.  The 
turbulence  environment  is  analyzed  in  terms  of  gust  velocity  primary  peaks,  amplitude  samples,  rms 
values,  maximum  values  and  derived  equivalent  gusts.  Mathematical  expressions  for  turbulence  spectra 
and  scale  length  statistics,  and  correlations  between  atmospheric  gust  velocities  and  meteorological  and 
geophysical  phenomena  are  shown.  The  most  predominant  characteristics  of  these  data  are  the  increases 
in  rms  values,  scale  lengths,  Taylor  micro  scales  and  maximum  gust  values  as  compared  to  the 
LO-LOCAT  Phases  1  and  II  data.  These  differences  are  attributed  to  the  increased  speed  of  the  Phase  III 
airplane  which  is  providing  a  better  definition  of  the  longer  wavelengths  of  atmospheric  turbulence.  This 
report  consists  of  two  volumes.  Volume  I  contains  a  discussion  of  data  acquisition  and  quality,  along 
with  a  preliminary  analysis  of  turbulence  and  meteorological  datailnstrumentation  details,  calibrations 
and  checks,  data  processing,  gust  velocity  data  (power  spectra,  peaViSml  amplitude  count,  etc.)  and  test 
logs  are  presented  in  Volume  11. 
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